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THE EXCITATION OF THE ENHANCED SPECTRA OF 
SODIUM AND POTASSIUM IN A LOW 
VOLTAGE ARC" 

By PAUL D. FOOTE, W. F. MEGGERS, anp F. L. MOHLER 
ABSTRACT 


Variation of the spectra of Na and K with energy of electronic impact.—The quantum 
theory, as developed by Bohr and Sommerfeld, leads to the conclusions: that all the 
arc lines are emitted by neutral atoms and should be excited at voltages above the 
ionization voltages, 5.1 for Na and 4.3 for K, and that the enhanced lines are emitted 
by singly ionized atoms and, in the case of the alkalies, should be excited when (6) the 
valence electron is knocked out of its orbit, after (a) an electron from the outer ring 
(L-ring for Na, M-ring for K) has been removed. In the case of Na, the ionization 
(a) requires 35 volts if it is complete, giving La-radiation, or 30 volts if the electron 
goes to the La orbit; and since the displacement (b) requires a maximum of 14 volts, 
the enhanced spectrum for Na should be excited by two successive impacts at 35 
volts and possibly at 30 volts, and the highest frequency should correspond to about 
14 volts. The corresponding figures for K are: 23 volts for the excitation of the 
Ma-radiation, possibly 20 volts for Ma and the enhanced spectrum, and 11 volts 
for the energy corresponding to the highest frequency. Since these conclusions 
contradict the general belief that high potentials are necessary to excite these enhanced 
spectra, a series of spectrograms was made, with a large quartz spectrograph, of the 
light excited by low voltage electrons in a special cylindrical tube in which the distance 
between the central tungsten filament, which supplied the electrons, and the accelerat- 
ing grid was so short and the vapor pressure used was so low (0.1 to 0.2 mm of Hg) 
that practically all the electrons entered the force-free space surrounding the grid, 
with the same energy. The voltages used varied from 3.5 to 5000 for Na, and from 3.5 
to 40 for K. The results completely verify the theoretical conclusions. Tables and 
photographs are given which clearly show the three-stage development in the spectra 
of these elements: (1) Single pair stage: 2.1 to 5.1 volts (Na); 1.6 to 4.3 volts (K). 
(2) Arc spectrum stage: 5.1 to 30 volts (Na); 4.3 to 20 volts (K). (3) Enhanced 
spectrum stage: above 30 volts (Na); above 20 volts (K). With a heavy ‘current 
(1 ampere) at 7 volts in K, the pair \ 4642 becomes very prominent, though it has the 


t Published by permission of S. W. Stratton, Director, Bureau of Standards. 
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series notation (1s—3d). Its presence cannot be due to incipient Stark effect, and 
accordingly furnishes an exception to the principle of selection as applied by Sommer- 
feld. The wave-lengths in I.A. and the series notations, when known, are given for 
250 lines. 

The quantum theory of spectra attributes enhanced lines to 
ionized atoms. In the case of simply ionized helium the Bohr 
theory, as generalized by Sommerfeld, not only gives the correct 
wave-lengths of the observed enhanced lines, but their fine struc- 
ture as well. 

In applying the quantum theory to heavier atoms certain 
simplifying assumptions must be made in order to carry out the 
mathematical analysis. Thus Sommerfeld’ considers the case 
where the atom has one electron in an outer orbit and where as an 
approximation the other electrons may be treated as though lying 
in a circular coplanar orbit. ‘These electrons, as far as their action 
on the outer orbit is concerned, may be assumed equivalent to a 
ring of electricity of equal total charge. 

Evaluating the kinetic energy of the outer electron, and apply- 
ing the integral [ pdg=nh both in respect to azimuth and radius as 
independent variables, in accordance with the quantum theory, 
Sommerfeld derives the following expression for the variable term 
of a spectral series, identical in form with the familiar Ritz equation: 


N(E/e)? 
(2g-+n,+a—a(m,a)}? 


(m,a) = (1) 
where E is the excess in charge of the core over that of the ring, 
n, the azimuthal and nu, the radial quantum numbers characterizing 
the outer orbits, and a and a constants which are approximately 
proportional to E and are functions of m, but not of m,. 

An atom of an alkali metal has one electron in its outer orbit 
while that of an alkali earth has two, accounting for the monovalent 
and bivalent properties respectively For the neutral atom of an 
alkali metal (£/e)?=1, and equation (1) represents the variable 
terms of the arc spectra. If one of the outer electrons of an atom 
of an alkali earth is permanently removed, that is, if the atom 
is simply ionized, equation (1) likewise applies except that 
(E/e)? =4 and the values of a and a are altered in a definite manner. 


t Atombau (2d ed.), p. 506. 
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Hence if the enhanced spectrum of magnesium is due to the simply 
ionized atom, it should resemble the arc spectrum of sodium. 
Such is the case, even to the appearance approximately of 4. 
and a*=1.5 a instead of N and a in the spectral formulae represent- 
ing the enhanced lines.'. A similar relation exists between the 
spectra of the other elements of these two groups. From the above 
and further considerations Kossel and Sommerfeld? concluded 
that in so far as enhanced spectra arise in simply ionized atoms 
the enhanced spectrum of any element should resemble the arc 
spectrum of the element of next lower atomic number. Very 
likely higher types of enhanced spectra may be produced by atoms 
which have suffered multiple ionization. 

Now, as is well recognized, the complete arc spectrum of any 
monatomic vapor may be excited by electronic collision with a 
normal atom when the energy of the impacting electron is eV, 
where V is the ionization potential. Furthermore for all mon- 
atomic vapors not exhibiting metastable modifications’ the voltage 
V is related to the highest convergence frequency v of the arc 
spectrum by the quantum condition chy=eV. The highest fre- 
quency in the arc spectra of the alkali earths is 1S and of the 
alkali metals 1s, corresponding to ionization potentials from 4 to 
10 volts. Similarly in the enhanced spectrum of an alkali earth, 
the highest frequency is 1G, and accordingly at the corresponding 
voltage V* the complete simple enhanced spectrum should be 
excited. We have verified this experimentally for magnesium.‘ 

* Sommerfeld neglected to consider the shrinkage of the ring when the number 
of electrons in the ring is maintained constant while the charge on the core is increased 
by one unit. If this is taken into account we find for sodium and potassium a*=1.5 a 
instead of a*=2 a given by Sommerfeld. The former relation is well satisfied if we 
assign values of 1, 2, 3, etc., instead of 1.5, 2.5, etc., to min the (ms) terms. Cf. Journal 
of the Optical Society of America and Review of Scientific Instruments, January, 1922. 

2 Verh. d. Phys. Ges., 21, 240, 1919. 

} From spectroscopic and other considerations there appear to be two modifica- 
tions of helium, possibly one in which the two electrons are in crossed orbits and the 


other in which the orbits are coplanar. In polyatomic vapors the ionization potential 
may involve the work of dissociation in addition to the ionization of one or more atoms. 

4The enhanced spectrum of magnesium is produced by electronic impact of 
14.97 volts corresponding to the wave-number 1@= 121,270, in which case the impact 
occurs with a simply ionized atom, with the result that the second valence electron 
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These voltages again are all small, of the order of 10 to 20 volts for 
the elements in this group of the periodic table. 

The simple enhanced spectra of the alkali metals should resemble 
the complicated arc spectra of the rare gases. Since only recently 
attempts to correlate the latter in series' have been successful it is 
not surprising that nothing is known of series relations in the 
enhanced spectra of the alkali metals. 

It is possible, however, to estimate roughly the voltages at 
which the enhanced spectra of sodium and potassium should be 
excited. In the case of sodium we have the A-ring with 2 electrons, 
the L-ring with 8 electrons, and the outer ring with the valence 
electron. There are two direct ways in which the atom may be 
ionized. First, the valence electron may be removed requiring 
an amount of work equivalent to 5.1 volts;? then a second col- 
lision of XY volts may expel an electron from the L-ring. Secondly, 
an electron may be ejected from the L-ring requiring 35 volts,’ 
followed by a collision of Y volts expelling the valence electron. 
X must be somewhat greater than 35, but, since the orbit of the 
valence electron is much larger than that of the L-ring, it is not 
evident that the repulsive force of the outer electron can materially 
assist in the expulsion of an electron from the L-ring. Hence to 
doubly ionize sodium by the first method should require electronic 
impacts of maximum value La+é6La volts where 6La is a small 
quantity and 6 is a factor. ‘The total work of double ionization 
is then (La+6La+5) volts. In the second method of ionization 
accordingly Y =é6La+5 which must be less than 35 volts. Hence 
by the first method of excitation, the enhanced spectrum of sodium 
should appear at 35(1+6) volts and by the second method at 35 
volts, and in both cases should accompany its L-radiation. 
is removed. This paper, Phil. Mag., 42, 1002, December, 1921, is of especial interest 
since De Gramont and Hemsalech, C.R., 173, 505, 1921, have just published a report of 


experiments in which they found that a potential gradient of 500 volts//cm was required 
to excite the enhanced spectrum of magnesium. 

* For example, neon: Paschen, Annalen der Physik, 60, 405, 1910. 

2 Tate and Foote, Phil. Mag., 36, 64, 1918. 

3 Mohler and Foote, “‘Soft X-Rays from Arcs in Vapors,”’ Journal of the Optical 


Society of America, 5, 328, 1921. For sodium La=35 volts. 
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Similar considerations for potassium show that its enhanced 
spectrum should be excited by direct impact at about 20 volts 
and should accompany its M-radiation." 

It is possible to arrive at an approximate estimate of the value 
of 6 in the foregoing formulae. After expulsion of an electron 
from the L-ring, the ring shrinks, and provided the valence electron 
does not fall into the L-ring—a condition which would curtail the 
emission of the enhanced lines—it assumes a new orbit where the 
angular momentum is unity. This orbit represents a higher energy 
level. Following Sommerfeld’s method of derivation except that 
two coplanar rings of electrons are assumed and their shrinkage is 
given consideration, we find for sodium and potassium: 


gta 2+) — © 
r+y ~  (g—S;)? 
: (2) 
int. sank 
(Z—sz-9)? (9—5s)? 


where a* and a are the values of a in the Ritz formula, equation 
(1), for the enhanced and arc lines respectively, Z the atomic 
number and s, the nuclear defect of the ring, i.e., effect of the 
interaction of the electrons. To an approximation, omitting 
terms beyond a in the Ritz formula: 


_ 13-55 
""G-+e) (3) 
re 54 
; ~ (1-+a*)? (4) 


where V is the simple ionization potential and V* the voltage 
corresponding to the work required to remove the valence electron 
from its new quantized orbit, after the removal of an electron from 
the L-ring. V* accordingly corresponds to the highest convergence 
frequency in the spark spectrum. 

Equation (2) gives a*=1.5 a for sodium and potassium. Com- 
puting a from equation (3) and the known ionization potentials, 
we obtain from (4) the following: 

V*=14 for sodium 
V*=11.5 for potassium. 


* Mohler and Foote, Joc. cit., observed two M-limits for potassium one at 20 volts 
and one at 23 volts. 
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The total work of double ionization is accordingly: 


Total work = La+V*=35+14=49 volts for sodium 
= 23+11=34 volts for potassium. 
Hence for sodium 


La+é6La+ V =409 


sihe= © 
b= 0.3. 
Similarly for potassium 
o= 0.2. 


These values of 6 confirm the estimates made above from more 
qualitative considerations. It is of interest to note that in sodium, 
for example, although it requires but 14 volts to excite the enhanced 
spectrum, this spectrum cannot be produced at this voltage since 
as a prerequisite the displacement of an electron from the L-ring 
is essential, which necessitates electronic impact of higher energy 
value. This is in contradistinction to the phenomena observed 
in the case of the alkali earths. Here the removal of the valence 
electron at the simple ionization potential leaves the atom in the 
fundamental state for the emission of enhanced lines (just as the 
normal, unexcited atom is in the fundamental state for the emission 
of arc lines) and further impacts of even less energy value are 
sufficient to excite the “single-line’’ enhanced spectrum 15—28 
However, with the alkalis, the simple ionization does not leave 
the atom in the fundamental state for the emission of enhanced 
lines, but rather in the state for the emission and even absorption 
of La and other X-ray lines. 

In the X-ray spectra of elements of higher atomic number it is 
known that Aa cannot be produced except when the limit Aa is 
reached, under which condition all A-lines appear simultaneously. 
If a similar condition held for the L-radiation of sodium, it would 
be impossible to excite the enhanced spectrum below the voltage 
corresponding to the limit La=35 volts. However, although a 
single-line X-ray spectrum has never been observed, it must be a 
possible phenomenon for elements of low atomic number where 
the ring concerned is the outermost “‘saturated”’ orbit. Hence 
La=35—V =30 volts might be produced at 30 volts instead of at 
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La=35. Such being the case the enhanced spectrum could be 
excited by successive impact at the voltage corresponding to La, 
the first impact (or equivalent absorption of radiation) displacing 
the electron to the La orbit, where a second impact is more than 
sufficient to doubly ionize. Still further, Mohler and Foote 
observed a new series of critical 

VARIABLE FIELD 


potentials of lower value than 7 4 


La, for several of the elements 
of low atomic number, which 
may be related to the L-series. 
If so, the enhanced spectra of 


7O SPECTROGRAPH 


sodium might be excited by 





successive impact in an arc at 
the voltage corresponding to 
this new limit. 
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heretofore all investigators 
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requisite an exciting field of 
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many thousand volts for the 








production of the enhanced ii 
spectra of sodium and potas-_ fj) 
sium. the low values above 





estimated necessitated experi- 














mental verification. 





The form of the discharge _ 

tube employed is illustrated win i “a cali fe Aare 
by Figure 1. Its novelty con- pressed field. 

sists in the use of a grid (helical 

coil) mounted extremely close to a tungsten hot wire cathode (or in 
some cases an equipotential surface heated inside) and in metallic 
contact with a concentric hollow cyclinder of relatively large diame- 
ter. By properly regulating the temperature of the sodium or 
potassium vapor, a vapor pressure may be maintained for which 
relatively few electronic-atomic collisions occur over the short 
accelerating field between the cathode and grid, and the majority 
of electrons falling into the large force-free space between the grid 
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and plate have, at the instant of collision, a velocity corresponding 
to the impressed field. The gas pressure registered by the gage 
was less than 0.001 mm of mercury and the vapor pressure was 
of the order 0.1 to. o.2 mm of mercury. 

The arcs at various exciting voltages were photographed with 
a large Hilger Type C quartz-spectrograph. Plate II shows the 
results for sodium and Plate III for potassium. Ordinary Seed 30 
plates were used for all exposures except the 3.5 volt potassium arc 
in which case the plate was made red sensitive by bathing in 
dicyanin. The pair \ 7699-7664, although present in the spectra 
at the higher voltages, was not recorded by the “ blue-sensitive”’ 
plates used for the other spectrograms. ‘The lines on all the plates 
were carefully measured and reduced to wave-lengths by the 
method of the Hartmann-Cornu interpolation formula, using lines 
of a mercury comparison spectrum and known arc and spark lines 
as normals. 

SODIUM 

Plate II shows the three successive stages in the sodium spec- 
trum. ‘The single-line spectrum consisting of the pair D, and D, 
is excited above the resonance potential, 2.1 volts, and below the 
ionization potential 5.1 volts. From 5.1 volts to about 30 volts 
the arc lines appear alone. At 30 volts numerous spark lines 
begin to appear while above 4o volts the complete enhanced 
spectrum is excited. 

The lower spectrogram was obtained by connecting a small 
5000-volt transformer between the hot wire and grid. The vapor 
was too rare for a discharge to pass except when the cathode was 
heated to give a copious emission of electrons. Since the effective 
resistance of the discharge gap must have been small compared 
to the impedance of the transformer, the actual applied voltage was 
undoubtedly much less than 5000 volts. Very few enhanced 
lines not excited at 50 volts were observed in this source but their 
intensity was considerably greater. All of the spectrograms of 
sodium were made with a very narrow slit and much of the detail 
readily apparent on the negatives is obscured in the reproduction. 
On this account reference must be made to Table I for a more 
careful consideration of the excitation of the enhanced lines, and 
verification of the foregoing statements. 
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TABLE I 


SPECTRA OF SODIUM 
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TABLE I—Continued 


‘ 
Wave-Length Volts | Volt Volts Volts Volt "te 
78 51 I- I 
75.20 5 3 2 I 
77.52 I I 
79.92 7 4 2 I 
2954.44 0 4 2 I 
3007.71 2 2 I 
I5.o1 3 y I I 
29.70 3 3 I I 
37-29 4 3 I I 
50.45 3 I I 
53 Q4 4 4 I af 
59.34 7 5 3 2 
O11. 75 I- 
04 S > I I- 
00 74 > I 
74.03 4 3 I I i" 
7o.51 5 4 2 
$0.55 2 
Q2.QI I ’e) 0 2 
3104.0 , , 
07.0 I I 
24.603 I 
29.57 5 + 
35.05 j 3 I 
396.11 I I 
40.07 I 
49.49 t 3 I 
04.11 
75-390 I I 
79.19 
* 
3189.94 a 
3212.42 3 2 
10.5 I 
20.1 
35.1 ; I 
51.3 I I 
55.30 7 j I 
00.4 I I 
74-27 3 
85.76 8 ” } 2 q 
3302. 33 3 4 3 4 fe) I 3p 
02.93 } 2 3 3 re) I 3 Pa 
On: 9... peal 2 : 
18.2 3 
20.7 “a I— I ‘ 
! 
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TABLE I—Continued 


Wave-Lensth | 00° | 100 50 30 TO | Notation for 
ave-Leng Volts | Volts | volts | Volts | Volts | — Arc Lines 
baie | : 
3327-9 ° 5 - | 
2 | one — 
3400. 2 I I Boxe , 
27.1 I | I | I | 4 7 1s—3d 
3402.55 2 2 I | : 
3533.08 I2 8 5 3 
3031.31 5 | 4 : I 
34-3 I x | 
3711.15 4 | 3 I 
| 
3745.0 [— Soom | 
1 — 


24.0 acne 
4141.0 I—,d 2p—18d 
48.0 1—,d | 2p—17d 
56.0 I—,d 2p —16d 
65.0 I— | 2p.—15d 
67.8 I— | 2pi:—15d 
77.4 I— | 2p2.—14d 
859.2 — I— 1 2p:—14d 
g2.¢ I— 2p,—13d 
4195.5 I I I 2pi1—13d 
4212.4 I | 2p.—12d 
15.4 I— I— I 2— 2p,—12d 
22.8 1,d 2p —115 
38.¢ I— I } 2 2— 11d 
41.6 I— I I I . 2pi1—11d 
49.4 [- 2p2a— 1058 
52.4 I I [= I= 2p1— 108 
3.2 I— | 2 2p.—10d 
4270.3 I I I I Dom 2p1— 10d 
07.5 I- I- I— I- 2pza— 9 
4290.6 I I I I } 1 | 2Pr— 9S 
4321.1 I I I I 3 2p.a— 9d 
24.3 2 2 2 2 4 2pi:— od 
41.3 | I I I I 2pa— 
44.5 : I I I I 2 | 2pr— oS 
go.14 3 2 2 2 4 2p.— od 
4303.45 4 2 2 3 6 2pi— 8d 
4419.04 [— I I I 2— 2pa— 75 
23.31 I 2 2— 2- 3 2p:— 75 
94. 266 5 4 4 4 5 | 2p,— 7d 
907-724 ‘A 6 5 5 5 7 | 2p:— 7d 
4541.671 ; 2 I I I 2 2p.— 6 
4545.215 3 2 2 | 2 4 2pr— OS 
46064.858 —e 7 4 4 4 7 2p.— 6d 


4005 507 aa re) 6 | ‘| | 2pi— 6d 
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TABLE I—Continued 
—— ———————— | — — oe = $$ me —— ——— 
eo 5000 } 100 50 30 10 Notation for 
Wave-Length | Volts | Volts | Volts | Volts Volts Arc Lines 
4748.016..... | 4 2 2 | 2 3 2pr— 5s 
4751.81. 5 2 2 | 2 4 2pi— 558 
4978.608....... 9 6 6 | 9 2pa— 5d 
4982.864....... 10 8 8 9 IO | 2p:— 5d 
§I49.090....... 2 2 2 2 5 | 2p2— 48 
5153.645.... 3 3 3 3 6 2pr— 45 
6662.675....... I I I I | 2 2p.— 4d 
S000. 222. ...... I I I I 2 | 2p:— 4d 
SBO0.005.....+ 15 10 fe) 8 15 | Is —2py 
g60s.982....... 10 8 8 6 | 10 Is —2p, 
ee6.684....... I Roeataay | I 2p.— 38 
6360.725....... I Jost eeeees | | Pees I 2pr— 38 
DATA FOR TABLE I 
| 7; ; 
Exposure Number Volts Milliamperes | Menthol 
DR a eos sitnlecins 5000 eo 4 hours 
Se ee ; 100 150 | 30min. 
ee a 50 | 40-140 50 min. 
Deere eae icews 30 65 65 min. 
Wea ekereies aio Guat oniecatereud 8-10 2000 125 min. 


The wave-lengths in Table I are in international angstrom 
units. The spark lines given by Schillinger,* having been cor- 
rected from the Rowland to the international scale, are represented 
by six significant figures while the values for the arc lines from 
d\ 4390.14 to \6160.725 are those recently measured by Datta.? 
The remaining values are given to four or five figures and are the 
means of several determinations from the quartz prism spectro- 
grams; they are probably not in error by more than one unit in 
the last place. Approximate values for the wave-lengths of 
many subordinate series lines were given by Zickendraht* and 
our numbers are in fair agreement with his for the higher terms. 

It is seen from this table that the critical voltage for the excita- 
tion of the enhanced spectrum of sodium is about 30 volts, in the 
neighborhood of the values for La and La. 

t Schillinger, Sitzungsberichte der Kais. Akad. d. Wiss., 118, 605, 1900. 

2 Datta, Proc. Roy. Soc. (A), 99, 69, 1921. 

3 Annalen der Physik, 31, 233, 1910. 











Td 














WAISSVLOG dO Valo 
| 
eW O09 
SLIOA OV 
. ', Ef YS : 
3 3 2 & | om Q 
| j | | ; 
‘eUW OOE 
SLIOA G2 
” 
WOALIFIS = OFINVHNI 5s) 
“ 
eW ODN!) 
SLIOA L 
/\ 
o¢ WAALIIAS Dav 
$3 
ew OL 
SLIOA L 
38 ¥: a ce —— 
© VY © ¢¢ <¢ - 5 
re ne $F nm m * 


SLIOA 
se 





WOALITAS ANI? FTON/C 
Ill ALV Id 





Na AND K IN A LOW VOLTAGE ARC 157 


POTASSIUM 

Plate III shows the three stages in the development of the 
potassium spectrum. Table II contains the wave-lengths and 
estimated intensities of the lines at various exciting voltages from 
3.5 to 40 volts. From the resonance potential, 1.6 volts, to the 
ionization potential, 4.3 volts, the single-line spectrum is emitted, 
consisting of the pair \ 7699 and \ 7664. Above the ionization 
potential the complete arc spectrum is excited while in the neigh- 
borhood of 20 volts (viz., near Ma and Ma) the enhanced lines are 
prominent. <A wider slit was used for this series of exposures than 
for sodium and the results are therefore more favorable for repro- 
duction although less advantageous for precise measurement of 
wave-lengths. The values for the wave-lengths in Table II are 
compiled from the measurements of Shillinger and Datta. Only 
the stronger spark lines are listed, as these serve sufficiently to show 
how the intensity of the enhanced spectrum rapidly increases when 
the critical voltage is exceeded. A great many faint spark lines 
were observed in addition to those given by Shillinger, a fact also 
noted by McLennan,’ but their wave-lengths are not known with 
sufficient accuracy to justify publication. 

While the main issue of the present paper concerns the excita- 
tion of enhanced spectra, the second and third spectrograms 
of Plate III, showing the arc spectrum of potassium, are of con- 
siderable interest. Both of these exposures were made at 7 volts, 
the only difference being that an ionization current of 70 milli- 
amperes was employed for the upper spectrogram and 1000 milli- 
amperes for the lower. In the latter case the pair \ 4642 and 
d 4641 appears as one of the most prominent lines in the spectrum. 
This pair has the notation y=1s—3d, representing an interorbital 
transition where the change in azimuthal quantum number is 
according to Sommerfeld two units.? It accordingly contradicts 
the Bohr or Rubinowicz principles of selection as applied by 
Sommerfeld to spectra of the non-hydrogenic type While the line 


t Proc. Roy. Soc. (A), 100, 182, 1921. 
2In a paper by Foote, Mohler, and Meggers, to appear in the Philosophical 
Magazine, this and other exceptions to the principle of selection are considered in more 


detail. 
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TABLE II 








SPECTRA OF POTASSIUM | 
| 
s 4 : I ¢ ; Notation for 
Wave-Length | volts Volt Volt Volt Volts Volt Arc Lines 
9002.22. .... I [— I SPr,2 | 
3034.6. ; I I I I I 7 Pry2 
3002.44. 3 : 
3102.03 , i ; I Op: 
02.25/°°"*| . , : . I »p 
05.31 2 I I 
29.39. I I I 
69.59. I I I 
67.9! 2 I I 
3190.42 2 I I 
3201.96 | I I- 
17.01 : P _ I 5 
1 sO 4 5 5 5 5 : ct 
: j 
20.80 I I 
3290.93 2 2— I 
£329.03... 2 I [= 
45-7- 5 4 y 
55.33 I I 
60.16 I I 
63.25 
8) = 2 
04.50 . - 
73.84. 2 2 } 1 
50.97 2 | 2 I 
3385.2 | ¢ 2 I 
3404.57 5 2 
22.20. | I I I 
33.01. | I ja 
40.30. ..<% | 2 I i— 
"799 Iis— p 
40.722) 5 6 8 8 8 ore 
3447-701 IS— 4) | 
3530.83. | § | 5 3 I 
3609. 20. 2 I io 
18.43. 3 | 2 I 
68.91. I 
| 
76.01. 2 } I 
3681. 52 | 4 3 2 
3716.51 I I 
21.93. I I 
39.13. I [— 
44.62. Ref Soe ; 
-\- 7 
07.37. | 3 2 I 
3753.19. | 4 3 2 | [= 
3800. 00 3 2 I 
17.99. a 6 5 3 I } 
61.87. 4 2 I 
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TABLE II—Continued 


a 40 25 2c 16 6 3.5 Notation for 
, —— Volts Volts Volts Volts Volts | Volts Arc Lines 
73.78. I I 
3878.31 } 1 | 
3897.86. 7 5 3 2 
3923.65. 3 2 I- 
42.86 3 2 I 
| . | 7 
Se. 27. f Z 2 I 
60.75. 2 I I 
72.0%. 2 I 
3995.05 ‘ 2 I 
| 
4001. 20. } S 3 2 I } 
12.07. 
oe I 
44.140... 6 7 8 8 8 IS—3p: 
47.201. | 4 5 0 18) 0 IS—3p2 
57-95. 3 I es 
4114.96. a x 4 2 I— 
34-72. |} 5 4 } 2 I 
49. 24. 5 4 2 I 
4180. 23. 9 5 4 2 
4210. 21 2 I I 
. } - = . 
me ew - | | 
44.93 5) | 7 
25.60 5 | 4 | 2 I 
4263. 32.. o 5 2 t 
4305.01. 5 4 2 I 
o9.668.....]. 9 6 4 2— 
4388 .13.. 5 4 2 I 
4423.72.. I 
4400.60. | 2 | 
4505.39. 4 : I | 
4008. 31. 4 2 
1.585 | | f1s—3d’ 
4+ 5°35 2 4 : 4 8 3d 
42.172 1s—3d 
4659.58.. I 
4741.6... ie I— | 2pz—10s 
cf. ae a I I I ..| 2pa—11d 
.6..%. I— 2pi1—10s 
| | | | 
a I I I I 2p:—11d 
86.80.... I I | I I | 2p2— OS 
4791.08.. I I I 2 | 2p,—10d 
4800.16. I I I | 2p:— 9s 
05.19 | I I I 2 | 2p:—10d 
25.900 | 5 3 } 2 | ' 
49.88 | I— I I 2 2 | 2pa— 3s 
56.03 I 2 2 3 3 | 2pa— od 
63.01 [= | I I 2 2 2pi— 8s 
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TABLE II—Continued 
| ee a E aa : PY 
— - 40 25 20 16 | 6 | 3.5 Notation for 
Wave-Length | Volts | Volts | Volts | Volts | Volts | Volts | Arc Lines 
|) -_—— ee 
4869.70.... .| I 2 2 3 3 2p:— od 
4941.964....| I 2 2 2 |} 2 2p.a— 75S 
43.02.....| 1 eer me errr Migweledawee 
50.850... 2 3 , 3 |} 3 Jo seseee | 2pa— 8d 
56.043 I 2 i 3 2 2 | 2p:— 75 
4965 .038 | 2 | 3 3 3 3 | 2p:— 8d 
5005.34 | 2 & ie 
55.78 | I eee Geese Gennes Gentnnnn Gennnnnn Geet n nnn nnn 
84.212. 2— | 33> 3— 3 4—- | 2p.— 6s 
97-144 . | 3 3 3 4 2p.— 7d 
5099. 180.... = | ioe - 3- | «4 2p:— OS 
6It2. 204... 2 3 3 s 4 2p1— 7d 
5323.228.... , 3- 3— 4> 4> 2pa— 55 
39-670....) 2— : hae a 4- ‘~ 2pi— 55 
42.974....| 2 | 3 3 4 4 2p.— 6d 
5359-521. ..| 2 | 3 3 4 4 oes 118 2p, —6d 
7664.94..... cra amlene see A, ee ee Date eis 6 IS— 2p; 
7699.01.....| Pn dares ; ; : , 4 IS— 2p, 
DATA FOR TABLE II 
Exposure Number Volts Milliamperes P cade 
cecilia leat ea re aren Ge eee eee Suet tea ES 
2 40 500-600 Io min. r 
7 25 0-300 32 min. 
heared - ; 20 100-700 15 min. 
eiemees ita 16 o-500 20 min. 
ae en 6 1200 15 min. 
"NER eee 24.5 12 7 hours , 





is well known, its excitation has been thought to be occasioned by 
the presence of the electrostatic field—an incipient Stark effect. 
However, in the present case, with the special form of discharge 
tube employed, the radiation was observed only in the space shielded 
from the applied electrostatic field—itself only 7 volts. On 
referring to Table I and Plate II, it will be noted that the corre- 
sponding line \ 3427 appears in sodium. Its intensity increases 
very rapidly as the ionization current is increased, independently 
of the voltage, provided this exceeds the ionization potential. 


t Sommerfeld, Annalen der Physik, 63, 249, 1921. 
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The importance of this result has been pointed out in the paper 
referred to above, in which we suggest that the explanation may 
involve a reconsideration of the method whereby single azimuthal 
quantum numbers have been assigned to each of the s, p, b, and d 
terms. 

The following table summarizes the values of the critical 
voltages characterizing the three-stage development in the spectra 
of sodium and potassium. 





Element Volts Notation | Spectrum 
ee | 2.09 | ms—2p | Single pair 
5.11 Is | Arc lines 
| CE PP ae | Enhanced linesand L- 
radiation, La=35, 
| La=3o0 
Potassium......... | 1.60 1s—2p | Single pair 
4.32 Is Arc lines 
| ~e eCateiocte ubatea Enhanced lines and 


M-radiation 
Highest frequency in enhanced spectrum: sodium=14 volts; potassium= 
Ir volts. Total work of double ionization: sodium=49 volts; potassium = 34 volts. 


BUREAU OF STANDARDS 
December 1, 1921 








THE VACUUM-SPARK SPECTRA OF THE METALS! 
By EDNA CARTER 
ABSTRACT 

Vacuum-s park spectrograms of Ca, Mg, Cd, Ti and Fe were obtained with compar- 
ison arc and spark spectra, on portrait films, using for the most part a 20-cm concave 
grating spectrograph. The Fe spectrum was also photographed with a dyed film from 
\ 4000 to A6600 A. The new spectra show no striking new characteristics. In the 
case of Ca and also of Mg (see Plate [V), the vacuum-spark spectrum is practically 
identical with the spark spectrum in air; in the case of Cd (see Plate V), the arc lines 
are somewhat more intense; with Ti this tendency to strengthen the arc lines is even 
greater, and the vacuum-spark spectrum of Fe is more like the arc spectrum than the 
spark. In general these spectra resemble the luminescence spectra produced by 
cathode-ray bombardment in a high vacuum, and it is probable that the conditions of 
emission are very similar in the two sources. 


I. INTRODUCTION 

The production of a brilliant spark between electrodes very 
close together in a high vacuum was observed by Rowland? and 
also by Wood. The nature of the discharge was studied by 
Loving.* He concludes in regard to the spectrum produced that 
it is characteristic of the anode, the cathode having no effect, 
and that it is not analogous to either the spark or the arc. The 
spectra which he obtained were not sufficiently intense for an 
exact determination of their character. 

Millikan’ developed independently a method of producing this 
discharge; and, by using it as a source in a vacuum spectrograph, 
succeeded with his co-workers in extending the spectra of various 
elements into the extreme ultra-violet. 

It was of interest to see whether a careful study of this source 
would disclose any peculiarities of metallic spectra in the usual 
photographic region. A comparison of the vacuum-spark spectra 
with the arc and the spark spectra of typical metals was therefore 
undertaken. 

t Contributions from the Mount Wilson Observatory, No. 219. 

2 Rowland, Physical Papers, p. 574. 3 Physical Review, 5, 1, 1897. 

4 Astrophysical Journal, 22, 285, 1905. 

5 [bid., 52, 47, 286, 1920; 53, 150, 1921; Physical Review, 12, 168, 1918. Science, 
19, 138, 1919. 
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2. EXPERIMENTAL METHOD 


The vacuum sparks were produced in a glass bulb about 12 cm 
in diameter provided with a quartz window, which in the case of 
the more easily vaporized metals was removed to a distance of 
30 cm on account of the sputtering. The iron electrodes consisted 
of cylinders about one cm long and one cm in diameter, tapering 
to blunt points. These cylinders were screwed on the ends of 
aluminum rods, which projected to the outside through side tubes 
and were sealed in with sealing wax. The iron electrode was 
retained as cathode throughout the experiments, while anodes of 
the other metals were made by wedging a piece of the metal into 
a hole in the end of the aluminum rod. The distance between the 
electrodes was about 1.5 mm in the case of iron, but it amounted 
to 3 or 4mm with some of the other metals, which were largely 
vaporized during the exposure. 

Evacuation was effected with a Langmuir diffusion pump, the 
mercury vapor being kept away from the discharge by a trap 
immersed in solid CO, and ether. 

The discharge was produced with a large induction coil used 
with an electrolytic interrupter. A current of from 15 to 20 amperes 
was passed through the primary. To secure the necessary capacity 
across the secondary, four parallel plate condensers, consisting of 
copper plates about 56 by 61 cm separated by 5-mm glass plates, 
were connected in series. The whole was immersed in oil to 
prevent leakage. 

About eight hours were required to obtain a spectrogram of 
iron, because of the tendency of the discharge to pass over into an 
ordinary gas discharge, with the liberation of gas from the elec- 
trodes, and because it was necessary to keep the sealing wax joints 
cool. The total duration of the spark discharge was estimated 
at about thirty-five minutes for iron and not more than ten minutes 
in the case of cadmium. The iron required a much higher potential 
gradient between the electrodes for the production of the spark 
than the other metals. 

A spectrograph with concave grating of 20-cm radius and ordi- 
nary portrait films were used in obtaining most of the spectrograms. 
The iron spectrum was photographed in the region 4000-6600 A 
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by using an aesculin filter and films sensitized with a Wallace 
three-dye solution. Fairly rich spectra of iron and titanium 
were obtained also with a one-meter concave grating spectro- 
graph. 

The vacuum-spark spectra of magnesium and cadmium with 
comparison arc and spark spectra are reproduced in Plates IV and V. 
3. DISCUSSION OF RESULTS 

The vacuum-spark spectra of calcium and magnesium are prac- 
tically identical with the spark spectra in air. 

The cadmium vacuum-spark spectrum retains arc lines in 
somewhat greater relative intensity than they possess in the 
ordinary spark. The line \ 3261 is an example. This line came 
out also with marked intensity in the luminescence spectrum pro- 
duced by cathode-ray bombardment." 

Titanium shows an even greater tendency to retain the arc 
lines. The spectrum is plainly an enhanced one, however, as is 
evidenced by the brightness of the lines \\ 4164 and 4172, for 
example, the ratios of whose intensities in the spark and the arc, 
as given by Exner and Haschek, are 20:2 and 15:1, respectively. 

The vacuum-spark spectrum of iron, on the other hand, seems 
to resemble more closely the arc spectrum. The low-temperature 
lines found in the arc are missing, and the lines which are stronger 
at the poles and in the core of the arc are relatively more intense 
in this source. Even in this spectrum a close examination reveals 
the presence of the enhanced lines, such as AA 4508, 4515, 5018, and 
5316. Ha comes out very strongly and sharply in the photographs 
of the iron spectrum made in this region. 

The varying tendency exhibited by these different metals 
toward the production of enhanced lines suggests a resemblance 
of the vacuum-spark spectra to the luminescence spectra produced 
by cathode-ray bombardment in a high vacuum. It is probable 
that the conditions of emission are very similar in the two sources. 
The difficulties in the way of obtaining rich spectra are con- 
siderably less in the case of the vacuum-spark. 


Mount WILSON OBSERVATORY 
September 1921 


t Astrophysical Journal, 44, 303, 1916; 49, 224, 1919. 
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THE MASSES AND DENSITIES OF THE STARS! 
By FREDERICK H. SEARES 


ABSTRACT 


Mean magnitudes and masses of visual binaries and dwarf stars of the different 
spectral types—A statistical comparison of the absolute magnitudes (M,) of 105 


visual binaries, the hypothetical magnitudes (M_,) of these and 400 other visual 
binaries (computed by Jackson and Furner, assuming the combined mass, u=y:+ 4a, 
of each binary to be 2), and the absolute magnitudes (M,) of 1152 single stars (includ- 
ing some binaries) by Kapteyn and Adams, gives the mean difference (M,—M,) 
= (M,—M,)—0.3*0.1 for each type; that is, M,=M,—0.3+0.1. This indicates 
that the selection used is sufficiently homogeneous, also that the variation of the 
mean magnitude of binary star components with spectral type is the same as that 
of single stars. The values of M, (Table IV and Fig. 1) for dwarf stars of types 
Bo, Ao, Fo, Go, Ko, and Ma are —1.60, +0.70, 2.40, 4.35, 5.90, and 9.80 respectively. 
The geometric mean mass for each spectral type of visual binary was next computed 
from the formula log u=log 2—0.6 (M,—0.3—M,). The values of u, slightly cor- 
rected for the effect of selection as determined by a study of the binaries of known 
parallax, give, when multiplied by 4 to reduce to single stars, a first approximation to 
the geometric mean masses of single stars along the dwarf branch of Russell’s diagram, 
the values for types Bo, Ao, Fo, Go, Ko, and Ma, coming out 10, 6.0, 2.5, 1.0, 0.68, 
and 0.59, respectively, in terms of the sun’s mass (Table IV). The data for visual 
binaries of known parallax indicate that the distribution of log u is approximately 
Gaussian and that the dispersion in mass is small, probably half the stars of each 
type having masses between 0.7 and 1.5 times the mean (Table VIII). 

Mean velocity and energy of dwarf stars of different spectral types—The mean 
square space-velocity corresponding to the mean magnitude was determined from the 
available data for each type, and these values were multiplied by the corresponding 
mean masses. The mean energy for each type so obtained is practically the same for 
all types from Ao to Ma, in spite of a 10 to 1 variation in mass (Table IX). The 
logarithm of the mean kinetic energy constant is 3.270.036 (average deviation per 
type) in terms of sun’s mass and km/sec., corresponding to 3.5 X10 ergs. The mean 
radial energy is also approximately constant for all except the B-type stars, which 
seem to be abnormal in this as in other respects. Whatever the cause may be, the 
distribution of energy among these stars seems to approximate equipartition. 

Masses of stars as functions of magnitude and spectral type were then computed 
from the values of the mean velocity for each magnitude and type, assuming the mean 
energy independent of magnitude as well as of type (Table X). In Figure 2, curves 
are given showing M as a function of type for various constant values of the mass. 
To get more reliable data for the giants, the masses of 28 Cepheids were computed from 
the period, magnitude, and surface brightness of each, using a formula based on the 
pulsation theory (Table XVI), and Figure 2 was revised to agree with these values, 
giving Figure 3. If the masses computed in accordance with Eddington’s theory of 
Cepheids (Table XIX) are adopted, a further slight revision will be necessary. The 
constant mass lines would be less irregular if bolometric magnitudes were used instead 
of visual. For the Cepheids, the probable dispersion in mass for a given M and type 
is only about 20 percent. This result taken with the data for visual binaries indicates 
that in general the luminosity of a star is very closely determined by its spectral 
type (temperature) and its mass. 

Effective temperature, color index, and surface brightness of stars of different spectral 
types.—The effective or black-body temperatures of the giants as determined from 





* Contributions from the Mount Wilson Observatory, No. 226. 
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Wilsing’s spectrophotometric measures of c/T for 199 stars, vary from 2890° for Mc 
to 10,500° for Bo, while those of the dwarfs as determined from Mount Wilson color 
indices (provisional values, hitherto unpublished, also given for the giants) reduced to 
the same system, vary from 3330° for Ma to 6080° for F5 (Table XII). The 
difference between dwarfs and giants of the same type is too great to be neglected, 
the later-type dwarfs being considerably hotter. The surface brightness 7 as a function 
of temperature may be computed (1) from the difference between visual and bolometric 
magnitudes as determined by Eddington (j=M—M,—10 log T+ const.) and (2) 
from Hertzsprung’s formula derived from Planck’s radiation law (j=2.3 (¢/T)°% 
+const.), both formulae giving practically the same results (Table XI). Combining 
these values with the temperatures of the stars of various types, we get the surface 
brightness as a function of spectral type, the results varying from 4.81 for T= 2890° 
to —2.28 for T=10,500°, expressed in magnitudes, sun=o (Table XII). As a check 
on these values, the surface brightness of three stars, a Orionis, a Boétis, and a Scorpii, 
whose diameters have been determined directly by Pease, were computed and found 
to agree with the interpolated values within +o™2 in each case (Table XIII). 

Diameter and density of stars as a function of magnitude and spectral type were 
readily computed from the surface brightness, absolute magnitude, and mass (Table 
XIV). For dwarfs the diameter decreases from 6.8 to 0.54 (sun=1) and the density 
increases from 0.045 to 5.4 (water=1) as the type changes from Bo to Ma, while for 
giants of zero magnitude the corresponding changes of diameter and density are from 
3.2 to 66 and from about 0.3 to 10-5. When M is plotted as a function of type for 
various constant values of the density (Fig. 2), the equidensity lines are straight 
and nearly parallel and quite regularly spaced, so that the equation: log p=—1.22S 
+o.57M+1.11, where S is the type number beginning with o for Bo, represents the 
results quite closely for a wide range of density. For dwarfs of types B to F, the 
values agree satisfactorily with the values found by Shapley for the densities of 
eclipsing variables of these types. The giant part of the diagram, however, is less 
reliable, and is later modified in Figure 3 to agree with the values for 28 Cepheids 
(Table XVI) computed from their periods by the use of the theoretical formula: 
log p=—2 log P+const. Adoption of Eddington’s formula for Cepheids would lead 
to only a slight further modification (see Table XVII). * 

Relation of the observed variation of mean mass and of M with spectral type to the 
problem of stellar evolution.—For the dwarf branch, the decrease in mass from 10 for Bo 
to 0.6 for Ma may be partly, and perhaps wholly, accounted for by selection, the 
earlier and brighter stars being observed throughout a larger volume; further, only 
massive stars can attain high temperatures. The possibility that mass may decrease 
with loss of energy by radiation is also suggested, a possibility which agrees with the 
theory of relativity and does not necessarily conflict with Newtonian mechanics. 
For the giants, the observed change in M with spectral type (mass constant) is diffi- 
cult to account for. Probable changes in mean atomic weight are insufficient. 
Changes in the opacity are doubtless involved, and likely also the source of energy, 
which almost certainly is not wholly gravitational and probably not independent of 
the time; for instance, the probable degree of ionization involves an amount of energy 
greater than seems to be available from gravitational forces. The intersection of 
equal-mass lines with the dwarf branch at large angles (Figs. 2, 3), however, does 
not conflict with the gravitational contraction hypothesis. 

Eddington’s radiative equilibrium theory of giants.—Although this involves the 
assumption that the ratio of the energy flow to the gravitational acceleration is constant 
throughout a star, whereas the ratio probably decreases considerably from center 
to outside surface, the calculated variation of magnitude with mass for a given 
spectral type is in fair agreement with observed values (Table XXII), but for a 
given mass the absolute magnitude is not constant as required by theory. For 
Cepheids, the theory gives the formula: log p=—2 log P—2 log (ya)i+const., and 
this combined with the empirical relation, mass is proportional to radius, which 
Eddington’s results suggest as a condition for Cepheid variation, gives: M= 
—s5 log P+j—5 log (ya)+const., which agrees with the observed variation of M 
with P for P> 2 days, quite closely (Table XVIII and Fig. 5). This shows the 
consistency of Eddington’s theory of Cepheids with his theory of radiative equilibrium 
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and tends to confirm the relation of mass to radius, which is equivalent to the condi- 
tion that the average heat content per unit mass is constant. Various consequences 
of this condition are given. The theoretical values of mass and density given by 
the pulsation theory are close to those indicated by Figure 3 (Table XVII). 

Ionization and mean atomic mass inside stars.—Theoretical computations indicate 
that the temperature and pressure at various points inside a star change in such a way 
as to keep the degree of ionization of the same order throughout (Table XX); hence 
the mean atomic mass is also about constant. When we consider giant stars of 
constant mass but of increasing temperature, we find that while the densities change 
very greatly, the high degree of ionization in the earliest stages is maintained with 
but little change (Table XXI), so that the mean atomic mass should decrease but 
slowly as development proceeds. 

Comparison of spectroscopic parallaxes of visual binaries with the hypothetical 
parallaxes of Jackson and Turner shows systematic differences for their Table II 
(Table ITI). 

Information on the masses of stars has its origin in what can be 
learned of the masses of binary systems, spectroscopic and visual. 
Any discussion of spectroscopic binaries from the standpoint of 
mass is complicated by lack of definite information on the inclina- 
tion of orbital planes, while in the case of visual binaries the lack 
is that of accurate knowledge of stellar distances. Approached 
from either standpoint, the problem, with rare exceptions, must 
be treated by statistical methods whose application presupposes 
an abundance of data. Any solution is therefore involved in 
some uncertainty, which is not lessened by the fact that a transfer 
of results thus obtained to the stars at large involves assumptions 
of comparability: we do not know, for example, how the mass of 
an average binary is related to that of an average single star of 
the same type and luminosity, or even that the relation is the same 
for all classes of stars. The uncertainty, however, is probably not 
great; at any rate, we can proceed only by making such assumptions. 

The following discussion is based on the visual binaries. These 
lead to values of the masses of stars along the dwarf branch of the 
Russell diagram, including those of B-type spectra. The average 
mass decreases continuously with advancing spectral type. The 
combination of these results with the space velocities of large num- 
bers of stars recently published by Adams, Strémberg, and Joy 
shows that the mean kinetic energy of translation of the A’s and the 
dwarfs is approximately constant. The principle of equipartition 
is then used to derive values of the masses of stars of other luminosi- 
ties. The distribution curves thus obtained for mass as a function 
of spectral type and absolute magnitude are a first approximation 
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in so far as they concern stars of highest luminosity, but in the 
vicinity of the dwarf branch they should be reliable within a 
moderate percentage. 

Valuable information on density has at various times been 
obtained from the eclipsing variables, but it has been difficult to 
correlate these results with other stellar characteristics, except 
loosely, because here again the distances, and hence the luminosities, 
of individual stars are inadequately known. Density, however, 
depends in a simple way on mass, total luminosity, and surface 
brightness. The surface brightness can be computed from radia- 
tion formulae, the results being checked in part by their agreement 
with Pease’s measures of angular diameters. Hence with mass 
once known as a function of spectral type and absolute magnitude, 
the values of the density follow as a matter of course. The results 
thus found for mass and density are subject to control and readjust- 
ment with the aid of the Cepheid variables, whose relative densities 
and masses can be calculated directly from their periods, spectra, 
and surface brightness. 

The lines of equal mass for the giant stars in Russell’s diagram 
are irregular curves, considerably inclined to the axis of zero 
absolute magnitude. If therefore a giant star in its evolutionary 
development follows an equal-mass line, its luminosity does not 
remain constant during the giant stage as has often been supposed. 
The curves of equal density are approximately linear and parallel 
to the line of maximum frequency for the dwarfs. The correlation 
of mass with spectral type along the dwarf branch can be explained, 
in part at least, as a consequence of selection, although other 
factors may also enter. The concluding sections of the paper dis- 
cuss the relation of these results to certain phases of Eddington’s 
theory. 

I. THE METHOD 

The mean masses of the visual binaries of different spectral 
types may be found by a statistical comparison of the absolute 
magnitudes corresponding to their hypothetical parallaxes with the 
absolute magnitudes of stars of known parallax. The relation 


a3 
a ey (1) 
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leads at once to 
log u=k—0.6M (2) 
where 
k=3 log a—2 log P+0.6m+3 


and where the absolute magnitude is defined by 
M=m-+5+5 log z. (3) 


For a hundred systems of known orbital elements Jackson and 
Furner’ have used equation (1) to calculate the hypothetical 
parallax m, corresponding to an assumed mass y=2. For about 
450 systems of unknown elements they have calculated hypothetical 
parallaxes from a formula depending on changes in position angle 
and distance, the assumption for the mass again being y=2. 
The absolute magnitudes M, corresponding to the hypothetical 
parallaxes for the two groups of stars are given in Tables I and II 
of their paper. 

For any system of known elements, by (2), 


0.30=k—0.6 M, 
and by combination with (2) itself, 
log u=0.30—0.6 AM (4) 


where AM=M-—M,. Hence for a group of stars? 
log u=log p=0.30—0.6 " 


AM=M_M, (s) 


in which p represents the geometrical mean mass. 

Equations (4) and (5) are the basis of the present discussion. 
Since the hypothetical parallaxes derived from changes in position 
angle and distance involve the mean inclination of the orbit-planes, 
the individual values of 7, thus found will differ from those obtained 
with the rigorous relation (1). Equation (4) therefore holds for 
individual stars only in the case of systems of known elements. 
Formulae (5), however, are valid for any representative group of 

t Monthly Notices, 81, 2, 1920. 

2 Upon re-reading Russell’s address, ‘Relations between the Spectra and Other 
Characteristics of the Stars,”’ Publications of the American Astronomical Society, 3, 22; 


Popular Astronomy, 22, 275, 331, 1914, I find that'in principle the method is essentially 
one used by him. 
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stars, whether their elements are known or not, and can therefore 
be applied to objects in both the lists of Jackson and Furner. 

To determine the mass from equations (4) and (5), the values of 
M, or at least of 7, must be known for the various spectral types. 
For a certain number of systems direct determinations of distance 
are available. Hence, M is known, and the use of (4) gives at 
once values of uw for the individual stars. These results will be 
affected by large uncertainties, but the means for groups of stars 
should be accurate within a small percentage. For most of the 
binaries, however, only the hypothetical absolute magnitude M, is 
known; no other approximation for M is available, and we must 
proceed indirectly. This we might do by using values of M 
derived from single stars of known parallax whose selection with 
respect to the stars as a whole presents the same characteristics as 
those of the binaries themselves. Similarity of selection in the 
two lists of stars—binaries and single stars—is, however, difficult 
to attain in practice. It will be sufficient if we can find a list of 
single stars whose differences in selection as compared with the 
binaries are the same for all spectral types. 

Let M, represent the mean absolute magnitude for a given 
spectral type of the single stars of known parallax, and let M be the 
corresponding magnitude of the binaries of the same type. We 
may then write 

M=M,+6M (6) 
where 6M indicates the influence of difference in selection. If this 
difference is the same for all types, 6M will be constant, and we shall 
then speak of the selection as homogeneous. 

The second of (5) then becomes 


AM=M,—M,+6M =AM,+6M (7) 


where M, is the mean hypothetical absolute magnitude of the 
binaries of the spectrum under consideration. The value of 
AM,=M,—M, is found by comparing the two lists of stars for 
which it is presupposed that the selection for different spectral 
types is homogeneous, while the constant 6M is given by 








5M =AM—AM,=M—M,—AM. (8) 
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applied to the systems of independently determined parallax 
referred to above. M represents here the magnitude corresponding 
to the measured parallax of a binary whose hypothetical absolute 
magnitude is M,. In other words, we derive 6M by (8) from the 
hundred or more binaries of measured parallax, calculate AM, for 
each type by comparing the five hundred and fifty-odd stars in 
the lists of Jackson and Furner with the homogeneously selected 
single stars of known parallax, form AM by (7), and, finally, solve 
for » by means of the first of equations (5). Since the question of 
homogeneity of selection can be put to a direct test, the only 
assumption underlying this procedure is that the relation of mass 
to luminosity is the same for both binaries and single stars, which 
may be accepted as plausible at least. 

The values of M, and M,, whose difference is AM,, are readily 
found by a graphical process. Plotting the values of M, against 
spectral type, we find for the binaries a frequency diagram similar 
to that first given by Russell. The late-type giants are not numer- 
ous, but the dwarf branch is well defined. The points of maximum 
frequency for the dwarfs define a curve which joins smoothly with 
the line of modal values of the B- and A-type stars. This curve 
expresses the variation of M, with spectral type, the giants of 
late type being disregarded altogether for the present. For M, a 
similar curve was found from the data in the list of spectroscopic 
parallaxes,’ supplemented by those for about 430 helium stars whose 
parallaxes have been determined by Kapteyn.? These two curves 
are shown in the central part of Figure 1. Their ordinates and 
differences in ordinates are given in the second, third, and fourth 
columns of Table IV. The data for the M, curve could have been 
greatly increased by including results from trigonometric paral- 
laxes. Little would have been gained in precision, however, for, 
as it is, the number of stars entering into the M, curve is more 
than twice that for the M, curve. 


t Adams, Joy, Strémberg, and Burwell, Mt. Wilson Contr., No. 199; Astro- 
physical Journal, 53, 13, 1921. 


2 Mt. Wilson Contr., Nos. 82, 147; Astrophysical Journal, 40, 43, 1914; 47, 146, 
255, 1918. : 
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Fic. 1.—The curve M, is the line of maximum frequency of M for the B stars 
of Kapteyn and the dwarfs from the list of spectroscopic parallaxes. M_, is the similar 
curve for the hypothetical absolute magnitudes of Jackson and Furner. The differ- 
ences in the ordinates, corrected for zero point, are shown by the curve below (AM). 
The remaining curves refer to the discussion in section 16. 
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If we may assume that the selection for different spectral types 
is homogeneous, the differences in the ordinates of the two curves 
represent the values AM, as a function of spectrum. Irregularities 
of selection, or a progressive change in selection with type, will 
produce deformations or a tilt of the curve for AM,. The curve 
for AM, formed directly from the curves for M, and M, and 
corrected for 6M, is shown in the lower part of Figure 1 (the points 
are explained later). Since the curve is smooth, we conclude that 
there are no marked irregularities of selection, or at least that 
their influence has been minimized in drawing the curves for 
M, and M.. 


2. INFLUENCE OF SELECTION ON THE STATISTICAL COMPARISON 


To examine the question of selection more closely, consider 
the numbers. of binaries and single stars of known parallax in 
each interval of apparent magnitude. Excluding the late-type 
giants and a few scattering stars at the ends of the sequence of 
apparent magnitudes, we find the data in the first two divisions 
of Table I. If the selection for the two lists were the same (not 
merely homogeneous), the ratios of the corresponding numbers 
for each type would not vary with the magnitude. Difference in 
selection for the binaries and single stars of any type implies a 
change in ratio with increasing magnitude; and, if the selection is 
homogeneous, the ratios for different types will show the same 
change. Absolute values of the ratios are of no significance, 
and for ease of comparison the ratios for each type have been 
multiplied by a constant factor which reduces that for median 
magnitude 6.5 to unity. Further, the numbers for the Gs5—K5 
binaries have been combined with the mean numbers for single 
stars of Fo-Fg and Ko~Kog types. 

The results in the third division of Table I indicate a satis- 
factory degree of homogeneity. The irregularities affecting the 
bright Gs5s—Ks5 stars are unimportant, since the individual groups 
include only two or three stars each. The only deviations requiring 
comment are those shown by the faint Gs—Ks stars and the faint 
B’s. ‘The latter arise from the fact that Kapteyn’s lists of B stars 
include none fainter than the 6th magnitude. The binaries 
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involved are ten in number, all of the 7th magnitude or fainter, 
and all but one or two of types B8 or Bg. Any error in the corre- 


TABLE I 


NUMBERS OF STARS—INFLUENCE OF SELECTION 


| Median Apparent Magnitude 
| - Total 
Spectrum | 2.5 3.5 4.5 ss 6.5 7.5 ° >.5 we 

} Binaries—Jackson and Furner 
dinininl | | ; — 
Oes—Bo.........- | I I 7 7 7 5 2 33 

| 
Ao “ee | 4 3 12 22 42 | 25 15 I 134 
Fo -Go(M>-+1) Oo 6 Is 20 66 78 4 2 257 
Gs -K5(M>-+3) fe) I 2 | 4 ra) 24 39 I Str 

Total Tere Peers : oe: © ere eer bi. . 505 
Single Stars*—Kapteyn and Adams 
l l l 
Cer-Bo...025.52.) 23 | 42 137 | 104 30 | ° fe) Oo} 426 
ee | oO 7 | 11 11 | 8 | 3 2 I 43 
Fo -Fo(M>-+1) 2 14 | 48 | 105 63 50 29 3 314 
Go -—Go(M>-+3) ° 5 10} 35 57 54 40 8 | 209 
Ko -Ko(.U>-+4) ° | I 2 14 26 | 43 62 12 160 
——— 
: ee | 1.263 


Ratio—Binaries/Single Stars 


| | 
} | 





eee 0.1 0.2 0.2 1.0] © x 
| A eee! See O.I | 0.2 0.61 1.0 | 1.6 | 1-5 0.2 
eee seer 0.4 6.21 6.3 L.o1s.5 * 2.8 0.7 
Gs-Ks5 } 1 | 3 | 08 1.0 | 2.2 3.5 0.5 

Ratio—Jackson and Furner/ Pickering 

| | | 
— Es ere: are 7 | 2 | 1.0] 0.52 0.08 | 0.00 
EEE ET CEP 3 2 1.0 | 0.29 | 0.0 0.00 
; : ‘ 
— rere ere 3 | 2 1.0 | 0.38 | 0.10 | 0.01 
Se 3 | 2 1.0 | 0.80 | 0.54 | 0.02 
| 

Ratio—Kapteyn and Adams/ Pickering 

| 
SEE rec mer earores aeene | 34 «| 13 1.0 | 0.00 | 0.00 | 0.00 
eee 7 117 | 4 1.0 | 0.18 | 0.03 | 0.02 
ae ore errs Tae ies 6 | 1.0] 0.25 | 0.05 | 0.01 
tees tran bea soaaca chen wt as | Be wes 3 4 | 1.0 | 0.33 | 0.08 | 0.02 
RRR ere ae Seen oo 0.4 Es I © | 0.46 | 0.34 0.12 

| | 





| 
| 





* The “single stars” include the brighter members of a few binaries. 


sponding AM, must be small, for there is no pronounced irregu- 
larity in the curve. The value of M, for Kapteyn’s earlier B’s 
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can scarcely have been affected by the limitation in apparent 
magnitude, for these stars are all so luminous that the limitation is 
of little consequence. 

Among the stars fainter than m=7.0 there is an excess of 
Gs5-Ks binaries which represents a considerable fraction of the 
total number of Gs5—Ks5 stars. The corresponding M, may be 
slightly too large. But, if appreciable, this effect must also be 
small, for as the stars actually observed in the intervals of fainter 
apparent magnitude are successively added, the value of M, 
increases very slowly. A direct test is afforded by the binaries 
of measured parallax considered in a later paragraph. 

The number of A stars in the spectroscopic list is small, but 
this does not seriously interfere with the determination of the 
curve of modal values because there are well-determined points 
for the early B’s and for the F’s. It does suggest, however, the 
advisability of an independent test of the selection of the much 
larger group of A’s among the binaries, whose numbers can be 
considerably changed without greatly affecting the ratios in the 
middle section of Table I. To cover this point and as a further 
general test, both groups of stars have been compared with the 
counts of W. H. Pickering.t Using means of his values for the 
galactic poles and the galaxy itself (except for the B’s, for which the 
galactic counts alone have been employed), we find the ratios in the 
last two sections of Table I. There is here some progression with 
type due to the fact that Pickering’s counts include the late-type 
giants; but the ratios for the A-type binaries are in general 
agreement with those of the adjacent classes, which answers the 
question as to the selection of the A’s. The large ratios for the 
bright B’s in the last section of the table show merely what appears 
directly from the counts, namely, that Kapteyn’s lists are not 
complete beyond the 6th magnitude. There seems also to be a 
defect of fainter A’s in the spectroscopic list, a conclusion appar- 
ently confirmed by the fact that the points for these stars lie above 
the smooth M, curve shown in Figure 1. It seems clear that none 
of the questions raised is serious and that we may accept the 


* Publications of the Astronomical Society of the Pacific, 33, 140, 1921. 
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values of AM, found in the manner described and proceed to the 
determination of 6M. 


3. ZERO POINT OF THE CURVE FOR MASS AND SPECTRUM 


It will be noted that the values of AM, determine the rate of 
change of mass with spectral type, while the constant, 6M, deter- 
mines the zero point. For the calculation of 6M only the spectro- 
scopic parallaxes have been used, since these are based on a homo- 
geneous system and are sufficiently numerous to afford an excellent 
determination. In forming values of M-—M, for individual stars 
to be substituted into (8), two points must be borne in mind: 
First, the values of M given in the spectroscopic list correspond 
to the most probable values of the parallaxes, and, because of the 
logarithmic relation connecting absolute magnitude and parallax, 
are not the most probable magnitudes themselves. To obtain 
the latter, a small constant correction of +0.08 mag. must be 
applied to the tabular values. Second, the values of M, given 
by Jackson and Furner usually represent the total luminosity of 
the system. To make them comparable with values of M in the 
spectroscopic list, they must be reduced to the brighter component. 

A summary of the results for 6M is given in Table II, the two 
lists of Jackson and Furner being treated separately. The values 
of AM, required for equation (8) were interpolated from the fourth 
column of Table IV. The adopted result is 6M =-—o.3, whence 
equation (7) becomes 


AM =AM,—0.3. (9) 





It will be noted that there is no marked progression in the values 
of 6M with type, which confirms, in part at least, the foregoing 
conclusion as to homogeneity of selection. This is also shown by 
the lower part of Figure 1, in which the curve is that representing 
equation (9), while the points correspond to the values of M—M, 
in Table II. With the exception of the constant term, the curve 
depends wholly upon the statistical comparison, while the points 
are derived from values of M and M, referring to the same star 


t Mt. Wilson Contr., No. 199, p. 15; Astrophysical Journal, 53, 27, 1921. 














THE MASSES AND DENSITIES OF THE STARS 177 


and are uninfluenced by any possible dissimilarity in selection. 
The agreement of the points with the curve is, however, within the 
uncertainty affecting their positions. 


TABLE II 


ZERO POINT OF CURVE FOR AM; 
































J anv F, Taste I } J anv F, Taste II 
| -_ : 
Sp. |M-M,| AM, 5M No. Stars}| Sp. | M—M, AM, 5M No. Stars 
| | | 
Az.....| —0.7 | —o.5 | —0.2 3 || Ag....| —o.8 | —0.4 | —0.4 II 
F2. —0.6 | —0o.2 0.4 14 | Fs....1 ©.0 0.0 0.0 rf 
Fs |} —o.r | 0.0 o.1 13 || Go....] +0.3 | +0.4 ] —o.1 
. | 3 ; 9 
Go -| 0.0 | +0.4 0.4 18 ! K3....| +o.3 | +0.8 | —o.5 9 
G4. | +o0.4 | +0.6 0.2 10 | 
K2. +0.4 | +0.7 0.3 7 || 
K8....| +0.5 | +0.8 | —o0.3 4 || 
| | 
’ a 
Mean and total..... —o. 29) 69 Mean and total...... —0.27] 36 











A further detail bearing on the question of selection is the 
circumstance that the constant 6M arises almost wholly from the 
use of the tabular values of M, for the derivation of the curve for 
M,.. Since these generally represent total luminosity, the reduction 
to the brighter component appears in 6M, and in the mean is nearly 
equal to the value found for this constant. The effect upon AM, 
of the obvious differences in selection shown in Table I is therefore 
very small. Consequently, the influence of small deviations 
from homogeneity in selection in passing from type to type may 
be regarded as altogether negligible. 


4. SYSTEMATIC DIFFERENCES IN THE HYPOTHETICAL PARALLAXES 


Although the two mean values of 6M given in Table II are 
practically identical, this was not the case when the calculation 
was originally made, the difference then amounting to several 
tenths of a magnitude. This discordant result has been traced to 
a systematic difference between the parallaxes of the two lists of 
Jackson and Furner, the character of which may be seen from the 
comparisons shown in Table III. For systems of known elements 
(J and F, Table I) the agreement of the mean hypothetical paral- 
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laxes with the spectroscopic results is excellent; the one large 
difference +0’014 is only 7 per cent of the corresponding 7, and 
would be reduced to +0”004 by the rejection of a single star, 
Bu 2100. 

TABLE III 


COMPARISON OF SPECTROSCOPIC AND HyPOTHETICAL PARALLAXES 





} : 
ian ae i s | J ann F, 
| J anv F, Taste I | J anv F, Taste II | Taste I—Taste Il 
| | 
. lc, | neeen | | | 
Mean x=, | No Mean | 7, No Mean , en ee 
al | ™ | | af | s 
| | 
” ” . a a . | 5.8 . | 
<0.020 eer O-.O16/+0.00I] 7 ©O.O15}/—0.002} 5 | O.O15}/—0.001] 6 
© .020-0.049 . fe) 030|—o 001} 28 | 0.036]/—0.009] 16 | 0.024/—0.007| 12 
©.050-0.099...... ©.0605|/—0.003} 22 | 0.064/—0.012] 13 | 0.043/—0.012] 8 
>o.100 , ..| O.198|/-+0.014| IO | 0.209)/—0.071| I a ee oe 
—— _ —$—$$$_—___—__—_ —_—__—— — — San, es, —— _ a 
* The intervals for x, are =0%020, 0%021-0%029, and S0%030. 


: 


For the stars in the second list, however, the agreement for the 
first group alone is within the uncertainty; for the others the 
differences are systematic and an important percentage of the 
parallaxes themselves. Further, the same difference is shown by 

a direct comparison of the hypothetical parallaxes of twenty-six 
stars which occur in both the lists of Jackson and Furner. The 


means, arranged in three groups, are in the last division of Table ITI. 


nn aan 


For half of these stars no spectroscopic parallaxes are available, 
and to that extent the second comparison is independent of the first. 
The difference, which is about 20 per cent of the parallax, corre- 
sponds to a difference of 0.4 mag. in the two values of 6M and ' 
fully accounts for the discordance. For the final determination 
of 6M, the hypothetical parallaxes from the second list were reduced 
to the spectroscopic system in accordance with the indications 
of Table III, in order that all the results may refer consistently 
to this system. The resulting value of 6M is that given in the 
Hl right-hand division of Table II.' 





tT am greatly indebted to Mr. Adams for the spectroscopic parallaxes of nineteen 
binaries in the first list of Jackson and Furner (six stars also occur in the second list) 
which have been obtained since the reductions described in the text were finished. 
These reductions have not been revised, but I have satisfied myself that the inclusion 
of the new data would not change the results in section 4 or the following sections 
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5- PROVISIONAL MASSES—-EFFECT OF SELECTION 


The derivation of values of the masses now requires only the 
substitution of AM from (g) into the first of (5). The results for 
log w are in the fifth column of Table IV. The first and last 


TABLE IV 


ADOPTED GEOMETRICAL MEAN Mass (©=1), SPECTRAL TYPE AND ABSOLUTE 
MAGNITUDE 
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K5 7.10 | 6.30 0.80} 0.00 |+ 3 0.03 £8 9.79 | 0.62 
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values depend upon extrapolations, and it is possible that those 
for Bs and Ao are also affected by considerable uncertainty. So 
far as the statistical comparison is concerned, the remaining values 
should be dependable within 5 or ro per cent. 

These results can be accepted, however, only with some reserva- 
tion, for it is easily seen that the visual binaries are a selected 
class of stars in the characteristic of mass itself, and especially 
is this true of the more distant objects. Since motion is detected 
only with difficulty, or not at all, in case the period is very long, 
equation (1) shows that the selection of large masses is favored to 
an increasing degree as more and more distant stars are considered. 
For a given value of a the factor a*/m increases rapidly with 
increasing distance; the smaller masses are therefore excluded in 
succession because they correspond to periods so long that motion 
has not been observed. For large values of 7 a wide range in mass 


in any essential particular. Thus the first two differences in the first division of 
Table III would be +0%002 and —o’o0o02 for seventeen and thirty-seven stars, respec- 
tively. The differences in the third division would not be changed by the inclusion 


of the six stars also occurring in the second list. 
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will be included, but since there is a practicable lower limit to a, 
only the most massive among the more distant stars will appear 
in our lists. 

Another and perhaps even more important form of selection 
arises from the fact that only the most luminous of the distant 
stars get into our lists; the others are too faint to be seen. This, 
too, excludes in succession the smaller masses with increasing 
distance, and is related to the correlation of mass with spectral 
type discussed in section 16. 

The influence of distance on the selection of mass seems to be 
shown by the data now available. To test the matter, consider 
the sixty-nine binaries of known parallax’ in the first list of Jackson 
and Furner. For these, the individual masses, yp, can be calculated 
by (4). Comparing these values with the mean mass for stars 
of the same spectrum, as interpolated from Table IV, we obtain a 
series of residuals A log u=log w—log uw. When grouped according 
to parallax, these lead to the mean systematic deviations shown in 
Table V. A progression with increasing parallax, in the direction 


TABLE V 


DEPENDENCE OF MAss ON PARALLAX 


®. Mean Sp. | No. Stars Alogu | Curve 
“errr F2 10 —0o.08 | —o.II 
.026.. aie F6 10 —0.12 —0.08 
co ee G1 10 —o0.08 —0.07 
Ms ak sbineatere ofl F7 10 +0.03 | —0.03 
ce ee ie. 10 | 70-04 | 0.00 
Ses cia ais miei Gi 10 | +0.04 +0.04 
CBIR Civics G8 re) | +o.12 | +0.12 








anticipated, is evident, the negative sign of the differences associated 
with the more distant stars indicating that the corresponding 


* The thirty-six stars of known parallax in the second list are not included because 
of the systematic difference in the parallaxes of the two lists. In the preparation of 
Table II, where only a mean result for all the stars was finally used, the systematic 
difference can be corrected with sufficient precision. Here, however, we are con- 
cerned with individual stars, or at most with small groups, and it is not certain that 
the mean correction is applicable. Then, too, the fact that the second list of Jackson 
and Furner is based on a statistical formula, renders these data unsuitable for the 


present discussion. 
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mean masses in Table IV are too large. The character of the pro- 
gression is best seen from the smoothed values of the differences 
given in the last column of Table V. Considered as corrections, 
these differences are such as would reduce the masses to values 
representing the selection for t=o%06. Even at this moderate 
distance the smallest masses will still be excluded; hence some 
further correction, a constant, is required. This cannot be accu- 
rately specified; but in general it may be said that, if we regard 
the selection for the two nearest groups in Table V as approximately 
representative of the visual binaries as a whole, the data of this table 
would indicate that the mean masses of the binaries actually 
observed at distances corresponding to t=o"%o2 are roughly 50 
per cent larger than those of a representative collection. 

This result at once suggests that the variation of log wu with 
spectral type found above has been influenced by selection, for in 
general the mean parallax of the earlier-type stars, which are all 
highly luminous, is less than that of the late-type dwarfs. Forming 
the mean parallax for each of the spectral groups in the first half 
of Table II, we find in fact a progression in distance as shown in 
Table VI. If the data of Table V may be accepted, the values of 


TABLE VI 
CORRECTIONS TO LOG yu 


(To Eliminate Selection in Mass) 











Sp. ts No. Stars ~~ to 
Rite a Dds ee (07133) 3 —0.03 
Re errr 0.051 14 —0.02 
_. EE he rere 0.038 13 —0.05 
BS ba cia wamsamuiee 0.055 18 —0O.O1 
Ra eat Se ae ee 0.069 10 0.00 
Deca cwidksVoKwewmer O.114 7 +0.06 
ea re 0.229 4 +0o.1I 














log « in Table IV must be corrected as indicated in the last column 
of Table VI, in order to obtain the true variation of log uw with 
spectral type. To make the results fully representative, the 
additional constant correction referred to above must also be 
applied. This seems to be of the order of —0.08. 
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Unfortunately, however, the results in Table V are not as 
reliable as could be wished. The individual values of the masses 
are affected by large uncertainties, and the number of stars in 
the separate groups is small. Further, the residuals log u—log yu 
depend upon the ordinates of the curve for AM _ shown in Figure 1. 
Any error in the curve, in the nature of a tilt with respect to the 
axis, must therefore have introduced a progression in the residuals 
which only partially disappears when they are grouped according 
to parallax, because such a grouping, in part at least, is also a 
grouping with respect to spectral type. This uncertainty is 
probably not very serious, for Table V shows that the relative 
effect on the corrections to log u, for the two extremes of distance, 
is the differential error in the ordinates of the AM curve for F2 
and G8. Since these ordinates are not widely separated, the 
influence of incorrect slope in the curve cannot be great. 


6. ADOPTED MASSES 

On the basis of the evidence in Table VI, the corrections given 
in the sixth column of Table IV were finally adopted. These are 
the smoothed values of the results in Table VI, extrapolated to 
Bo and slightly modified in the case of type K because of the small 
number of these stars. The small constant correction referred to 
above is ignored, because of its uncertainty and relative unimpor- 
tance. It merges with the uncertainty involved in the assumptions 
of comparability which must be made in extending the results for 
the masses of the binaries to the stars at large. 

The differential correction, affecting the relative values of the 
masses as it does, is of more interest. The present determination 
is to be considered as only a first approximation; the results are at 
least in the right direction, and it is believed that they do not 
represent an over-correction. 

We have, finally, as adopted values of the geometrical mean 
masses of the visual binaries, the data in the seventh and eighth 
columns of Table IV. These results refer mainly to dwarf stars 
and correspond to the mean absolute magnitudes given in the 
second column of the table. 
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Owing to the small number of late-type giants in the lists of 
visual binaries, it is not possible to use the statistical method for 
a determination of the masses of these stars. Some indication 
of their values is afforded, however, by the application of equation 
(4) to such of these objects as occur in the spectroscopic list. For 
14 Gand K stars, with a mean absolute magnitude of 0.0, 


log w=0 .47+0.16 


whence 4 =3.0 with an uncertainty of about 4o per cent. 

The use of the geometrical mean makes it difficult to compare 
the masses thus found with other results,’ for which the arith- 
metical mean has usually been employed. For a Gaussian distribu- 
tion of the logarithms of the mass, which seems to be that indicated 
by the available data, the conversion formula would be? 


log w=log w+2.53 7? 


where yu is the arithmetical mean, and r the probable dispersion, 
i.e., the median value of the differences log u—log uw, arranged in 
order without regard to sign. The value of ris difficult to determine. 
For the dwarf stars it is certainly small, probably not larger than 
0.15 (see page 185), whence, approximately, w=1.14 uy. 

Applied to the present data, this relation would lead to arith- 
metical mean masses a little larger than those previously found, 
but close comparisons are not possible. The change in mass with 
spectral type is of the general character of that resulting from 
other investigations and is certainly well established, although the 
variation given here is not as large as that derived by Ludendorff? 
from a discussion of spectroscopic binaries. 

Strictly speaking, the present results refer only to the particular 
group and class of stars upon which the conclusions are based. 

‘ Russell, loc. cit.: Publications of the American Astronomical Association, 3, 327; 
Popular Astronomy, 25, 666, 1917. Aitken, The Binary Stars, pp. 206 ff., 1918. Van 
Maanen, Publications of the Astronomical Society of the Pacific, 31, 231, 1919. Luden- 
dorff, Astronomische Nachrichten, 189, 145, 1911; 211, 105, 1920. 


2 Adams, Strémberg, and Joy, Mt. Wilson Contr., No. 210, p. 17; Astrophysical 
Journal, 54, 25, 1921. Although derived for use in another connection the formula 


is generally applicable. 


3 Loc. cit. 
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Their applicability to any other class of stars presupposes some 





assumption. Thus in passing from the combined masses of the 
visual binaries to those of single stars, we can make use of the 
average mass-ratio of the two components and associate the result- 
ing average values for the central stars with that for the non- 
binaries of the same spectra; but it is by no means certain that the 
stars so associated are strictly comparable. For the further dis- 
cussion we make this assumption. Presumably the uncertainty 
is not large, for we thus arrive at a mean mass for the solar-type 
stars which is the same as that of the sun itself. To give this con- 
clusion force it must be added that the average deviation of the 
mass of a dwarf star from the mean mass for the type in question is 
small (see next section). 

There is some evidence that the mass-ratio varies with the type. 
Adopting, however, the mean value 0.75 for the ratio of secondary 
to principal component, we find for the stars at large the geo- 
metrical mean masses (J) in the last two columns of Table IV. 


7. DISPERSION IN MASS 


The residuals in log uw for the sixty-nine stars discussed above 
throw some light on the dispersion in mass. Freeing the individual 
residuals from the influence of selection in mass, we find the fre- 
quencies shown in Table VII. In spite of the small numbers of 





TABLE VII 
DISTRIBUTION OF A LoGu 



































| 
| No. RESIDUALS TOTAL 
INTERVAL OF ALOG yw | 
| + | - | Obsd | Theor 
©.00-0.00...... | 8 4 12 17 
©. 10-0. 80... 2.2. ; 10 8 18 15 
8 ee ee 6 7 13 12 
| ne 5 7 | 12 10 
SS er 5 5 Io 6 
fk eee ° I I 4 
©:0B-O.660. 2. ccecss I I 2 3 
ee, ° ° ° I 
0.80-0.89.... I ° I ° 
era 36 | 33 69 | 68 





? Aitken, The Binary Stars, pp. 206, 216, 1918. 

















THE MASSES AND DENSITIES OF THE STARS 185 


stars, it is seen from the second and third columns that the frequency 
curve is approximately symmetrical. Further, the last two 
columns show that, within the uncertainty of the data, the values 
of log uw follow a Gaussian distribution. 

The probable dispersion (dwarfs Fo to M) is 


r(A log w) = +0.22. (10) 


This, however, is not a true measure of the dispersion in log p, 
for the result has been derived with the aid of equation (4) and 
includes, therefore, the influence of errors in M and M,. Those 
in M, arise from the factor a3/P? and are relatively small. Bearing 
on the accuracy of M are the data in Contribution No. 199, p. 15, 
from which for the spectral types and luminosities involved, 
r(M)=+0.36 mag. By (4), the contribution to the calculated 
dispersion in log yw arising from errors in M is therefore +o.22, 
which is the same as the calculated value of r(A log y) itself. Since 
the dispersion in the masses presumably is not zero, the adopted 
value of r(M) is probably too large. This unfortunately makes it 
impossible to gain any definite information as to the dispersion, 
and we can only indicate, as in the second column of Table VIII, 


TABLE VIII 


PROBABLE DISPERSION IN MASS 











r(M) ae Limits for Prob. =4 
0.35 +0.07 0.85u-1.184 
0.30 0.13 0.74 -1.35 
0.25 ©.10 0.09 -1.45 
©.20 0.18 0.66 -1.51 

+0.15 +0.20 0.63 -1.58 











the probable values corresponding to various assumed values of 
r(M). The last column of the table gives the limits within which 
it is an even chance that the mass of any visual binary will be 
found. Thus, if we assume that the true probable error of the 
spectroscopic absolute magnitude is +0.25, the masses of one-half 
of the stars should be between o.7 and 1.5 times the mean mass. 
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Apparently the range in the masses of the dwarf stars from 
Fo on must be small. Further evidence is given on page 214. 
These results are in general agreement with the earlier findings of 
Russell and Eddington.* 


8. EQUIPARTITION OF ENERGY 


A possible analogy between the dynamical behavior of the 
stars and that of the molecules of a gas has frequently been men- 
tioned.2, Applicability of the kinetic theory would imply an 
equipartition of the energy of translation such that 


ee see = MnUn (11) 


where the individual values in v;, include all the velocities associated 
with a definite class of masses u,. In a gas this final state is 
brought about very quickly; in a collection of stars the period of 
relaxation, if the term has any meaning at all, obviously must be 
enormously long. 

As a matter of fact, the theory for a collection of stars is much 
simpler than that of a gas, because the occurrence of collisions 
and close encounters is so rare that they may be neglected.’ The 
effective agency for the transfer of energy is therefore only that 
arising from the attraction of the great mass of stars as a whole 
upon the individual members of the group. But this simplification 
of the theory and the consequent restriction of the mechanism for 
the redistribution of energy so enormously lengthens the interval 
required for the attainment of statistical equilibrium that Edding- 


t Monthly Notices, 77, 30, 1916; 604, 1917. Russell, in this connection, makes 
the important suggestion that the real dispersion in mass may be largely concealed 
in a discussion based upon spectroscopic data. This would be the case if the spectral 
characteristics used for the determination of a star’s absolute magnitude are functions 
of its surface brightness (temperature) and density alone. For Russell’s discussion 
of this interesting point see the following note on page 238. 


See for example Poincaré, Hypothéses Cosmogoniques, p. 257, 1911; Halm, 
Monthly Notices, 71, 634, 1911; Eddington, Stellar Movements, pp. 159, 247, 1914; 
Jeans, Problems of Cosmogony and Stellar Dynamics, pp. 224 ff., 1919. 


3 Jeans, op. cit., p. 229. 
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ton’ rejects altogether the analogy with a gas and proceeds on the 
principle that the motions of the stars do not mutually interfere. 

As far as the present state of the universe is concerned, Jeans? 
concurs in this opinion, but points out that the conditions of an 
earlier stage in its history may have been entirely different. Fol- 
lowing his conjecture that the system has evolved from a sotating 
nebula of dimensions much less than those of its present configura- 
tion, he shows that the earlier stages of development would be 
favorable to a rapid redistribution of energy, and that the system 
may have proceeded far on the way toward equipartition before its 
gradual expansion reduced the rate of transfer to its present negli- 
gible amount. The fact that readjustments are now inappreciable 
does not therefore demonstrate the absence of equipartition, which 
must be put to an observational test. 

Although a certain type of systematic motion is not incompatible 
with a steady state, the phenomenon of star-streaming shows 
that equipartition has not been completely attained. Within 
the individual streams, however, there may be an approach to this 
state, just as in the case of two intermingling streams of gas. 
Further, if the stream motion is small as compared with the peculiar 
velocities of the stars, it will have little influence on the statistical 
averages in (11). For a decade or more some correlation of mass 
with velocity of the type hereby implied has been suspected, 
for it has been known that in general the most massive stars are 
the most luminous, and the most luminous stars have on the average 
the smallest velocities.s Now, however, we are in possession of 
data which admit of a more definite test. 

For the B and A stars, whose extreme range in luminosity is 
small, we have the determinations of mean radial velocity by 

* Eddington, Joc. cit., p. 254. See also Astronomische Nachrichten, Jubilaums- 
nummer, p. 9, 1921. 

* Op. cit., p. 237. 

3 Campbell, Stellar Motions, p. 207, 1913. Kapteyn, Mt. Wilson Conir., No. 45; 
Astrophysical Journal, 31, 258, 1910. Russell, Observatory, 37, 174, 1914. Kapteyn 


and Adams, Mt. Wilson Comm., No. 1; Proc. Nat. Acad. of Science, 1, 14, 1915. 
Adams and Strémberg, Mt. Wilson Contr., No. 131; Astrophysical Journal, 45, 


293, IQI7. 
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Campbell.* For the remaining types, within which the differences 
in luminosity are much greater, we have the recent Mount Wilson 
results,? giving both radial and space velocities for each of these 
types as a function of absolute magnitude. Since the values of 
the masses in the last column of Table IV are probably representa- 
tive of the stars at large, they may be combined with the mean 
velocities of the different types, in conformity with (11), to deter- 
mine the degree to which the condition of equipartition is satisfied. 

The values of the velocities have not been corrected for stream 
motion and thus subject the hypothesis of equipartition to a severe 
test. The Mount Wilson results for radial and space velocities 
include the A-term. For the late-type dwarfs, this is small, and 
negligible for the present discussion. Campbell, on the other 
hand, has freed his results from this term, which for the B stars is 
two-thirds of the mean radial velocity itself. If one prefers to 
regard the A-term as a real velocity, and not merely a systematic 
correction, its contribution to the energy can be found as indicated 
below. 

It has been shown that the distribution of the space velocities 
for types F to M is not in accordance with Maxwell’s law,’ but 
agrees well with the assumption that the logarithms of the velocities 
follow a Gaussian distribution. Nevertheless, Campbell’s theorem 
to the effect that the arithmetical mean space velocity equals twice 
the mean radial velocity, holds for the giants with close approxi- 
mation. We assume that a similar state of affairs exists for the 
B and A stars, and thus obtain for these types the arithmetical 
means of their space velocities. The square of the arithmetical 
mean space velocity could have been used directly for a test of 
equipartition; but to evaluate the influence of the A-term, con- 
sidered as a velocity,’it is simpler to use the mean-square space 
velocity. 

t Lick Observatory Bulletins, 6, 101, 1911; 7) 19, 1912. 

2 Adams, Strémberg, and Joy, Mt. Wilson Contr., No. 210; Astrophysical Journal, 
54, 9, 1921. These results are vital for the remainder of the discussion. Without 
them it could not have been undertaken. 

3 Kapteyn and Adams, Mt. Wilson Comm., No.1; Proc. Nat. Acad. of Sciences, ty 
14, 1915. Adams, Striémberg, and Joy, op. cit., p. 15. 

4 Ibid., p. 13. 
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For a Gaussian distribution of the logarithms of the space 
velocities, the relation of the mean square velocity to the arith- 
metical and geometrical means is expressed by* 


log v?=2 log v+ =2 log v+0.148 


I 
2h? Mod. 


I 
Po ie 2 log y+o.296 


(12) 
log v?=2 log v+ 


where / is the modulus of the distribution function. It is assumed 
that the value of / for types F to M given in Contribution No. 210 
applies also to the B and A stars. 

The mean radial velocities taken from Campbell’s papers? are 


B3, 6.2; B8.5, 6.7; Ao, 10.3; A2, 11.4; A5, 12.6. 


The corresponding mean square velocities were obtained by 
doubling the radial velocities, as indicated above, and substituting 
the results into the first of (12). For the later types the values of 
log v were derived from the data for v given in Contribution No. 210. 
These were plotted—v against corresponding M—thus giving four 
curves, which refer to the approximate mean types Fs, Gs, K3, 
and Ma+. The mean absolute magnitudes of the dwarf stars of 
these types were then interpolated from Table IV and used as argu- 
ments to read from the curves values of v corresponding to the 
average dwarf of each of the four types. Graphical interpolation 
and an extrapolation to Fo, followed by substitution into the 
second of (12), gave finally the required values of log v* for every 


«The second equation is found by applying the mean value theorem to the 
frequency function 
h Mod. —p 
. 


Tr 


F(v) = (log v—log v)? 


The first results from the combination of the second with 
I 
log v=log v - 

5 6 Tie Mod. 


which is given by Adams, Strémberg, and Joy, op. cit., p. 17. 


2 Lick Observatory Bulletins, 6, 127, 1911; 7, 28, 1912. 





Igo 


half-interval of spectral type from Fo to Ma. 
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These mean-square 


velocities are not affected by any smoothing from type to type, 


since the interpolation card was drawn through the points derived 


from the four curves giving v the function of M. 


The results, including those for log Mv’, 


Table IX. 


columns of 


are in the first six 


With the exceptions of the B stars, the 


TABLE IX 


THE ENERGY CONSTANT 

Sp VM Log 4%} Lo Log M Residual Log @ | Log 4@ Residual 
Bz —0.60 | 0.95 34 (3.29) (+0.28) | 1.74 69) 0.07) 
B8.5 +o.4 | 0.81 | 2.40 | (3.21) + .36) | 1.62 13) tr .33) 
Ao | ©.7 | 6.78 | 2.978 3.56 + .O1 2.00 78 — .02 
| ae og See eo 3.57 CO =| 2.10 2.50 = .O4 
As.. ..| 1.5 | 0.60] 2.95 | 3.55 + .o2 | 2.19 2.79 — 03 
Fo 2.4 1 0.40 | 3.11 3.81 + .o6 2.36 “6 00 
5 | 3-3 | 0-19 | 3.36] 3.55 + .02 2.56 | 2.75 r .O1 
Go | 4.4 90.909 3.62 ke - .04 2.76 | 2.75 4 OI 

| aR | 5.2 | 9.88 | 3.78 | 3.66 — .09 2.041 2.92 - .O 
Ko 5.9 | 9.83 | 3.80 3.63 - .06 | 2.84 67 + 09 
5. | 7-1 | 9.79 | 3-74 3-53 tT .O4 2.92 71 05 

| 9.81 9.77 | 3 78 | 3.55 +o.02 3.06 2.83 —0.07 
Means, average deviations ... 2.87 +0.036 | 2.76 +0.036 





values of the energy show close accordance. For a range of 10 to 1 


in the masses the average deviation in Jv’ is g per cent, and, by a 
moderate amount of smoothing, is easily reduced to 3 per cent. 
The large divergence for the B stars seems to be real. It cannot 
be supposed that the adopted values for their masses are in error 
by 50 per cent, nor does it seem permissible to assume that the 
velocities could by any chance be increased by the 2 or 3 km 
necessary to bring the energy up to the mean for the other types. 
It is more likely that the velocities are already too large, because 
of undetected spectroscopic binaries among the stars combined in 
the means. 

The inclusion of the A-term as a velocity removes only a small 
part of the divergence. Thus, for any radial direction, the resultant 
mean-square velocity r? would be 


=v+K?—2vK cos a=v?+K? 


r? 
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where a is the inclination of the space velocities to the given line of 
sight. Extending this to all radial directions, we have for the 
entire group of B3 stars, 


r?= (14 .8)?-+(4.7)?=241.1; log Mv? =3 .33, 


an increase of only 0.04 in the logarithm of the energy. 

The foregoing results depend upon the assumptions that the 
mean space velocities of the B and A stars may be found by doubling 
the mean radial velocities, and that the modulus / applies equally 
to all types. We may avoid these assumptions by restricting 
the test to the mobility in the line of sight. For this comparison 
the A-term has been included throughout. The mean radial 
velocities from Campbell’s data thus become’ 


B3, 7-45; B8.5, 6.50; Ao, 10.1; A2, 11.3; As, 12.5 km/sec. 


For the remaining types Contribution No. 210 was again used, 
the values being read from curves for radial velocity @ as a function 
of absolute magnitude as in the case of v. The values of log M6? 
(eighth column, Table [X) show about the same degree of accordance 
as those of log Jv’, the average deviations, without smoothing, 
again being 9 per cent. The inclusion of the A-term as a part of 
the radial velocity has brought the B3 stars into agreement with 
the mean for types A to Ma, and only the B8—Bo stars are now 
divergent.?, But this improved agreement is not to be taken very 
seriously, and perhaps the most satisfactory point of view is that 
since the B stars are peculiar in many particulars, the abnormality 
shown by Table IX should not be looked upon as surprising, 


t Tbid., 6, 108. The means for the B stars are those found by correcting Camp- 
bell’s values V; by 4.7 (Bo—Bs5) and 4.1 (B8-Bo). For the A stars the mean of V3, 
Lick Observatory Bulletin, 7, 20, has been used without correction. 


2 Freundlich, Physikalische Zeitschrift, 20, 561, 1919, has attempted to show that 
the K-term may be identified with the displacement of spectral lines toward the red 
required by the theory of relativity. This displacement may be involved, but results 
given later in this paper indicate that other factors must also play a part. See foot- 
note p. 212. 
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especially when it is remembered that they form a local aggregation' 
and that their masses are among the largest known. 

There are suggestions of systematic divergencies in both sets 
of residuals, but by forming the means of the two series of values 
of the energy constant, these disappear, and the average deviation 
is reduced to 5 per cent. 

Were the condition of equipartition exactly satisfied, we should 


have 


| 


2= 3, 


fos) 


For a Gaussian distribution of the radial velocities 


where @ is the square of the arithmetical mean radial velocity 
used in Table IX. For unit mass, log v? = 3.57; and since log 7/2= 
0.20, log # =2.96. The mean-square velocities therefore appear 
to have a ratio of 4 to 1 instead of 3 to 1. It is known, however, 
that the radial velocities do not follow a Gaussian distribution 
closely. 

Approaching the matter from another direction, we have for 
a random distribution in direction, irrespective of the form of the 
frequency function of the magnitudes of the velocities, 

v=20. 
Substituting into the first of equations (12), we find the condition 
log v7? =log #+0.75 

where @ denotes the square of the arithmetical mean. The differ- 
ence in the two members of this equation is 0.06 in the logarithm, 
or 11 per cent, a much better agreement than before. The remain- 
ing difference is doubtless to be attributed to the fact that the 


velocity distribution is not random; stream motion produces a 
preferential drift parallel to the galactic plane. 


t This result follows directly from an analysis of star counts. Seares, Publications 
of the Astronomical Society of the Pacific, 30, 114, 115, 1918. Shapley in recent papers 
has accumulated much additional evidence: M?. Wilson Contr., Nos. 157, 161; Astro- 
physical Journal, 49, 311, 1918, and 50, 107, 1919; Mt. Wilson Communications, 
Nos. 54, 64; Proceedings of the National Academy of Sciences, 4, 224, 1918; 5, 434, 1919. 
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In spite of the known departures from Maxwell’s law, there is, 
with the exception of the B stars, a surprisingly close approach to 
a state of equipartition. That the deviations are no larger is 
remarkable,’ and suggests that for the stars considered (dwarfs) 
there is little difference between the different spectral types in the 
average contribution of stream motion to the total energy of 
translation. Since the unit of mass—that of the sun—is 1.9 X 10% 
grams, and since the velocities are expressed in kilometers per 
second, the mean kinetic energy of translation is of the order of 
3.5 X 10% ergs. 


Q. DERIVATION OF MASS FROM THE ENERGY-CONSTANT 


The constancy of the energy is apparently such as to justify 
an attempt to evaluate the masses of stars having luminosities 
and velocities different from those appearing in Table IX. Thus 
by reading from the curves referred to above the values of v and 6 
corresponding to different values of M and combining them with 
the energy-constant, more or less hypothetical values of the corre- 
sponding masses can be obtained. 

This procedure assumes that equipartition holds for the stars 
of high luminosity as well as the dwarfs, for which alone a direct 
test has been possible. The behavior of the B stars suggests 
caution in an extension of the principle to other giants; but since 
the translatory energy of the B’s is of the same order as that found 
for the dwarfs, it is perhaps not too much to expect that the use of 
the energy constant will lead to masses for later-type giants which 
will also be of the right order of magnitude. 

The results obtained in this manner are given in Table X. The 
values of log % are the means of the two values found with the 
aid of v and @ from the two energy constants, log Mz? and log Mé. 
To indicate the agreement between the masses derived with the aid 
of the two velocities, the differences log M¢,—log J, are also given 
in the table. In general these are reasonably small, although 
for the group of Ma stars the extreme values correspond to a two-to- 
one ratio in the masses. The data of Table X express, with more 


t In this connection see Jeans, Problems of Cosmogony, p. 238, 1919. 
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detail and precision, the correlation between mass and luminosity 
found by van Maanen.' 
TABLE X 


Mass AND LUMINOSITY 


I Gs K3 Ma 
M 
og K Diff log Dif. | log %] Diff log Diff 

mee | ; ; 0.74 | —0.08 pate ; 0.95 0.30 
—2 0.78 ( 62 02 } - 62 
J 74 51 t 0d | 0.72 >. 10 47 t 10 

° 68 - 38 | + 10 45 t 05 3 24 
+1 62 tT .02 28 | + .14 28 t 14 08 r -3° 

2 51 | 0 17 rT «32 17 I 

SC eceee 30 0.00 of O7 T 5 

4 0.01 —-0.04 | 9.96 | — .02 | 0.01 } 

«. S.37 f. . 9.55 = 27 9.94 - 

_ ee Pas |e : g.dI =O. 21 9.59 — Jee 

: entree: (ee ee: Se: .-| 9.84 — .06 

Meiipala wie chain poacs re sibiemeteterenrs Qg.5E | —0.05 

PP sakdeeven 9g .05 0.14 


The general relation of these results to spectral type and absolute 
magnitude is best seen by means of the diagram in Figure 2. 
This was prepared by plotting the data in Table X and reading 
from the resulting curves the values of M _ coresponding to 
integral values of the mass J%. ‘These values of M were then 
plotted in the diagram of Figure 2, along with similar points on the 
dwarf branch, whose co-ordinates were derived from the second 
and last columns of Table IV. Through these points were drawn 
the full-line curves of the figure, which thus represent lines of 
equal mass. 

The curves are highly irregular, the maximum near Go followed 
by the minimum at Ko~-Ks5 being particularly noteworthy. Much 
of this disappears, however, when the ordinary visual absolute 
magnitude is replaced by the bolometric magnitude, which measures 
the total radiation and is therefore the quantity of physical signifi- 
cance. The reductions to bolometric magnitude (see following 
section) of interest here are 


Go, —0.09; G5, —0o.30; Ko, —0.73; K5, —1.30; Ma, —1.60 mag. 


t Publications of the Astronomical Society of the Pacific, 31, 231, 1919. 
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The application of these corrections would practically remove 
the minima from the curves for all but the largest masses, and these 
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Fic. 2.—Distribution of mass (full-line) and mean density (broken curves) 
derived from the principle of equipartition. The full-line curve running downward 


to the right is the line of the maximum frequency of the dwarfs. 








. 


196 FREDERICK H. SEARES 


are relatively uncertain because they depend more or less upon 
extrapolations of the data in Table X. 

To explain the derivation of the lines of equal density, which 
are also shown in Figure 2, it is necessary to consider the question 


of surface brightness. 
10. SURFACE BRIGHTNESS 

To derive values of surface brightness we may start from 

Stefan’s law 
E=oT‘4 (13) 

where E represents the total energy radiated per unit area of the 
surface of a star whose effective temperature is 7. Let J be the 
integrated visual intensity, referred to the same unit of surface, 
and write J =aE, whence 


2.5 log E=2.5 log J—2.5 loga 


and, by substitution into (13), 


2.5 log J—2.5 log a=2.5 log a+ 10 log 7. 


Further, let Mg and / represent the total energy and surface bright- 
ness, respectively, expressed in magnitudes. The former is the 
bolometric absolute magnitude, while = — 2.5 log (J & Area) is the 
ordinary visual absolute magnitude. Since a is also the ratio of 
the visual energy to total energy received from the entire surface 
of the star, 
—2.5loga=M-—Msz, 
and 
j=M—Ms-—10 log T+ const. (14) 

where the constant involves logo and an additional constant 
depending on the zero point chosen for Mg. Eddington’ has calcu- 
lated values of the difference between the visual and bolometric 
absolute magnitudes M/— Mz, for a range of temperatures including 
the normal classes of spectra. Through these values, which are 
based upon Nutting’s measures of visual sensibility for the normal 
eye, 7, by (14), becomes a simple function of 7. 

* Monthly Notices, 77, 605, 1917. Nutting’s measures of visual sensibility, which 
are the basis of Eddington’s calculation, are practically identical with the A.L.E. 


visibility-curve. 
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Results found by this equation may be checked by an equivalent 
formula derived by Hertzsprung in 1906. Starting from Planck’s 
law and the visibility measurements of Langley, Abney, and others, 
he found’ 


jg=+2-3()° 34 const. (15) 


in which ¢ is the second constant in Planck’s formula, for which 
the value adopted is 14300. Both (14) and (15) involve quadra- 
tures, and, starting as they do from entirely different observational 
data, independently discussed, provide an excellent numerical, con- 
trol. How close the agreement is may be seen from Table XI. 


TABLE X1 


SURFACE BRIGHTNESS AND TEMPERATURE 











T M— Mp j | JH jI-jH 
2540°. +2.59 +6.32 | +6.32 0.00 
3000 . +1.71 +4.72 +4.68 + 4 
ME a sei pus wine +0.95 +3.17 +3.16 + I 
OS eee +0.35 +1.60 +1.58 oa 2 
ere 0.00 0.00 0.00 ° 
7EO0 ...2cc0s-] FO.02 —0.95 —0.93 - 2 
GOOD ..ccccreel FO.E8 —1.64 —1.62 - 2 
EOSOD .... ccc.) SO. Zt —2.12 —?2.09 _ 3 
$9600 ....2..+--) FO.§3 —2.40 —2.45 —0.03 

















The values of 1/—Mz have been taken directly from Eddington’s 
table. The maximum difference between the results from the two 
formulae is only 0.04 mag., which is ample. 

To connect these results with spectral type, Wilsing’s spectral- 
photometric measures of the reciprocal temperature, c/T, of 199 
stars? were plotted against their respective types. The ordinates 
(c/T) of the curve thus obtained are in the second column of 
Table XII. The corresponding effective temperature, 7, and the 
surface brightness, 7, calculated by the foregoing formulae and 
expressed in magnitudes (© =o.0o), are in the third and fourth 


columns, respectively. 


t Zeitschrift fiir wissenschaftliche Photographie, 4, 43, 1906. 
2 Potsdam Publikationen, No. 74, 1919. 
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With few exceptions, the stars in Wilsing’s list are all giants 
whose mean absolute magnitudes, with sufficient approximation, 
may be assumed to be zero. To extend the results to other lumi- 
nosities, use was made of the Mount Wilson color indices,’ whose 
values are in the fifth and last columns of Table XII, those for the 


TABLE XII 


SPECTRUM, EFFECTIVE TEMPERATURE, AND SURFACE BRIGHTNESS 


GIANTS, M =o DWARFS 
Se 
/tT | 1 rey 1 7 M re 
Bo 1.36 | 10500 2.28|—0. 22 
Bs I.43 | 10000 2. 14/—-0.17 
Ao 1.55 9230 1.89} 0.00 
As 1.76 8110 1.45/+0.19 
Fo 2.04 | 7000 0.88) 0.38 
Fs 2.35 6080 0.26) 0.62] 2.35 | 6080 26/+23.2 +o. 62 
Go 2.70 5300 +o.44| 0.86) 2.48 | 577 0.00] 4.4 0.72 
Gs ; 3.10 4010 .. 2a :.¥e 60 5500 ».25 <2 9.82 
Ko 3.70 3860 2.41 1.45) 2.93 | 4550 ©.gO 5-9 99 
K5 4.37 3270 3.71 1.84) 3.47 | 4120 I. of : Ea 
Ma 4.65 3080 4.25] 1.88] 4.30 | 3330 3.58 .8 r.7( 
Mb 4.82 2960 4.58/+1.88 
M« 4.94 2890 +4.81 


dwarfs corresponding to the mean absolute magnitudes in the 
column immediately preceding. The color indices were reduced 
to Schwarzschild’s absolute system,’ whose relation to effective 
temperature is known. ‘This gave the relation between the Mount 
Wilson values and the temperature for both giants and dwarfs. 
Owing to various necessary assumptions the temperatures thus 
found for the giant stars differ slightly from the adopted values 
given in the first half of Table XII, and this systematic difference 
was made the basis for reducing the calculated temperatures for 
the dwarfs to the adopted system of Wilsing. The results, in 

t These results, which have not hitherto been published, are still to be regarded 
as provisional. They take account of the important difference in color between 
giants and dwarfs, which cannot be disregarded in any discussion of the surface 
brightness of late-type stars. See Mt. Wilson Comm., No. 59; Proc. Nat. Acad. of 
Sciences, 5) 232, 1919. 


2 Gottingen Aktinometrie, B, 29, 1912. 
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the sixth column of Table XII, were then used to compute the 
corresponding values of 7 as before." 

The difference in surface brightness between giants and dwarfs is 
large, and cannot be neglected. Unfortunately, the change in 7 
with absolute magnitude is not well determined, for this must be 
inferred from the behavior of the color indices, and the data avail- 
able at present do not extend to stars of the highest luminosity. 
This means that extrapolation will be necessary, and, in view of the 
consequent uncertainty, a simple linear relation of surface bright- 
ness to absolute magnitude has been used in the calculation of 7 
for all values of M not appearing directly in Table XII. 

The reliability of the foregoing results depends upon the possi- 
bility of representing the energy distribution in the spectra of stars 
of different types by means of the black-body radiation formula. 
It should be noted that it is not at all a question of agreement 
of the effective temperature, calculated from the radiation formula, 
with the true mean temperature of the radiating layers in the 
star’s atmosphere. In general these will not agree. The stars 
are not perfect radiators; further, general and selective absorption 
deforms the energy-curve to such an extent that effective tempera- 
ture cannot be an accurate measure of true temperature; but in 
the visible region the resultant curve seems to agree? approximately 
with a black-body curve, although not with the curve corresponding 
to the true mean temperature. That being the case, the surface 
brightness is immediately calculable by the foregoing formulae, 
provided the effective temperature corresponding to the deformed 
curve is employed for the computation; and it is of course just 
this temperature which is given by Wilsing’s observations. 

t The values of j in Table XII may be compared with Opik, Astrophysical Journal, 
44, 296, 1916; Russell, Publications of the Astronomical Society of the Pacific, 32, 307, 
1920. Related questions are discussed by Wilsing, Potsdam Publikationen, No. 76, 
1920; Astronomische Nachrichten, 214, 185, 1921. 

2 See, for example, Wilsing, Potsdam Publikationen, No. 74, pp. 17, 18. The mean 
systematic deviations of the observed values of the energy at different points in the 
visual spectrum from Planck’s formula are small. Further, Milne, Monthly Notices, 
81, 378, 1921, finds that the theoretical spectrum for non-selective absorption is 
practically identical with that corresponding to black-body radiation, the principal 
difference being a bodily displacement of the energy-curve, corresponding to a reduc- 


tion of 3 per cent in the constant of Wien’s displacement law. 
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The possibility of systematic error in the values of the temper- 
ature is more serious. ‘Those for the B stars given in Table XII 
are almost certainly too small. The early-type stars in Wilsing’s 
list are not numerous, and the dispersion in the values of c/T is 
large. In fact, the data would be almost equally well represented 
by assigning a temperature of 12,000° to Bo, which would change 
j from —2.28 to —2.67. The value for B5 would also be slightly 
changed, but that for Ao can scarcely be modified without violence 
to the data.’ 

It will be noted that 7 for a Go giant is nearly 500° lower than 
that for a dwarf of the same spectrum. Some question as to the 
reliability of Wilsing’s measures has been felt because of the low 
value found for Go stars as compared with the probable temper- 
ature of the sun. But the objects observed by him are giants; 
with proper allowance for the sun’s absolute magnitude (4.85), 
the corresponding temperature from Wilsing’s data would be about 
5820°, which is not far from the values usually assigned, and 
indeed, agrees well with that derived by Abbot? from the solar 
constant (5860°). 

A partial control of the values of 7 in Table XII is afforded by 
Pease’s recent measures of the angular diameters of stars with the 
interferometer. In this connection, and for later use, it is con- 
venient to insert at this point the following formulae, all easily 
derived, expressing relations between the various physical constants 
of a star: angular diameter (d), linear diameter (D), mass (JA), 
density (p), apparent (m) and absolute (7) magnitude, and surface 


brightness (7). 


log d=o0.2 (j—m)—2.061 (16) 
log D=o.2 (j—m)—log r—0.030 (17) 
log D=0.2 (J—M)+0.970 (18) 

D=107.4d/r (19) 

J=5 log d+m-+ 10.30 (20) 
log p=log M—3 log R+0.14 (21) 
log p=log M+0.6 (M—j)—2.77 (22) 


t See also Rosenberg, Abhandlungen der Kaiserl. Leop.-Karol. deutschen Akademie 
der Naturforscher, 150, No. 2, 1914. 
Op. cit., D. II4. 
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These formulae give d in seconds of arc, D and J% in terms of 
the solar diameter and mass, and 7 in magnitudes, the value for 
the sun being 0; the unit for p is the density of water. The 
numerical constants depend upon the following data for the sun: 
D=865,000 miles, d=1920”, distance from earth 92,930,000 miles, 
m= — 26.72, M=+4.85. Equations (16) and (17) are the equiva- 
lent of formulae already given by Russell." 

Proceeding now to a comparison with the results of Pease’s 
measures, we have the data in Table XIII. The “cbserved” 

TABLE XIII 


COMPARISON WITH MEASURED ANGULAR DIAMETERS* 








a Orionis a Bodtis | a Scorpii 











ak taitie bent nakses a Ma Ko Map 
EG HEE SP hak. | -+0.9 +o.2 +1.2 
uty is a eRe a we ail Area y | ©7020 07095 070085 
ae rane cielo a ack ra —2.6 +o.1 —4.2 
eS ee ee ree | 0%047 07022 07040 
Obs. 7, equation (20)....... +4.6 +2.2 +4.5 
ee 8 eee +4.4 +2.4 +4.5 
oy Re eae ee ee ae —0.2 0.0 
Obs. D, equation (19)...... 252 25 506 
Cal. D, equation (18)....... 235 27 513 
Approx. J (Fig. 3)........ 8 3 5 
Approx. log p (Fig. 3)...... —6.2 —3.8 —6.6 





* Michelson and Pease, Mt. Wilson Contr., No. 203; Astrophysical Journal, 
53, 249, 1921. Pease, Publications Astronomical Society of the Pacific, 33, 171, 
204, 1921. The parallax of a Scorpii assumes it to be a member of the Scorpius 
group. The other two parallaxes are weighted mean values. 


values of 7 derived from the angular diameter by (20) are in the 
sixth line of the table. Immediately below are the calculated 
results interpolated from Table XII. It seems unlikely that the 
mean systematic difference for O—C will exceed one- or two-tenths 
of a magnitude. The effect of darkening at the limb has not been 
taken into account, either in the measures of d or in the calculation 
of 7. Nor has the question of a possible differential effect from 
this source upon the observed and calculated values of 7 been con- 
sidered, for it is evident that this cannot be large. 

It is of interest, further, to compare the values for the linear 
diameters calculated from the measured angular diameter by (19) 
with those derived from (17) or (18) with the aid of the adopted 


* Publications of the Astronomical Society of the Pacific, 32, 307, 1920. 
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values of 7. The results are in the third and fourth lines from the 
bottom of table. On the whole, Pease’s measures afford an excellent 
confirmation of the surface brightness found for late-type stars of 
high luminosity, and at present there seems to be no observational 
evidence which would justify a modification of the values adopted.’ 


It. DISTRIBUTION OF DENSITY 


The densities are now easily calculated by equation (22) with 
the aid of values of 7 from Table XII. The results for stars on the 
dwarf branch, including the B’s, are in the fourth column of Table 
XIV, which also gives values for the mean linear diameters D found 


by (18). 
TABLE XIV 


LUMINOSITY, Mass, MEAN DENSITY, AND DIAMETER 


(Units for 4 and D are values for sun; for p, that of water) 








DWARF BRANCH Giants, M =o GIANTs, J =10 
oP. 

M M. p D M p D M ? D 
Bo —1.60 | 10 0.045| 6.8 3.2 —1I.2 >.08 6 
Bs —0.20 8.3 0.20 | 3.8 8.2 >. 26 iz —1.8 03 8 
Ao ©.70 6.0 0.36 | 2.8 7.0 | 0.16 3-9 | —2.4 | 0.008 12 
As +1.50 4.0 40 | 2.4 5.6 | 0.071 4.8 | —3.1 | 0.002 | 20 
Fo +2.40 2.5 0.40 | 2.0 4.3 | 0.025 6.2 —3.6 | 0.0004 | 32 
Fs +3.35 I.5 0.39 | 1.8 3.2 | 0.0078 8.3 —4.0 | 0.0001 | 52 
Go 4.35 1.0 65 1.20 2.0 | 0.0025 Ir.5 —3.9 0.00002 } 33 
Gs +5.20 °.7¢ Ss 0.90 2.8 ©.00057 17 —3.3 0. OOOO! 102 
Ko +6.00 0.68 | 1.3 0.89 | 3.0 | 0.00018 | 28 —3.§ | 0.00 | 107 
Ks 7.20 o.62 '.6 0.83 | 2.6 | 0.000026 | 51 —2.2 | | 182 
Ma +9.50 | °.59 4 0.54 2.0 0 .0000096 | 60 at *) | 288 


For the determination of the equal-density lines (dotted) shown 
in Figure 2, values of M were read from each mass line in the 
diagram for every half-interval of spectrum. The use of (22) then 
gave the values of log p corresponding to each of these points in the 
sequence of spectral type. Plots of these results—log p against 
spectrum—then led to a series of curves, one for each integral 
value of the mass (also one for 4%=1.5), from which were read 
spectral types corresponding to integral vaiues of log p. The 


t Mr. Pease has kindly placed at my disposal provisional results for three other 
stars. The observed and theoretical values of j and d, or D, are here also in excellent 
agreement, although the accidental deviations are larger than those shown in Table 
XIII. 





~ 
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corresponding points were marked on the diagram and connected 
by the dotted lines shown in the figure. 

Within the limits considered, and with the exception of that for 
log p=-—1, the curves of equal density are sensibly linear; the 
regularity of their spacing is remarkable. Including the line for 
log p=—1 from As5 onward through the later types, the mean 
density corresponding to any given spectrum and absolute magni- 
tude is represented by the formula 


log p= —1.22S+o0.57 M+1.11) 
(S=o for Bo, 1 for Ao, etc.) j 


with a probable error of about 15 per cent, over a range in p extend- 
ing from 0.1 tO 0.000001. 

The question now arises as to what means can be found for 
testing the validity of these results. Some control is obviously 
necessary, for with the exception of the masses and densities along 
the dwarf branch, the whole scheme of mass and density distribu- 
tion depends upon the assumed applicability of the principle of 
equipartition. In this connection, the eclipsing variables, many 
of which are known to be of exceptionally low density, immediately 
suggest themselves as suitable for a test. The weak point, however, 
is that we do not know their absolute magnitudes, except in a 
very few cases. Any attempt to use hypothetical absolute magni- 
tudes fails, because it involves an assumption for the mass, and 
this, as shown by (22), is equivalent to a direct assumption for 
the density. 

Some general indications, however, can be obtained. With 
slight revisions of spectral classification, Shapley’s data" for the 
densities of eclipsing variables give the following mean values: 


Sp. p No. Stars 
B6 0.12 15 
Ao+ 0.20 48 
F (dwarfs) 0.35 9 


For these groups of B and A stars the values of log p are —o.9 and 
—o.7, respectively. From Figure 2 the corresponding values of M 
are seen to be —o.6 and +0.2. -These are of the order to be 


* Contributions from Princeton University Observatory, No. 3, p. 124, 1915. 
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expected, and are obviously possible values. The value p=o.35 
for the F-type dwarfs agrees closely with that found here for Fs, 
namely, 0.39. For the dwarfs the dispersion in absolute magnitude 
is so small that lack of knowledge of the luminosities gives rise to 
no appreciable error. 

These results are all satisfactory, as far as they go, but they 
do not reach that part of the diagram which is involved in the 
greatest doubt, namely the region of very high luminosity and 
exceptionally low density. On the other hand, if the validity of 
the mass and density distribution can once be established, the 
densities of the eclipsing variables should give a reliable determi- 
nation of their distances. Fortunately, the uncertainty can be 
removed to a large extent with the aid of the Cepheid variables. 


I2. REVISION OF MASSES AND DENSITIES BY MEANS OF 
CEPHEID VARIABLES 


There are still difficulties connected with the pulsation hypothe- 
sis as a means of explaining Cepheid variation, although there 
seems little doubt that, in some form or other, it must be used to 
account for the changes in brightness of these stars. That being 
the case, the mean density must vary inversely as the square of 
the period. Hence’ 

log p= —2 log P+ const. (24) 


Further, the period-luminosity relation? 


M=f (P) ( 


i) 
wn 


t The exact form of the relation of density to period as applied to Cepheids is 
possibly open to some question. Emden, Gaskugeln, p. 451, 1907, finds for ellipsoidal 
oscillations of a polytropic gas sphere 


Pa«(m p)7! 


where m is the order of the harmonic oscillation; the important vibration corresponds 
to m=2. Eddington, Monthly Notices, 79, 2, 1918; 177, 1919, discusses pulsations 
which are symmetrical about the star’s center and finds a result equivalent to (24), 
with the exception that the constant is replaced by a slowly varying function of the 
specific heats and the ratio of radiation-pressure to gravitation, the latter increasing 
with the mass. For a first approximation we may adopt (24). 


2 Shapley, Mt. Wilson Contr., No. 151; Astrophysical Journal, 48, 89, 1918. 
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determines the absolute magnitude. With M and the spectral 
type as arguments, log p and % can be interpolated from the 
mass-density diagram. Then, 

a) The interpolated values of log p may be compared with 
those calculated by (24), the constant being adjusted to give the 
best representation of the densities of the entire group of stars. 
Since the range in p as determined by (24) is from about 0.1 to 
0.00001, the control on the densities should be effective. 

b) Having determined the constant in (24), we may write (22 
in the form 


log M=—2 log P—o.6 (M—j)+ const. (26) 


which expresses the mass of a Cepheid variable as a function of its 
period and surface brightness, for, by (25), M is a function of the 
period alone. The values of % interpolated from the diagram may 
be compared with those calculated by (26). This checks the 
distribution of mass. 

c) Finally, with spectral type and log p calculated from (24), 
we may interpolate M from the diagram. The absolute magnitudes 
thus found may then be compared with the period-luminosity 
relation (25). This is not an independent test, but exhibits the 
results of (a) in a different form. 

Because of the enormous range in the density, the application 
of these tests will mean little except in the case of stars of accu- 
rately determined spectra. This makes trouble at once, for the 
type of a Cepheid, in so far as determined from the hydrogen lines 
at least, varies with its brightness, as was first clearly shown by 
Pease’s' observations of RS Boétis. Isolated determinations of 
spectrum based on the appearance of these lines are therefore 
unsuitable for the comparisons. Adams and Joy,’? on the other 
hand, have shown that many spectral characteristics change little 
or not at all as the light of the star varies; their ‘“estimated”’ 
spectral classes are practically independent of phase. 


* Publications of the Astronomical Society of the Pacific, 26, 256, 1914. 


2 Mt. Wilson Comm., No. 53; Proceedings of the National Academy of Sciences, 4, 


129, 1918. 
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The list of spectroscopic parallaxes' includes twenty-nine 
Cepheids for which estimated types are available. One of these, 
RY Boétis, although usually classed as a Cepheid, does not show 
Cepheid spectral characteristics and has an absolute magnitude 
at variance with that derived from the period-luminosity relation. 
Further, three cluster-type Cepheids, SU Aurigae (8.7), SU Draconis 
(9.2), and SW Draconis (10.0), and the long-period Cepheid 
V Vulpeculae (8.6) are so faint that the types assigned are less 
certain than the others. The classification of faint stars, observed 
with low dispersion, is likely to be influenced to some extent by 
the hydrogen lines. The case of V Vulpeculae, moreover, is 
peculiar. The light-curve is irregular; and the period-luminosity 
relation leads to a value of M which not only differs widely from 
that derived from the spectrum, but falls outside the limits of the 
mass-density diagram. ‘These five stars were therefore excluded,’ 
leaving a total of twenty-four. 

These data are supplemented by Shapley’s* results for twenty 
Cepheids, based on systematic observations of spectral characteris 
tics that do vary with the changes in light. The median value of the 
type, which was used for the discussion, is about 0.3 of a spectral 
interval earlier than the estimated spectrum. 

This difference will lead to different values of the constants 
in the foregoing formulae, unless the value of WM and the spectral 
type are referred to the same phase. Thus, the estimated type 
agrees with the hydrogen-line type at, or very near, the time of 
minimum brightness (estimated minus hydrogen-line type at 
minimum = —o.1 spectral interval). Hence for the comparisons, 
the minimum, and not the median M (the period-luminosity relation 
gives the latter), should be associated with the estimated spectra. 
The simplest procedure, however, is to use the median M and 
subtract 0.3 of an interval from the estimated spectra before 

* Mt. Wilson Contr., No. 199; Astrophysical Journal, 53, 13, 1921. 


?It was not noted until after the comparisons had been finished that, on the 
criterion of brightness, W Serpentis (9.0) should also have been excluded. 

3 Mt. Wilson Contr., No. 124; Astrophysical Journal, 44, 273, 1916. In general 
the means of the limiting values of the type as given by Shapley have been adopted. 
Deviations from this rule occur in a few cases where the distribution of the observa- 


tions is such that the mean does not correspond to median brightness of the variable. 
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making interpolations from the diagram. Sixteen stars are com- 
mon to the two lists, but owing to the different methods of classi- 
fication used, the duplicate observations have been treated as 
separate stars. 

The results of the application of tests (a), (6), and (c) are 
summarized in Table XV, which gives the range in log p, log %, 
TABLE XV 
REPRESENTATION OF MASSES AND DENSITIES 
EQUATIONS (24) TO (26) AND FIGURE 2 














AVERAGE DIFFERENCE 
TEST RANGE 
Est. Sp. H.L. Sp. 
(a) logp...... —1.1 to —4.8 0.15 +0.18 
(b) log M..... 0.46to 1.28 0.15 0.16 
(c) M........] —0.3 to —4.7 +0.30 +0.32 

















and M over which the comparisons were made, and the average 
differences resulting therefrom. Thus, on the basis of estimated 
spectra, formula (24) represents the density distribution in Figure 2 
over a range in log p from —1.1 to —-4.8, with an average difference 
of 0.15. Similarly, the values of M interpolated from the diagram 
in the interval —o.3 to —4.7 agree with those from the period- 
luminosity relation within +o.30 mag. We conclude at once 
that the masses and densities derived with the aid of the principle 
of equipartition are of the right order of magnitude, at least as far 
as type G3, which is the latest occurring among the Cepheids used. 

A classification of the differences in log p and log J% according 
to M, which is also more or less a classification according to spec- 
trum, revealed, however, the presence of systematic divergences. 
These run from about —o0.18 (at M=—o.7, A8) to +0.18 (at 
M = —3.8, F8) in the logarithms of both quantities. The algebraic 
signs mean that the range of mass in Figure 2, for a given interval 
in M, is too small, that the masses of the most luminous Cepheids 
are larger, and those of less luminous stars smaller, than the values 
given by the diagram. The density lines, on the other hand, are 
too closely spaced. 

These differences were made the ‘basis of a readjustment of the 
The percentage divergence in p and J% 


lines of mass and density. 
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is the same, as it must be by (22), but since the scale for M% is much 
larger than that for p (in passing from one curve to the next the 
mass changes by one unit while p changes tenfold), the adjustments 
were made in the mass-curves, the densities being computed by 
(22) as before. The masses and densities along the dwarf branch, 
of course, remain unchanged. 

The result of several trials is the diagram shown in Figure 3. 
A repetition of tests (a), (b), and (c) gave the differences shown in 
the last six columns of Table XVI, the first column under each of 
the headings log p, log %, and M referring to estimated spectra, the 
second to those based on the hydrogen lines. The average differ- 
ences have been reduced to about two-thirds their original values, 
while the run in the differences has disappeared. ‘Thus for four 
groups of eleven stars each, the mean systematic values are: 








M —1.0 —2.1 —2.7 —3.6 
Sp Fo F6 F7 F8 

log p - —0.01 —0.03 +0.02 +0.02 
log M ; —0.02 —0.O1 —0.04 0.00 





The comparison of the interpolated values of M with the period- 
luminosity relation is shown graphically in Figure 4. 

In this connection it will be noted that values of M interpolated 
from the diagram with log p and the type as arguments are sensitive 
to uncertainties in the spectra. The relations are such that the 
average difference 0.2 in M (Table XVI) corresponds to 0.1 of a 
spectral interval (see equation [23]). The normal uncertainty of 
classification therefore accounts for nearly all the deviation shown 
by individual stars in Figure 4. 

The principal change in the diagram has been the depression 
of the mass-curves near Go and the change in their curvature in 
the region of the early types. Since the data for the Cepheids do 
not extend beyond G3, the lines through types G to M have been 
left much as in the original diagram. As will be seen later, there is 
some evidence favoring the reality of the minimum near Ko. 
If the reductions to bolometric magnitude given on page 194 be 
applied, the variation of mass with spectrum becomes fairly 
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regular, the course of the topmost line (4¢=10) being indicated by 


the dotted curve in Figure 3. 
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Fic. 3.—Distribution of mass and mean density, revised with the aid of Cepheid 
variables. The dotted curve at the top indicates the reduction to bolometric absolute 


magnitude for late-type stars of mass ro. 
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The curves of equal density are little changed. They still 
remain sensibly linear, but are no longer so closely parallel. 
Formula (23) has not wholly lost its applicability, for it still 
expresses approximately the relation of density to spectrum and 
absolute magnitude, though less accurately than before. 

It will be noted that the densities along the dwarf branch are 
consistent with the values demanded by the density-distribution 
defined by the lines log p=o, —1, —2, etc. In other words, if we 


-05 LoqgP 0.0 +0.5 +1.0 +1.5 
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Fic. 4.—Absolute magnitudes of Cepheids interpolated from Fig. 3 with spectral 
type and log p from equation (24) as arguments. Circles, ‘‘estimated” spectra by 
Adams; points, hydrogen-line spectra by Shapley. The curve is Shapley’s period- 


luminosity relation. 


interpolate, say, for the value log p=—o.4 (p=o0.4) along the 
line J6=6, we find the point where it belongs, namely, on the 
dwarf branch. ‘This detail is of some importance because it shows 
that the constants in the formulae for log p and log % have been 
correctly determined. The spacing of the lines is fixed by the 
periods of the variables; their position as a group, by the condition 
that they must be consistent with the densities adopted for the 
dwarfs, whose values are independent of any assumption as to 
equipartition. 

Presumably the diagram in Figure 3 is to be preferred to that 
in Figure 2. Certainly that is the case when Cepheid variables 


a. 


are concerned; and it seems likely that this holds also for the 
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stars at large. The mass lines in Figure 2 pass through the points 
defined by the energy-constant and the velocities; no allowance for 
observational uncertainty by smoothing between the types was 
attempted; and with the exception of the terminal points on the 
dwarf branch everything preceding F5 is an extrapolation. On the 
other hand, an assumption of comparability is here involved, and 
it is perhaps worth remarking that if we reject the original diagram, 
we may thereby exclude relations between mass, density, lumi- 
nosity, temperature, etc., which differentiate stars of constant 
light from Cepheid variables. 

In view of the tentative character of the results, a detailed 
tabulation of the mass and density relations is unnecessary; 
Figure 3 will answer all practical purposes to which the results 
can be adapted. The data for the dwarf branch given in the first 
part of Table XIV, which are the most reliable, have already been 
referred to. To these have been added, by way of numerical 
illustration, results for the giants of zero absolute magnitude and 
for those whose masses are equal to ro.’ 

As a further illustration, the masses and densities corresponding 
to the spectra and absolute magnitudes of the individual stars 
measured by Pease with the interferometer have been given in the 
last lines of Table XIII. As a source of error, the values of M 
for these stars are probably as serious in respect to mass and density 
as the uncertainties of Figure 3 itself. Equation (22) shows that 
an error of 0.5 in M means an uncertainty factor of 2 in the density, 
whereas it seems unlikely that values from the diagram should be 
in error by so large an amount. In the case of an individual star, 
dispersion may of course enter, since the diagram gives mean 

t Since the displacement of spectral lines required by relativity is proportional 
to 4/R (R=radius of star) the results in Table XIV have a direct bearing on this 
question. The values of J4¢/R do not run parallel to those of the K-term. Hence 
this term cannot be wholly identified with the relativity displacement, although the 
latter may of course be involved. As an illustration, the K-term for the B stars is 


about four times that for the A’s. The ratio of the displacements on the other hand 


would be about 1.3. The most extreme deviation, however, is found by comparing 
the B’s with the giant M stars. Their K-terms are of the same order of magnitude, 
whereas the red-displacement of the lines in M stars should be of the order of 0.01 


that of the B’s. See footnote No. 2, p. ror. 
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values; but the probability of the occurrence of large deviations 
from the mean is small. 


13. ON DISPERSION IN THE MASSES OF CEPHEIDS 

From Shapley’s work it appears that, for the variables in 
globular clusters at least, the deviations in M from the values 
determined by the period-luminosity relation are of the same order 
as the observational errors in the apparent magnitudes,’ which 
of course must be used for any comparison. Practically, therefore, 
M is uniquely determined by P. Hence, by (26), the dispersion 
in the mass is related to that arising from such differences in 7 
as may be associated with the same value of P, and it is an observa- 
tional fact that the latter are small. 

Direct evidence is also afforded by the comparison of the values 
of % for the Cepheids calculated by (26) with those interpolated 
from Figure 3. The average difference in log J%% is 0.10, corre- 
sponding to about 25 per cent in 4%. But what part of this may 
be attributed to uncertainty in the spectrum ? 

With proper allowance for the systematic difference of 0.3 
between the hydrogen line and estimated spectra, the internal 
agreement is such as to indicate an average uncertainty of 0.05 
of an interval. Assume an adopted value of the type to be too 
late by this amount. From (26) 


¢ 


Mod. an = +0.6 dj. (27) 


For the region Fo-Go, j increases about 0.07 for the assumed 
increase in the spectrum. Hence, the calculated % will be 10 
per cent too large. On the other hand, we find from Figure 3, 
for F5, M = —2.5, as an average, that the same error in the spectrum 
leads to an interpolated % about 3 per cent too small. The 
resulting errors are in opposite directions; whence, in the difference, 
the uncertainty becomes 13 per cent. 

This average percentage error, combined with the dispersion 
and the errors of interpolation, must equal the average deviation 


* Mt. Wilson Contr., No. 151, p. 26; Astrophysical Journal, 48, 114, 1918. 
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of 25 per cent from the mean masses of the diagram. Neglecting 
the interpolation errors, whose inclusion would make the dispersion 
still smaller, we find the probable dispersion r(J¢)=19 per cent. 
As far as the evidence goes, we should therefore expect one-half 
the Cepheids to have masses lying between 0.8 and 1.2 times the 
mean mass. Since the Cepheids are a very special class of objects, 
it would be too much to assume that this result holds for all classes 
of stars (see Section 7). 


I4. MASSES AND DENSITIES OF CEPHEID VARIABLES 


The calculated values of the masses and densities used for the 
comparisons with the diagram of Figure 3 have been collected 
in the sixth and seventh columns of Table XVI, while their 
deviations from the adopted distribution are in the following 
columns. The formulae used are: 


log p= —2 log P—1.90 (29) 


log M=—2 log P—o.6 (M—j)+0.87 (30) 


in which the constants are the means for the two sets of data, 
the individual values differing by only 0.02. In case the spectral 
type is the median value derived from the hydrogen lines, the 
median value of M, given directly by the period luminosity rela- 
tion, is to be used. For estimated spectra, the minimum WM is 
required, although it is simpler to subtract 0.3 from the estimated 
spectrum and use the median M. 

The values of both density and mass depend upon the assumed 
constancy of the numerical term in (29), but this assumption 
apparently leads to results of the right order, although for stars 
of short period the percentage deviation may be considerable. 

For individual stars the calculated values are probably to be 
preferred to those interpolated from Figure 3. This follows not 
merely because of the dispersion in mass, but also because the 
differences for stars having determinations based on both estimated 
and hydrogen-line spectra are almost always of the same sign. 
This indicates real deviations from the mean distribution, in addi- 
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tion to those arising from the unavoidable errors of spectral classi- 
fication. In such cases the values from (29) and (30) are likely to 
be nearer the truth than those derived from the diagram. 

The results in Table XVI, both calculated and interpolated, 
may be compared with the masses and densities derived by Edding- 
ton' for certain Cepheids on the basis of his theory of radiative 
equilibrium. The agreement shown by Table XVII is remarkable 


TABLE XVII 


COMPARISON WITH EDDINGTON’S RESULTS FOR CEPHEID VARIABLES 


























Mass MEAN DENSITY 
STAR | ‘ . : —_ ey 

i Eq. (30)| Fig. 3 | Eddington | Eq. (29) Fig. 3 
, Gi ae ee | 13.0] 14.4 | 10 0.000080 |0.000044 |0.000032 
¢ Geminorum............... | 7.4 | 6.0 7.9 | .00046 .O00T 2 .OO0O014 
Sey | 6.1 | 7.8 | 6.9] .coo1g .00018 .00013 
We I os 6 Gv ccasacnnion | =< 5.4 6.5 | .coo4! .00022 .00025 
IN fais ache Woes cows | 5-2] 4.5 | 6.2] .00033 | .cc024 .00028 
eS eee eae | 5.1 | 4.7 | 6.4] .00025 | .00025 .00038 
eee 4.3 | 9.4 5-5 | .00027 00038 .00021 
NE ote "vhs Jo Wesaye ie Sesh 4.0 6.2 5.3 | .00052 .00044 .00045 
(ig es ee 3.6 5.1 5.1 | .00072 .00063 .0007 I 
I fe ot aw sae ateo's | 3.8 | 3.4 | 5.0] .0o160 .00083 ,OO13 
ee 3-0 | 4.9 4.8 | .oor10 . 0000 I .OO10 
i er 2.8] 4.5 4.5 | .0o15 .0013 .OOI! 
SU Cassiopeiae........ 2.31 3.0 5.0 | .0028 .0033 0045 
8 RE ee ere r.7 | 6.4 5.6 jo O15 0.038 0.033 














when it is considered that the two methods approach the problem 
from entirely different directions. There are a few cases of serious 
discordance, and some systematic divergence, but nothing so large 
as might have been expected. The agreement with (29) has been 
commented upon by Eddington himself. 

The geometrical mean mass of the twenty-four Cepheids in 
Table XVI is 6.2, which happens to agree almost exactly with that 
derived from the principle of equipartition. This can scarcely 
be more than accidental, however, for the velocities of only fifteen 
of these stars are available. 


t Monthly Nolices, 79, 5, 1918. 








216 FREDERICK H. SEARES 


I5. AVERAGE HEAT CONTENT AND CEPHEID VARIATION 


On the basis of the theory of radiative equilibrium Eddington’ 
finds for the central temperature of a star 


Tx pi'MiBm (31) 


where p is the mean density, 8 the ratio of radiation pressure to 
gas pressure, and m the atomic weight which is here assumed to 
be constant. The formula may also be written 
T° (32) 
R being the radius. R=D/2 is given by (18) or an equivalent 
expression; 6 is a function of %, and % in turn appears in his 
formulae as a function of the luminosity alone, which can be obtained 
from (25). 7, can therefore be calculated for any star. For 
thirteen Cepheids with periods ranging from 2 to 17 days Edding- 
ton? finds that 7,/8 is sensibly constant. There are accidental 
deviations, but no progressive change. Moreover, for two theo- 
retical stars with periods of 4.5 and 30.8 days the values of 7./8 
are practically identical. This suggests as a possible condition for 
Cepheid variation 
“p= const., or pM? =const. (33) 
Since %?/R is proportional to the total heat-content of a star,’ 
the physical significance of (33) is that the average heat-content 
per unit mass is a constant, or, in another form, that the gravi- 
tational potential at the surface of the star is constant for the 
Cepheids. Shapley‘ has discussed this condition from the stand- 
point of the relation of period to absolute magnitude and finds 
agreement with the observed period-luminosity relation (25). 
We may proceed somewhat differently by combining (24) and (26) 
in accordance with (33). This gives the theoretical period- 
luminosity relation 
M=-—s5 log P+j+const. (34) 


* Monthly Notices, 77, 601, 1917. 2 Ibid., 79, 5, 1918; 181, 1919. 
3 Jeans, Problems in Cosmogony, p. 191, 1919. 


4 Mt. Wilson Contr., No. 190, p. 7; Astrophysical Journal, 52, 79, 1920. 
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A comparison of the values of M from (34) for the twenty-eight 
stars in Table XVI with those from the observed relation, (25) 
above, gives results which are shown graphically in Figure 5. 
The agreement is not good, the divergence amounting to about a 
magnitude between log P=o.25 and 1.5, with even larger differ- 
ences for the cluster-type Cepheids, which are not shown in the 
diagram. 
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Fic. 5.—Absolute magnitudes of Cepheids based on the assumption of constant 
average heat content. The points are values for individual stars from equation (34). 
The circles, from equation (36), show the agreement with Shapley’s period-luminosity 
relation when the term involving the specific heats in Eddington’s expression for the 
density is included. 
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It will be noted that (33) is an empirical result which follows 
from an application to Cepheids of the period-luminosity relation 
(25) and the general theory of radiative equilibrium, but its 
derivation does not in any way involve the theory of Cepheid 
variation itself. If valid, it ought, however, in combination with 
the latter theory, to lead back to relation (25). But the compari- 
son, as made, does not give complete agreement. Since the 
theory of Cepheids has entered only through (24), the difficulty 
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presumably lies with that equation, which has been assumed to be 
applicable. If instead we use the corresponding relation derived by 
Eddington’ himself, the difficulty disappears. For (24) we must 
substitute 

log p= —2 log P—2 log (ya)!+ const. (35) 


where (ya)! is a function of 8 and of the ratio of the specific heats, 
and changes slowly with the mass. With this emendation, (34) 
becomes 

M=-—s5 log P+j—s5 log (ya)!+const. (36) 


The factor depending on the specific heats enters to the fifth power 
and thus becomes important. Its neglect in (34) causes the 
divergence shown by Figure 5s. 

The necessary corrections are easily calculated, however, with 
the aid of Table V of Eddington’s paper. The values of 1—8 
required for the interpolation of (ya)! are given in M.N., 77, 602, 
1917, under the heading “ Molecular Weight 2,” with the mass as 
argument. For J we may use the values in Table XVI, or more 
consistently with the present procedure, we may calculate % from 


log M=—o.2 (M—j)+ const. (37) 


which follows at once from (18) and the condition %«R given by 
(33). The latter method was adopted, the constant in (37) being 
determined so that the mean % from (37) agrees with that from 
the diagram in Figure 3. We are interested only in the variation 
of (ya), which is practically the same for all possible values of the 
ratio of specific heats. I'=1) was used, however, since Edding- 
ton’s results show that this agrees well with the data of observation. 
With the constant in (36) equal to 1.48, the values of —5 log (ya)! 
are as shown in the third column of Table XVIII. The sums of 
the remaining terms on the right of (36) are given in the second 
column under the heading /,. These are the co-ordinates of the 
dotted curve in Figure 5 and were read directly from a large-scale 
diagram. The calculated values of M from (36) (fourth column) 
when compared with those from (25) give the residuals in the last 


* Monthly Notices, 77, 15, 1918. 
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column of the table. For values of P between 3 and 30 days the 
agreement is very close. For P<2 days there is still a large 
divergence, but that is to be expected. For all but the cluster- 
type Cepheids, therefore, the theoretical relation (36) is practically 
the equivalent of the observed period-luminosity relation (25). 


TABLE XVIII 


THEORETICAL AND OBSERVED PERIOD-LUMINOSITY RELATION 











M 
Log P M, —5 log (ya)4 O—C 
(36) (25) 
0.25 —0.45 —0.45 —0.90 —1.08 —o.18 
0.50 —1.30 —0.29 —1.59 —1.58 +0.01 
0.75 —2.25 —o0.08 —2.33 —2.26 +0.07 
1.00 —3.34 +0.15 —3.19 —3.15 +0.04 
1.25 —4.43 +0.45 —$.90 —4.05 —0.07 
1.50 —5.54 +0.68 —4.86 —4.95 —0.09 




















Shapley’s discussion referred to above starts with the assump- 
tion of constant heat-content and leads to a theoretical relation 
which agrees well with (25), but does not involve Eddington’s 
theory of Cepheids because his comparison is based on equation 
(24) instead of (35). The close agreement found by Shapley is 
due partly to the fact that his final formula neglects 2.5 log (ya)? 
instead of 5 log (ya)’, which by itself would lead to a divergence 
from (25) equal to one-half that shown in Figure 5, and partly 
to his readjustment of the mean atomic weight. 

The matter is of interest from two or three points of view: 
First, because it shows the consistency of Eddington’s theory of 
Cepheid variation with the general theory of radiative equilibrium, 
in that it is necessary to include the term depending on the specific 
heats in order to work back to the observed period-luminosity 
relation which was one of the premises; and second, because it 
exhibits the empirical result (33) in another form and possibly 
gives weight to the suggestion that constant average heat-content 
may be a determining circumstance in Cepheid variation. 

The fact that the term which has been under discussion enters 
to so high a power subjects Eddington’s theory of Cepheid varia- 
tion to a severe test. Something of the sort is also true of the 
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empirical condition (33). Assume for a moment that the relation 
M«R is not exact, and suppose that all the variations which occur 
are thrown into the mass. The maximum systematic deviation 
in log %, by (18), (22), and (37), will then be equal to 0.4 of that 
shown by the differences in the last column of Table XVIII, 
which is well within the uncertainty of the comparison. 

This close agreement justifies the mention of other relations 
which follow at once from the assumption of constant average 
heat-content, combined with the following equations of radiative 
equilibrium :* 

Tex g(1—B) « detp!(1 —8) | 
T= Mfp'p f (38) 


1—Ba M85 


/ 


where 7, and 7, are the effective and central temperatures, respec- 
tively, g the acceleration due to gravity, and 1— 8 the ratio of 
radiation pressure to gravitational attraction. Assuming 4<R 
we find then 
M?p=c,; Mg=c.,; Rg«c;; px g*) fail 
T,.<B; T,.<T,Mix MIBx Tp} : 


where ¢,, ¢., and c; are constants. Further, with the aid of (35) 


dea (ya)'P; T.«T,(ya)'Pt« B(ya)*P?. (40) 


The relation T,< ‘8 is of special interest. Since 8 diminishes 
faster than J%<? increases, 7, must decrease with increasing J%. 
The change for the range in % shown above corresponds to an 
increase of about one spectral interval, which agrees well with the 
observed range in spectral type. At the same time M decreases 
from —1.2 to —5.1. 

The formulae of radiative equilibrium admit of no dispersion 
in the case of individual stars. This is also true of the relations 
(39) and (40) for the Cepheids. In the case of the masses and 
densities, we have already found reasons for believing that this is 
not strictly in accordance with the facts. Although the dispersion 
is small, it is in part undoubtedly real, for at present the theory 
necessarily neglects the modifications of gravitational acceleration 


* Eddington, Monthly Notices, 79, 180, 1919. 
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produced by the rotation of the stars on their axes, which must 
influence the relation of radiation to other stellar characteristics. 


16. RELATION TO RUSSELL’S FREQUENCY DIAGRAM—REMARKS ON 
THE DISTRIBUTION OF MASS AND DENSITY 

The justification for the discussion of Eddington’s theory of 
Cepheids in this place is the light thrown upon the distribution of 
mass and density. Equations (35) and (37), for example, afford 
a basis for a further revision of the mass- and density-curves; 
but an additional approximation seems scarcely justified at present. 
Table XVII shows already a surprisingly good agreement between 
his results and those in Figure 3. Further modification would 
affect mainly the stars of masses 1o or larger, and these are not 
numerous. Moreover, Eddington himself has repeatedly empha- 
sized the fact that his theory is but a first approximation, in that 
certain simplifying assumptions had to be introduced in order to 
make a beginning. Figure 3 is therefore allowed to stand as it is 
for the present, although the data in Table XIX are added to show, 


TABLE XIX 


Mass, ABSOLUTE MAGNITUDE, AND SPECTRUM, EDDINGTON’s THEORY OF CEPHEIDS 














| 
dM | M | Sp. M. | M Sp 
Oe ues aa —0.9 Fr.o | re eee —3.4 F8.0 
3 —1I.4 F4.0 Ws a catdres —3.6 F8.4 
4 —1.9 F5.2 Se —4.0 F9.0 
6 —2.5 F6.5 | Sere —4.6 F9.3 
eee | —3.0 F7.5 | a —5.1 F9.6 











The values of the masses are adjusted to give a mean value equal to that derived from Figure 3. 


in a general way, what would result, for the Cepheids at least, from 
rigorous application of the theory. The values of M and the 
spectrum may be regarded as co-ordinates for the location of corre- 
sponding values of % in a diagram similar to Figure 3. 

It is, of course, not to be supposed that the stars are uniformly 
distributed along the mass-curves of Figure 3. In general they 
cluster about certain lines of maximum frequency which, with one 


exception, bear little relation to the lines of equal mass. The most 
conspicuous of the frequency lines is that defined by the dwarfs, 
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to 
No 
No 


which continues through the A’s and joins smoothly with the line 
of maximum frequency of the B stars. Along this line we have 
found a regular decrease in the average mass with increasing 
spectral type and absolute magnitude. 

The number of A stars of very high luminosity seems to be 
relatively small.*_ If this is really the case, there must be a region 
of low frequency in the diagram between the most luminous B’s 
and the very bright giants of the later types. Beyond Fo the 
giants fall into two groups, the Cepheids and pseudo-Cepheids, 
and a large group of G, K, and M stars whose absolute magnitudes 
lie mainly between M=—o.5 and +2.5. The Cepheids and 
pseudo-Cepheids attain very high luminosities and show a rather 
definite correlation line running approximately from M = —o.4, F4, 
to M = —4, G8, which in general cuts the mass lines obliquely. For 
the Cepheids alone, its course is determined by the intersections 
of the mass lines of Figure 3 with those defined by equation (37), 
which, empirically at least, is a condition of Cepheid variability. 

The frequency line for the second and more numerous group of 
giants seems, on the other hand, to be definitely related to the 
distribution of mass. It has a clearly marked minimum at Ko 
and coincides approximately? with the mass line = 2 as shown in 
Figure 3. The chief difference is a steeper descent from Go to Ko 
and an even more pronounced minimum than is shown by the mass 
line. The coincidence is probably significant, for Eddington: 
finds on theoretical grounds that, among the giants, masses of this 
order of magnitude may be expected to have the highest frequency. 

t The frequency diagram of Russell for the helium stars of Kapteyn and the stars 
in the list of spectroscopic parallaxes is illustrated in Annual Report, Mt. Wilson 
Observatory, Year Book Carnegie Institution of Washington, 1921. 


2 The position of this correlation line and of the line of maximum frequency for 
the G, K, and M giants is shown in Fig. 1. This figure also shows the course of the 
equal-mass line 4=2. The correlation of spectral type with absolute magnitude for 
the most luminous stars of our own immediate system is essentially chat found by 
Shapley for the brightest stars in globular clusters; and it now appears from results 
not yet published that in some of the galactic clouds a similar correlation with apparent 
magnitude is to be found. It should be possible in such cases to estimate roughly the 
distances of the clouds. Provisional results indicate values of the order of 20,000 to 
50,000 light years. See Annual Report, loc. cit. 


3 Report of the British Association for the Advancement of Science, 1920, p. 43. 
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These circumstances seem to afford justification for the peculiar 
course of the mass lines for the late-type giants, which otherwise 
must depend solely upen the applicability of the principle of equi- 
partition. 

The minimum itself, or what remains of it after the reduction 
of the visual absolute magnitudes to bolometric values in the manner 
indicated above, must not be taken too seriously, however, for 
the residual irregularity in the mass lines is little in excess of the 
small systematic uncertainties affecting the absolute magnitudes of 
the different spectral types.' 

In this connection the result of page 183 for the masses of fourteen 
giant visual binaries may be recalled. Divided by 1.75 to reduce 
to the central component, the mean mass becomes 1.70.9. For 
the same stars we find by interpolation from Figure 3, with spectral 
type and absolute magnitude as arguments, the mean value 
M=2.7+0.2 The agreement is not good, but falls within the 
limits of uncertainty incident upon the small number of stars 
available for the comparison. 

The decrease in mass along the dwarf branch from %=10 at Bo 
to 4¢=0.6 at Ma raises a question of much interest, especially when 
it is recalled that the dispersion in mass is certainly small. Prac- 
tically, large masses are not to be found in our lists of dwarf stars. 
And yet the large-mass early-type stars are numerous. Why, 
therefore, do we not find more large masses among the dwarfs ? 

Several possibilities must be considered. Thus we might assume 
that the stars now on the dwarf branch (including the B’s, as 
usual) began their development at much the same time, and that 
those of largest mass have only partially run their course and are 
now in the B stage. The assumption would seem less arbitrary 
could we suppose that all the stars began their evolutionary careers 
at the same time; but this apparently is not the case, for numerous 
objects with only moderate masses of 2 or 3 are still in an early 
stage of development. The most serious objection to this explana- 
tion, however, is that it takes no account of a peculiar kind of 
selection. 

«Stromberg, Mt. Wilson Contr., No. 220; Astrophysical Journal, 55, 11, 1922. 


* Russell has also discussed this correlation., Observatory, 37, 173, 1914. 
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Let us therefore make the contrary assumption and suppose 
that stellar development has continued so long that all possible 
masses really occur among the dwarfs of late type. Among these, 
the large masses, say of the order of 10, will be of very low frequency 
as compared with the modal value for all types together. Kapteyn 
and van Rhijn place the maximum of the frequency-curve of 
absolute magnitudes at M@=+7.7. The modal value of the mass 
therefore should be about 0.7, and since the dispersion is small, 
the large masses must be very infrequent. 

Consider now only the stars in our lists of data. The dwarfs 
among these will all be included in a restricted region of space 
near the sun. The small masses will greatly predominate and the 
mean mass will be low. The giants, on the other hand, will be 
scattered through a large volume of space, and their average mass 
will be large, because only stars of large mass can attain high 
luminosity. The B’s among them will form the collection at the 
upper end of the dwarf branch whose mean mass we have found to 
be 10. The total number of these in our catalogues will be con- 
siderable, because they represent those present in a large volume, 
whereas the large masses among the dwarfs will be absolutely as well 
as relatively infrequent, because we can see only those which are 
near us. 

The stars of intermediate luminosity fill in the gap, and we 
must therefore expect to find a correlation of mass with absolute 
magnitude and spectral type very similar to that actually observed. 
The average mass decreases with advancing type; at the same 
time the number of stars of large mass in our catalogues decreases 
rapidly and the smaller masses appear in succession and in increas- 
ing numbers because they can attain the temperatures correspond- 
ing to the later spectral types. 

One other possibility should also be borne in mind, namely, a 
decrease in mass through loss of energy by radiation. From the 
standpoint of relativity" any change dE in the total energy of a 
system implies a corresponding change of dE/c? in the inertial 
mass, where c is the velocity of light. A similar conclusion follows 


t Einstein, Annalen der Physik, 18, 639, 1905¢ 
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also on the basis of the Newtonian mechanics" from the phenomenon 
of radiation pressure, for, unless we attribute mass to radiation, 
there are certain cases, at least, in which neither the energy nor 
the motion of the center of gravity will be conserved. Since the 
principles of mechanics are inductions, we cannot be sure that they 
apply to physical phenomena with the degree of precision implied 
by assuming that mass actually changes with changes in the 
energy. But considerations of consistency are sufficient to give 
the matter importance, and it derives added interest from the fact 
that a decrease in mass has been invoked as a means of accounting 
for the enormous unknown supply of energy which is now generally 
recognized as an essential feature of stellar phenomena. The sug- 
gestion that the packing of electrons and the nuclei of hydrogen 
atoms to form nuclei of heavier elements, with an accompanying 
loss in mass which reappears in the form of great quantities of 
radiant energy, attaches itself naturally to the relation between 
mass and energy that harmonizes with fundamental mechanical 
principles, whether of Einstein or of Newton. 

Of these three factors which may affect the correlation of mass 
and spectral type, selection certainly enters in the manner indicated, 
but its quantitative evaluation is at present impossible. This 
leaves the operation of the other two factors entirely speculative. 
If, however, there is any appreciable change in the mass of the 
stars during their development, they will not follow the mass lines 
of Figure 3, but cut across them obliquely; and this leads to a 
peculiarity of these lines that immediately attracts attention. 

The equal-mass lines do not meet the dwarf branch tangentially, 
as one might expect on the basis of a simple gravitational con- 
traction with its accompanying changes of temperature, but inter- 
sect the frequency line at large angles. A comparison of Figures 2 
and 3 shows that in this particular there is a large element of uncer- 
tainty; and yet the evidence seems to exclude the possibility of 
tangential intersections, unless we are prepared to abandon the 
principle of equipartition at the point where it seems to be best 
justified, namely, in the vicinity of the dwarf branch. 

‘ Tbid., 20, 627, 1906. This result seems also to have been pointed out by 
others—-Planck, for one, I believe—although I have been unable to locate the 


source, 
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The space-velocities of stars of any given type, say F5, show 
a definite correlation with absolute magnitude, the mean velocity 
increasing continuously as the brightness decreases. If equiparti- 
tion holds for the stars of a given type which are near the frequency 
line of the dwarfs, as it does in the mean for the stars of different 
types scattered along the frequency line, there must be a continuous 
decrease in mass downward along the vertical lines of Figure 3. 
The data of Table X give a numerical expression to this relation- 
ship of mass to luminosity, which, quantitatively, should be most 
exact for values of MW near the frequency line of the dwarfs. Now 
on this line, for spectral type F5, as an illustration, we find a 
mean mass of 1.5. Above this point in the diagram, that is, for 
more luminous stars, the masses are larger; below they are smaller. 
If the decrease in mass downward is really continuous, it follows 
at once, since the mass also decreases continuously along the dwarf 
frequency line, that curves connecting equal masses cannot meet 
the frequency line tangentially, but must intersect it at a consider- 
able angle. In fact Figure 2 shows that the strict application of 
the principle of equipartition leads to angles in excess of 90°, 
counting as o the tangential junction required by gravitational 
contraction. The evolutionary path of any given star, however, 
when once it has joined the dwarf branch, must lie in the general 
direction of the frequency line. 

The explanation of the contradiction lies in the peculiar selection 
which determines the decrease in mass along the dwarf branch. 
Consider for example two dwarf stars of Type Fs5, one a large 
mass which in its development has followed the dwarf branch 
downward from some earlier type, the other a star of mass so small 
that it has been unable to attain the high temperatures of the 
very early types and joins the dwarf branch at Type Fs5. The 
star of large mass will be the brighter of the two. Hence the path 
of the small star must intersect the line of maximum frequency at 
some point a little in advance of F5 in order that it may swing 
into position below the larger star. The path of the larger star, 
on the other hand, will lie above the frequency line. Similar 
reasoning applies to stars of intermediate mass and it is evident 
that the entire group of F5 dwarfs will show a correlation of mass 
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with absolute magnitude similar to that found from the principle 
of equipartition. In their further evolution these stars will move 
as a group along paths approximately parallel to the line of maxi- 
mum frequency. But these paths do not and cannot coincide with 
the equal mass lines of Figure 3, for when the group in question is 
at Fs, there is at Go another group, whose luminosities show a similar 
correlation with mass, but whose masses, because of selection, are 
systematically smaller than those at F5. If the points in the dia- 
gram corresponding to equal masses in these two simultaneously 
observed groups be connected, the lines of equal mass will cut the 
frequency line at high angles as shown in Figure 3. Hence there 
is no conflict between the mass-distribution shown by this diagram 
and the gravitational theory of stellar development. The only 
point is that in the vicinity of the dwarf branch the stars do not 
follow the equal mass lines. 


17. IONIZATION AND MEAN ATOMIC WEIGHT 


The final equations of Eddington’s theory’ are 


b= (1—8) (41) 
1—B=0.0026 M?Bim4 (42) 


where c and G are constants. These formulae refer to the giant 
stage and assume that the mean atomic weight is independent 
of the temperature. Variations in m during the star’s development 
would affect the relation between luminosity and spectral type. 
The absolute value of m fixes the relation between mass and 
radiation-pressure 8, and hence determines the values of the mass 
which are most likely to occur. 

It is unnecessary to repeat here Eddington’s arguments in favor 
of a small value of m. . In general these seem to be only strength- 
ened by the results of recent investigations of ionization phenomena. 
Eggert,? for example, has considered the behavior of the iron atom 
when subjected to pressures and temperatures comparable with 
those at the center of a typical giant star and finds that 16 of the 


t Astrophysical Journal, 48, 208, 210, 1918. 


2 Physikalische Zeitschrift, 20, 570, 1919. 
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26 electrons would be stripped from the atom, thus reducing the 
mean “‘atomic’’ weight to 3.3. 

It is easy to extend this result to other parts of the star. Using 
the reaction isobar of Nernst and basing the calculation of the 
energy of dissociation on the dynamical relations of the Bohr atom, 
Eggert finds as the condition for the removal of 8 electrons from 
99 per cent of the atoms’ 

2.68+0.4 log p= — 2A 10'+log T; (43) 
in which the pressure P is measured in atmospheres. For the 
separation of 16 electrons the corresponding condition is approxi- 
mately the same as (43) except for the appearance of an additional 
factor of ro in the first term on the right. 

Equation (43) and its alternate may be used to calculate the 
values of 7; corresponding to the pressures at different points 
along the radius. The comparison of the ionization temperatures 
thus found with the corresponding stellar temperatures should 
then give some indication as to the variation in ionization through- 
out the star. 

The total pressure P is given by the ordinary gas equation, 


divided by £, the factor arising from radiation-pressure. 


pole R pT (44) 
B Bm 

T and p may be found from Emden’s’ formulae, or we may use 

the values given by Eddington’ for his typical star (J¢=1.5, 

Pm=0.002), taking care to modify those for 7, which contain the 

factor m, to make them conform with the degree of ionization 


expressed by (43) and its alternate. Since 7 contains both m 


* Professor Paul Epstein of the California Institute of Technology calls my 
attention to an error in Eggert’s expression for K (op. cit., p. 573 \ factor 8° has 
been omitted. It has little effect, however, on equation (43), the absolute term being 
increased from 2.32 to 2.68. 

2 Gaskugeln, p. 97, 1907. Eddington has called attention to the misprint in the 
formula for 63. The last factor in the denominator should be squared. 

3 Astrophysical Journal, 48, 213, 1918. Unless I am in error, Eddington’s values 
for the temperature have been computed with m= 2, whereas the text seems to imply 
that 2.8 has been used. If so, the tabular values should be increased by 40 per cent. 
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and £ as factors, the total pressure P depends upon neither the 
value of the radiation-pressure nor the degree of ionization. T7, 
however, depends upon both, and, moreover, 8 also depends upon 
m (see equation [42]). 

For the removal of 8 and 16 electrons, respectively, the values of 
m become 6.2 and 3.3. Corresponding to these, for M= 1.5, 
8=0.496 and 0.766. With these data we find the values of P 
and 7 given in the second, third, and fifth columns of Table XX 


TABLE XX 


IONIZATION AND STELLAR TEMPERATURES, IRON IN EDDINGTON’S TYPICAL STAR 








j 
| PERCENT- 








| 
| 
| —8 ELECTRONS m=6.2 —16 ELECTRONS m=3.3 
| toe 4 AGE OF 
r ee Cf eens | MAss 
rt | | | OursDE 
| r qT; 7 | T; r 
Divscsciteed 7.40 8.7 xX 108 8.0X 105 7.1X10° 2.6X 108 100 
Bid ode aes 7.43 7.4 6.97 | 6.4 2.4 88 
er io 5.0 4.9 4.1 | 2.0 48 
Bs cates 5.63 £.4 3.4 $3 | 1.6 18 
Bin cans 4.69 1.8 2.5 r.5 rs 5 
er ee 3.50 0.96 1.8 ©.79 | E53 0.7 
6 1.98 0.39 Es 0.32 | 0.8 0.005 
| 














Radius of star=6.9, %=1.5, Py, = 0-002. 


for various points along the radius. The unit for 7 is such that the 
distance from the center to the surface is 6.9. Substituting 
the values of P into Eggert’s formula, we find the ionization 
temperatures given in the fourth and sixth columns of the table. 
These temperatures, required for the removal of 8 and 16 electrons, 
respectively, may be compared with the corresponding stellar 
temperatures. The latter, which change little with the assumed 
change in m, are approximately ten times the temperature required 
to separate 8 electrons from the nucleus, and of the same order as 
the ionization temperatures corresponding to the removal of 16 
electrons. The decrease in pressure toward the surface compen- 
sates for the fall in temperature to such an extent that the ioniza- 
tion should be nearly constant throughout the star. It is only 
within a thin shell near the surface, which contains an insignificant 
fraction of the mass (see last column, Table XX), that m can rise 
much above the value 3.3. 
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Eggert’s formula is an extrapolation and involves assumptions 
which cannot at present be directly controlled; but its general 
indications are in close agreement with Eddington’s conclusions. 
Even admitting the various underlying assumptions, the result 
applies only to a star consisting wholly of iron; but one can scarcely 
doubt that the mean values of m for all the elements actually 
present must be small and nearly constant for all points within the 
star except those close to the surface. 

It seems worth while also to consider the relation of ionization 
to stellar temperature in giant stars of constant mass but differing 
spectral types, in order to gain some idea of the permanence of 
this relation as a star contracts under the action of gravitation. 

Consider stars of Eddington’s critical mass 4% =2, which, as we 
have seen, are those occurring most frequently among the late- 
type giants. It will be sufficient to calculate the central tem- 
peratures and pressures, and in view of the preceding results it is 
clear that we may use m=3.3. The corresponding value of @ is 
0.70. Using Emden’s formula as before (remembering that for 
the type of equilibrium considered p,=54.25 pm), we find the 
results for 7, and P, given in Table XXI. The values of M and 


TABLE XXI 
CENTRAL TEMPERATURES AND PRESSURES, AND IONIZATION TEMPERATURES FOR IRON 
(S=2, m=3.3) 











| | 
Sp. | M Log Pm | Log ie peak T, T; 

isin = Ee SS ES aaa eee are 
nes: 0.0 | —5.1 | 6.10 4.3 I.3xX10° I.2X10° 
EE ee re +0.7 —4.1 6.43 5.6 2.7 | r.6 
Eee +1.1 —3.1 6.77 7.0 5-9 | 23 
ene +o.8 | —2.7 6.90 7.3 7.9 | 2.6 
Be o.5 tienes +o.7 | —2.2 + .O7 8.2 EE. z.2 

, ee - df +2.0] —1.0 7.47 | 9.8 | 20.5 7.0 





log p» in the second and third columns were read from the diagram 
of Figure 3. The ionization temperatures for the removal of 


16 electrons are in the last column of the table. 

These data represent the sequence of changes from an early 
stage of development, on the basis of Russell’s theory, through to 
the dwarf stage, where the density is so large that Boyle’s law is no 
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longer accurately obeyed. The increases in 7,, p and P, are 
approximately 20, 10,000, and 300,000 fold, respectively. In 
spite of these extraordinary changes, the ionization temperatures 
remain of the same general order as those at the center of the star 
throughout its development. Initially, the two temperatures are 
sensibly equal, and the gain in 7, over 7; is so slight that any 
further decrease in m must be small. The complete removal of 
the third sheath of electrons, for example, would give m= 2.2. 

Even for stars of the lowest density, such as Antares or Betel- 
geuse, conditions are very similar to those illustrated by the first 
line of Table XXI. Low density means relatively low central 
temperature and pressure, but this factor is compensated by the 
large masses which characterize the very luminous giants, so 
that the relation of ionization to stellar temperatures is much the 
same as for the objects already examined. 

It may be remarked incidentally that the values of T, in 
Table XXI throw some light upon the behavior of the effective 
radiating layers in stars of different types. The increase in 7, 
is approximately ten times that in the effective temperature. 
But for the type of equilibrium in question the temperature gradient 
is independent of the density, that is, the temperature increase for 
homologous points along the radius should be constant. The 
explanation must be that the thickness of the radiating layer in an 
Ma giant, because of the enormous difference in density, is very 
much greater than that for an F5 dwarf. As development pro- 
gresses, radiation from the deeper and hotter strata is absorbed 
as a result of increasing density, and, in consequence, the increase 
in effective temperature is less rapid than would otherwise be the 
case. 

On the basis of the preceding evidence we must expect a high 
degree of ionization at the earliest stages of stellar development 
with which we are familiar, and comparatively little decrease in 
the value of m as the development proceeds. The Bohr atom 
affords the only foundation we have for estimates such as have 
been made, and it is not certain that the results are of the right 
order of magnitude. One consequence of their acceptance, how- 


ever, should not be overlooked, and that is the total amount of 
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energy required for ionization. For the removal of 8 electrons 
this equals 1.8107 calories per gram atom," corresponding to 
an ionizing potential of about 800 volts. The mass of Edding- 
ton’s typical star is approximately 2.9 X10*% grams. The removal 
of 8 electrons from all the atoms in this mass would therefore 
require 
2 33 2 7 hi ent 
va! OX 108X I ma 107X4.2XI10 24X10" cens, (45) 
while for the removal of the second sheath of 8 electrons a quantity 
of the order of ten times this amount, or 
U,=4X 10" ergs (46) 
would be necessary. 
Now the total energy generated by gravitational contraction 
of the star from a state of infinite diffusion is? 


3(y—1) GM? 
Q=' (47) 

5sy—-6 R 
where y is the ratio of specific heats, G the gravitational constant, 
6.66 X 107° cm3/grams sec.?, and R the radius, 7x10" cm; y lies 
between the limits 4/3 and 5/3, whence the fractional coefficient on 
the right has the limiting values 3/2 and 6/7. At a maximum, 
therefore, which also corresponds to Eddington’s model, 
10-7X9X 10% | 


he —_' I.3X10* ergs. (48) 


wee 


Of this total the etherial energy alone, namely, 
(1—B)Q=0. 234 Q=3.0X 10% ergs, (40) 


t Eggert, loc. cit. 

2 Emden, Gaskugeln, p. 125, 1907. 

3 Eddington, Monthly Notices, 79, 23, 1918. Anderson, with whom I have had 
much stimulating discussion during the preparation of this paper, calls my attention 
to the fact that ionization cannot be produced by the collision of electrons with atoms 
without violating the Nernst formula of equilibrium, which is well established on 
thermodynamical grounds. This equation depends upon the pressure, whereas the 
hypothesis of ionization by collisions alone leads to a condition for equilibrium which 
is independent of pressure. If the effective energy is assumed to be partly radiant 
and partly kinetic, the pressure would enter, but not in the manner required by 
Nernst’s equation. Since the internal kinetic energy of the atoms is not available, the 
radiant energy alone can be operative (as in the case of the photoelectric effect, for 
example). See also Milne, Observatory, 44, 269, 1921. 
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seems to be available for ionization. This is of the same order as 
U,, the energy of dissociation given in (46). Admitting, therefore, 
the validity of Bohr’s theory and the applicability of the Nernst 
formula. of equilibrium, the insufficiency of the gravitational 
theory presents itself in a new form. As an alternative, the 
question is perhaps raised as to whether the giant stars, even in 
their very earliest stages of development, are ever in a state other 
than one of high dissociation. 

The calculation of the stellar temperatures in Tables XX and 
XXI neglects altogether the depletion of the available supply of 
energy by ionization. Were gravitation the only source, we could 
expect neither the temperatures nor the degree of ionization 
indicated by the discussion; but, instead, an equilibrium state 
which would be established for something less than the extremes 
of temperature and ionization shown by the tables. 


18. COMPARISONS WITH EDDINGTON’S FORMULAE 


The assumption concerning m which underlies equations (41) 
and (42) has been discussed in the preceding section. These 
formulae also assume that the outflow of energy across any surface 
of unit area is proportional to the gravitational acceleration at 
that point (H<g), and that the mass coefficient of absorption, &, is 
a constant. The energy assumption can also be stated in the 
form that the outflow throughout the star per unit of mass is 
constant. The source of the energy is not specified; it may be 
derived from gravitational contraction or any other source capable 
of providing an adequate supply. Assuming further that the 
transfer of energy is by radiation and not by convection, the condi- 
tions for equilibrium are found to be the same as those for adiabatic 
equilibrium with a ratio of specific heats equal to 4/3. It then 
follows through (41) and (42) that the bolometric magnitude of a 
star remains constant so long as it behaves as a perfect gas." 

A giant star in its development, however, presumably follows 
closely one of the equal-mass lines in the upper part of Figure 1, for 


t As Jeans points out, Monthly Notices, 79, 319, 1919, this general result follows 
at once, by simple considerations, from the assumption H/g=const. The nature of 
the dependence of radiation on mass and other characteristics is, however, another 
matter, less easily disposed of. 
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the selection which enters along the dwarf branch has comparatively 
little influence on the average mass of the giants which come under 
observation. Observations indicate, therefore, that the absolute 
magnitude is not independent of the stage of development. The 
change is even larger than shown in Figure 3, for, as already stated, 
the visual absolute magnitudes must be corrected by the quantities 
given on page 194 to make them comparable with the bolometric 
magnitudes calculated by (41). 

The question now arises as to the particulars in which the under- 
lying assumptions must be modified in order to effect a closer 
agreement with observations. Such consideration as can be given 
here may be prefaced by the remark that Jeans,’ starting from 
the assumptions of radiative conduction, gravitational contraction, 
and constant & and m, finds that the radiation of a star will remain 
constant as long as Boyle’s law is obeyed. This is also the case 
if the star has sources of energy other than gravitational contraction, 
provided the rate of generation per unit mass is independent of the 
time. This is Eddington’s result without his assumption that 
H«g; but the functional dependence of radiation upon the mass 
cannot be specified. Jeans also finds that H/g will be constant 
throughout the mass of gas only when the mass is so large that 
radiation-pressure predominates and 8 may be put equal to zero. 
For the ordinary run of stellar masses, however, H/g will not 
be constant within the star, but for a given point will be independent 
of the time, provided the only source of energy is gravitational. 

Eddington’ himself points out that H/g cannot be rigorously 
constant within the star, for, assuming gravitational contraction 
through a series of homologous states and a temperature- and 
pressure-distribution corresponding to adiabatic equilibrium with 
y= 4/3, the outflow of energy per unit mass would increase about 
70 per cent from the surface to the center. 

In view of a probable variation of H/g of the order indicated, 
the original assumption gives close agreement with observations. 
The matter can be tested by comparing the calculated values of M 
corresponding to different values of J% with that shown by Figure 3 

t Loc. cit. 

2 Monthly Notices, 77, 599, 1917. 
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for stars of the same spectral type. By proceeding in this manner 
we avoid disturbances arising from variations in k and m, since 
for stars of the same type these quantities should be less subject 
to change. | 

The results are shown in Table XXII. This table gives the 
deviations of the observed values of M (Fig. 3) from those calculated 


TABLE XXII 


OBSERVED AND CALCULATED (EDDINGTON) VARIATION OF ABSOLUTE MAGNITUDE 
WITH Mass 
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by Eddington, which are in the second column of the table. Zero 
point corrections have been applied to the latter to make the mean 
O—C for each type equal to zero. There is a large progressive 
change for the F stars, but the agreement improves with advancing 
type and, with the exception of one difference, is practically perfect 
for the K and M stars. From this comparison it is clear that the 
error in the assumption H/g=const. is not very serious when 
we consider stars of a specified type. 

The lack of parallelism of the equal-mass lines with the axis 
M =o is difficult to account for because the assumption of constancy 
for both & and m enters here, as well as that concerning the flow 
of energy. In general the luminosity for any given mass decreases 
with increasing temperature. For mass 4, as an illustration, the 
bolometric absolute magnitude increases from —3.1 at Ma to +0.3 
at Fo. This is in the direction corresponding to decreasing values 
of m, and one might expect much of the change to be accounted 
for by increasing ionization. The evidence of section 17, however, 
is against any considerable variation in m. Even in the earliest 
visible stages of development we should not anticipate a value 
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much in excess of 3.3, and unless there is a large amount of nuclear 
disintegration, which seems improbable, the value of m cannot 
fall below 2. It is easily shown that a decrease from 4} 2 to 2 
would cause an increase in M of about one magnitude, which is only 
a third of the amount required. The admission of much larger 
initial values of m is insufficient to bridge the gap, for (42) shows 
that 1—8, and hence the luminosity, are not much affected as m 
is increased. 

The behavior of & is obscure. Although the conditions within 
the star which determine the degree of absorption must be analogous 
to those governing the absorption of X-rays, it should be borne 
in mind that the effective value of k is that corresponding to the 
outermost layers of the star where the temperature is very much 
lower than in the interior. Even near the surface the opacity is 
large, but we have no data as to its variations that are certainly 
applicable. 

This leaves the larger part of the increase in the absolute 
magnitude unaccounted for, and it is not clear as to how the dis- 
crepancy is to be explained. If the contraction is homologous 
and the energy gravitational, or gravitationai plus a source inde- 
pendent of the time, we should not expect variations in outflow 
with advancing development other than those caused by changes in 
k and in m. 

As an explanation, changes in m that would be sufficient seem 
to be excluded; those of & are in doubt. There remain, however, 
the assumptions underlying the outflow of energy. The source 
of energy almost certainly is not wholly gravitational, and if not, 
with almost equal certainty, is not independent of the time. The 
outstanding questions therefore lead directly into the realm of 
atomic physics, and it seems more than ever likely that the problem 
of the stars will find its solution in the solution of the problem of 
the atom. 

In conclusion, attention may be directed to the relation of mass 
to luminosity along the dwarf branch. Here Boyle’s law is in 
general no longer obeyed. It will be noted, however, that the 
density increases slowly until the late types are reached, which 
suggests that the departure from the conditions of a perfect gas 
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may affect the behavior of stars of different spectral types much 
alike, so that a formula giving the relation of luminosity to mass 
for the giants, with the exception of a constant, would also hold 
approximately for the dwarfs. 

The matter is easily tested by calculating % from (41) for the 
values of L corresponding to the absolute magnitudes given in the 
second column of Table IV. Choosing the constant factor so that 
the mean of the values thus calculated agrees with the observed 
mean mass, we have the results in Table XXIII. With the 


TABLE XXIII 


VALUES OF MAss ON DWARF BRANCH 
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exception of a spectral interval at either end of the series, where 
the deviations from the average density are large, the agreement 
is good. 

I must express my grateful acknowledgments to Miss Joyner 
and Miss Richmond, of the Computing Division, for their efficient 
assistance in the extensive tabulations and calculations involved 
in the preceding discussion. 

Mount WILSON OBSERVATORY 

December 1921 








ON THE CALCULATION OF MASSES FROM 
SPECTROSCOPIC PARALLAXES' 
By HENRY NORRIS RUSSELL? 


ABSTRACT 


Calculation of masses of stars from spectroscopic parallaxes.—It appears probable 
thi at the line intensities upon which spectroscopic parallaxes are based are functions of 
the temperature and density of the star’s atmosphere. If this were exactly true, all stars 
of the same surface brightness and dens ity would have the same _ troscopi ic abaslete 
magnitude, and the masses, computed from the spectroscopic parallaxes, would come out 


] 


the same for all the stars of such a group whateve r the ‘dis spersion among their actual 
masses) and equal to the geometrical mean of the latter. To obtain a r i ible measure 
of the dispersion in mass among binary stars, parallaxes must be determined in some 
other way. 

Spectroscopic and dynamical parallaxes.—It follows that the spectroscopic paral- 
laxes and the dynamical parallaxes (derived on the assumption that the mass of a binary 
system is equal to the mean value for stars of its spectral type and absolute mag — ») 
are systematically equivalent to one another, and really rest on the same physical relation- 
ships and assumptions. 

Dis persion of mass among visual binaries (dwarf stars).—By an indirect n ethod, 
depending on Strémberg’s « omp arisons of spectroscopic and trigonometric parallaxes, 


Seares’s conclusion that the dispersion is small is confirmed. The probabk error of 
lot 


dispersion of log J% appears to be less than +o0.2, but cannot be exactly determined. 

The very small values obtained by Seares in the preceding paper? 
for the dispersion in mass among dwarf stars suggest the following 
explanation. 

It appears very probable, from physical considerations, that the 
spectral type of a star is determined by the temperature 7 of its 
outer atmosphere, while the characteristics associated with the 
absolute magnitude depend also, and mainly, upon the density 
of the atmosphere, p’. If 7, is the absolute magnitude, determined 
spectroscopically, we may then write 

M,=F(T, p’). (1) 

It is also probable that the surface brightness / depends almost 

entirely on the temperature, and that the density of the atmosphere 


* Contributions from the Mount Wilson Observatory, No. 227. 
2 Research Associate of the Mount Wilson Observatory. 
3 Mt. Wilson Contr., No. 226. 
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above the visible surface is a function of the mean density p 
(though not necessarily proportional to it). We have then 


M,=f(J, p). (2) 


It may be that M@, depends to some degree on other variables, 
but it is probable that their influence is slight. The form of the 
function / cannot be predicted by theory, but can be found empiri- 
cally when sufficient data regarding stellar masses and densities are 
available. 

Assuming the relation for the moment to be exact, and intro- 
ducing it into equation (22) of Seares’s paper, we have 


log M;=log p+0.6j—0.6/(J, p)+2.77. (3) 


Hence J, the mass corresponding to the spectroscopic parallax is 
itself a function of pand/j alone. Like (22) this is a general result, 
holding good no matter how the mass is calculated, so long as the 
formula involves the absolute magnitude and our assumptions 
about the latter are correct. 

Stated in words, this signifies: Jf the spectroscopic absolute 
magnitude of a star depends only on its surface brightness and density, 
the masses compuled from the spectroscopic parallaxes will be identical 
for all stars which have the same density and surface brightness, no 
matter how different their actual masses may be. 

In still other words, the spectroscopic parallax (on our assump- 
tions) will be correct only in the case of stars having a certain 
average value of the mass. For larger or smaller masses there will 
be a systematic error in the absolute magnitudes, making the 
former too faint and the latter too bright, which will alter the 
spectroscopic parallax to just such a degree as to conceal entirely 
the real differences in mass. 

If the spectroscopic calibration curves are adjusted so as to give 
a correct value for the mean absolute magnitude of the stars, the 
“spectroscopic mass’’ will be the geometrical mean of the individual 
masses, since, if 7 and p are constant, M varies proportionally 
to log M. 

To determine the dispersion in mass among visual binaries, it is 


therefore necessary to have recourse in some way to parallaxes 
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other than spectroscopic. One way of doing this is through 
Strémberg’s' discussion of the errors of the spectroscopic paral- 
laxes—which he has determined by comparison of individual 
spectroscopic and trigonometric parallaxes, allowing for the errors 
of the latter. Let the probable error thus found for log 7 be =r. 
Three sources contributory to this error may be distinguished: 
(1) the error just discussed, arising from differences between the 
masses of individual stars of the same surface brightness and 
density; (2) the error arising from other physical causes, on account 
of which the functional relation (2) is not exactly true; and (3) 
errors of observation. If we call the amounts of these, expressed 
as probable errors of log 7, respectively, 7, 7», and r,, we will have 


r=r,trtr 


If, on the other hand, we compute individual masses as Seares 
has done, the error r,, will be without influence on these masses, 
while r, and r, will affect log M% by 0.6 of their amounts, and there 
will be an additional source of error, due to the imperfections of the 
orbits and also to errors in the assumed ratio of the masses of the 
components, which we may call =r,. If then 7’ is the probable error 
of distribution, which Seares calls r(A log uw), we have 


rtrr+r=r =(0.36) 


0. 26(r.+r)+r=r’= (0.22) 


This would indicate that the actual dispersion in log M% (which 
corresponds to the probable error +o0.6 7r,,) is very small—being 
substantially equal to the spurious dispersion which arises from 
errors in the orbital elements of the binaries. 

The observed values of r and r’ are, however, subject to some 
uncertainty, and it is possible that some other concealed correlation 
may still be making the agreement appear better than it should. 

It is obvious, however, that a better way of determining the 
dispersion in mass among visual binaries of the same spectral type 
will be a direct comparison of “‘hypothetical’”’ or dynamical paral- 


* Mt. Wilson Contr., No. 199, p. 15; Astrophysical Journal, 53, 13, 1921. 
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laxes, computed with a mean value of the mass, with trigonometric 
parallaxes. In this case there seems to be no reason to fear a cor- 
relation of the sort discussed above. 

It may be noticed that the dynamical parallax of a binary, 
calculated on the assumption that the principal star has the mass 
dM, derived from spectroscopic parallaxes of stars of similar kind, 
and the spectroscopic parallax itself, must agree exactly in every 
case (barring the effects of errors of observation and the small 
uncertainty arising in the estimate of the ratio of the mass of the 
system to that of the brighter component), so that it may fairly 
be said that spectroscopic and dynamical parallaxes are systemati- 
cally equivalent to one another, so long as the mass used in com- 
puting the dynamical parallaxes is taken as a function of both 
absolute magnitude and spectral type. The two then rest essen- 
tially on the same physical assumptions—though this is far from 
obvious at first sight. If this conclusion is correct, spectroscopic 
parallaxes are unfitted by their very nature to give us information 
about the differences in mass between systems which have similar 


‘ 


spectra (with regard to the . 


as other characteristics). 


‘absolute magnitude lines” as well 

When a direct determination of the dispersion in mass among 
visual binaries has been made—for example, from the trigonometric 
parallaxes—it will be possible to obtain a valuable control upon 
the theory which atttributes the spectroscopic phenomena in 
question to differences in density. 

Mount WILSON OBSERVATORY 

February 16, 1922 








ON THE UPPER LIMIT OF DISTANCE TO WHICH THE 
ARRANGEMENT OF STARS IN SPACE CAN AT 
PRESENT BE DETERMINED WITH SOME CONFI- 
DENCE’ 

By J. C. KAPTEYN? anp P. J. VAN RHIJN 
ABSTRACT 


If both ¢(M) and Nm, the luminosity-curve and the total numbers of stars for 
each apparent magnitude, were completely known, the star-density could be found for 
any distance from the sun. The paper deals with the question: What can be reached 
now that our knowledge of these quantities is limited ? 

The conclusion is that the incompleteness of our knowledge of ¢() is of little 
importance, but that the contrary is true for Nm. Even if we assume the law of the 
densities, formula (2), to hold for all values of m, the constants cannot be determined 
with the accuracy needed for the derivation of the densities at considerable distances. 
The limits of distance for which the density can be determined within 25 or 40 per 
cent of its amount were finally found on the supposition that the values of Vm are 
known to m=14 and m= 20, respectively, and are givenin Table XII. With the limit 
m=17, visual, for which we may expect complete values of Nm in the near future, 
the densities should become pretty reliable for the whole of the domain within which 
the density exceeds 0.1 of that near the sun. 

The conclusions of the present paper must be considered 
only provisional, mainly for the two following reasons: (a) The 
extinction of light in space has been assumed to be inappreciable; 
(b) No use has been made of a recent substantial improvement.in 
the average parallaxes of the more distant stars. 

A more definitive solution now being made at the Laboratory 
of Groningen will duly take into account these two points. The 
progress of such a solution being necessarily very slow, the present 
provisional treatment was undertaken mainly for the purpose of 
finding the most promising lines for conducting the further treat- 
ment of the sidereal problem. 

The points which will be successively discussed are as follows: 

I. If for any determined part of the sky we know completely 
(a) the luminosity-curve y=¢(M) and (6) N,,=the total number 
ofistars of magnitude m—1/2 to m+1/2 per ten thousand square 
degrees, then we can find the star-density A(p) at any distance p. 

* Contributions from the Mount Wilson Observatory, No. 229. 


2 Research Associate of the Mount Wilson Observatory. 


242 

















ARRANGEMENT OF STARS IN SPACE 243 


The solution is contained in Mount Wilson Contribution, No. 188, 
formula (13). The difficulty of the problem lies exclusively in 
the required completeness. 

II. Half of the difficulty may be safely ignored, for it can be 
shown that at present and for a long, long while to come we may 
certainly adopt the form 


H 
€ 


Vr 


—H+(M— Me)? — eb +qM+rM?2 (1) 


o(M)= 


as representing the luminosity-curve completely, that is, for all 
values of M from —o to +o. 

We will here adopt the values of the constants p, g, r, given in 
Contribution No. 188, modified in accordance with the changes in 
definitions to be explained presently. 

III. There thus remains only the difficulty due to our incom- 
plete knowledge of N,,.. We will show what limitation is thereby 
introduced in the solution mentioned in I, and we will further 
try to show how the limits can be extended by future countings 
of the fainter magnitudes. 

In this inquiry, proceeding from the simple to the more difficult, 
we will consider three cases: 

a) We will assume that the formula 


. h ; 
N,.= e—/?(m— me)? — o¢9-+bm+cm* (2) 


Vr 
holds not only for the magnitudes for which we have the necessary 
data but also for all fainter magnitudes. 

b) We will drop the latter supposition, that is, in accordance 
with the results found thus far, we will adopt the error-curve 
(2) for the actually observed N,,, but will assume nothing about 
the fainter magnitudes which have as yet not been investigated. 

c) We will consider the most general case, namely, that the 
actually observed values of N,, themselves are not distributed in 
an error-curve. 

IV. Having thus found the limits within which a fairly reliable 
determination of the arrangement of stars in space can be made, 
we will introduce what we think must be considered a plausible 
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supposition, by which these limits, at least for the lower galactic 
latitudes, can be somewhat farther extended. 

Before entering on a discussion of these points it will be neces- 
sary to premise the following remarks: 

a) Cluster observations by Shapley have fairly well proved the 
absence of any selective extinction of light in space, at least for 
distances like those considered in the present discussion. It does 
not follow that there is no appreciable general extinction of light, 
and the matter will be carefully considered in our definitive solution. 
Still, in view of Shapley’s observations, we feel less hesitation in 
making the assumption of no appreciable extinction the basis of 
this provisional solution. 

b) Up to the present, the luminosity-curve has expressed 
the number of stars of each magnitude contained in a specified 
volume of space in the neighborhood of the sun. It would certainly 
have been much more natural and convenient had it been given the 
form of a true frequency-curve, that is, one showing the fraction 
of all existing stars belonging to each absolute magnitude. This, of 
course, is not rigorously possible as long as our observations do 
not embrace all the stars, even the very faintest in existence. But 
since the luminosity-curve has been derived (in Contribution 
No. 188) over a range of more than twenty magnitudes and since 
its coincidence with a normal error-curve throughout this entire 
interval, which extends far beyond the maximum, is truly astonish- 
ing, it would seem very desirable to introduce the assumption that 
the best-fitting error-curve represents the luminosity-curve over 
its whole extent. In this way it becomes possible to define the 
luminosity function as a true frequency-curve. 

Henceforth we will assume, therefore, that the luminosity- 
curve, over its whole extent, is perfectly representable by equation 
(1), which differs from formula (11) in Contribution No. 188, only 
in being divided by the factor A =0.0451, the number of stars per 
cubic parsec in the neighborhood of the sun. If later on it should 
appear that this assumption deviates sensibly from the truth, 
the final results will still be unaffected. The ordinates of the curve 
will have only to be multiplied by a certain factor, while the 
densities will be divided by that same factor. 
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c) The parsec is again taken as a unit of distance, and in 
accordance with this definition the absolute magnitude of a star 
will equal its apparent magnitude as seen from a distance of one 
parsec. 


TABLE Ia 
LUMINOSITY-CURVE 


(Equation [5]) 
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d) Finally, the star-density A(p) now means the total number 
of stars from the brightest to the faintest in one cubic parsec. 

Under these stipulations we have in (1), in accordance with 
Contribution No. 188, 


H =o. 2818, M,=2.6093 (3) 

p= —2.413' g=+0.4278 r= —0.079044 (4) 
log ¢(M) = —1.049+0.1858 M—0.03450 WM? (5) 
M=m-+ 5 log r=m—s log p. (6) 


The numerical values of (5), that is, of the luminosity-curve 
according to the present definitions, have been tabulated in Table Ia. 
For convenience we give in Table III the values of the con- 
stants of Contribution No. 188, p. 13, changed to agree with our 


TABLE Ib 


log A (p) (FORMULAE OF TABLE IV) 


| Galactic Latitude 








| 
log p | ——_—_————_ —_— —————___—__—. - — - 
°° 30 60 ) 4 ®) 
1.0. 5.05—10 8.65—10 8.65—10 6.05—10 8.65—10 
I.2. 05 05 05 605 605 
z.g. 60S 05 | -05 0s 04 
1.6. 05 05 | .65 65 61 
1.8 65 .65 65 65 57 
2.0 05 54 5° 52 49 
s.38 65 .40 42 45 °) 
2.4 65 33 25 26 .25 
2.6 a 10 | 00 7-95 03 
A 41 7.93 7.67 51 7.93 
3.0 21 04 25 0.95 32 
3.2 | 7-97 | 31 | 0.75 27 6.385 
Le ' 60 6.92 | IS 5.40 31 
3.0.. +] 35 45 5.48 4.53 5.70 
es +4:9-6 9.5 6 97 5.98 4.73 3.48 Ol 
4.0.. | 54 44 } 3.55 2.30 4.24 
eo | 07 4.84 2.95 00 3.40 
4.4.. 5.54 .Ig | 1.93 2.45 
4.6... onl 4.97 | 3.49 0.84 1.49 
cz | "2 ? 
oS) See 2.35 | 2.73 | 0.42 
a ‘f 3.08 Sine Aer 











‘ Differing from the corresponding value of p in Contribution No. 188, p. 15, by 


I 
the amount Mod log (0.0451) = 3.099. 
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modified definitions, that is, with (4) instead of with the values 
previously given for , gq, r. 
The values in Table III lead to the formulae in Table IV. These 
have been tabulated in Table Id. 
With these details premised we will now discuss the several 
points separately. 
TABLE III 








| °°? 30 60° } 90° 40° 90° 

ue 
kh... ...0..4-| —8.929 | —8.292 | —12.320 —17.418 | —10.984 
Sea | -+5.705 +5.481 + 9.320 +14.092 + 8.120 
err eee —1.366 —1.508 — 2.441 | — 3.542 — 2.171 








TABLE IV 


Gal. Lat. Density and Distance 
a eile a ooo log A(p) = —3.878+ 2.478 log p—o. 593 (log p)? 
Dates eG santa eas = — 3.60242. 381 —0.655 
Be adhe arn Aetoale wees = —5.351+4.048 —1.060 
eas sansa aaa aha cts wane =—7.565+6.120 —1.538 
oats nie a bia lana ie bard = —4.770+3. 526 —0.943 
POINT I 


As already mentioned, the solution of the problem involved in 
this case is contained in formula (13) of Contribution No. 188. 
Modified in accordance with what precedes, this formula becomes: 


* 00 


Nim = 3.046 p’A(p)o(m—5 log p)dp. (7) 


° 


From this equation it is at once clear that given V,, and (1) for all 
values of m and M, we can find A(p), and this constitutes in reality 
the solution of our problem. The derivation involves the solution 
of the integral equation (7). For the case in which ¢(/) has really 
the form (1) and N,, the form (2), Schwarzschild has given an 
elegant solution, which, quoted already in Contribution No. 188, 


t For this formula, in a slightly different form, see Astronomical Journal, 24, 27 


(No. 566), 1904. 
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but modified in accordance with our present definitions, is contained 
in the following formulae: 





A(p) dais ehtk log p+/ (log p)? (8) 
— 256" ) 
 7—C 
G?=—l—25r 
par O-D 
5r 
vuamsateel . (9) 
- rm 
logG RB 
h=a—p+—2— ——~ 2.520 
p . 4G? 5203 
=e oe ‘“ 
u Mod 6.9078. } 


Admitting that (8), which for brevity we will call the Schwarz- 
schild formula, holds for all distances from p=o to p=, it is, of 
course, easy to determine all such quantities as those which follow 
(the area is always assumed to be 10,000 square degrees) : 


N»= total number of stars of appt. mag. m 
(N,)? =the same within the distance p 
N =total number of stars of all magnitudes together 
N,=the same within the distance p 
M=median = distance within which lie one-half of all the stars 
W =(Nm)o/Nin 
Z=N,/N. 


We find the formulae: 


Ne =~: soon! ™ G(e+7) 1? a eotn (10) 
MU 
(Ns =WNe (11) 
gt 7GMog o—8) 
W=— { e~*dz (12) 
VF —o 
28G?= —R—10rm (13) 


yG?=h+ p+qm+rm? (14) 
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12.433 (e+5)’ 2 

, . - 

N= 29 gh— - ($=6.9078). (15) 
/ =f u 


As a check on the computations of /, k, 1, N we have, expressed 
in terms of a, b,c, 


- } ba 
N => ; A . Cn x (16) 
a 
k+3 ; 
—- 1 (3=6.0078) (17) 
N,=ZN (18) 
a I V —/(log p—log M) 
ae” { e~*dz. (19) 
7 J—« 


Introducing the constants of Table III into (12), (13), and (14), 
we obtain 


Gal. lat. =o° Gal. lat.=40° to go° 
G=1.833 G=2.039 (20) 
B=1.560+0.1182m B=1.550+0.0956m 


with the aid of which equation (12) gives at once an insight into 
the distribution in distance of the stars of each magnitude m. The 
results have been summarized for galactic latitudes o° and for 
40° togo’ in Table II. For the medians we find 
Gal. lat. 0° M = 39,000 peer 


(21) 


40° to go” 2,900 parsecs | 


Since the Schwarzschild formula does not hold for the smaller 
distances, results for these have been omitted from the table. 

But the supposition that we know ¢(M) and N,, completely, 
that is, for all values of M and m, is far from being justified. In 
the present state of science this incompleteness of our knowledge is 
very serious in the case of V,,, but of little account in the case of 
¢(M). The latter part of this statement constitutes the second 
point. 
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POINT II 
The practical sufficiency of our knowledge of the luminosity- 
curve is due to the fact, already alluded to, that it has now been 
found for so large a range of magnitudes, viz., from the very 
brightest down to about absolute magnitude M =-+o9. 
According to (6) a star of apparent magnitude m will be of 
absolute magnitude +9, when its distance is such that 


log p=? ° 


If, therefore, we write the equation (7) in the form 


m~9o 


New log p=—— 
3.046 -{ > p?A(p)o(m— 5 log p)dp+ 


ogp=—-x 


log p=+00 
{ fas p’?A(p)o(m—5 logp)dp (22) 
) 
5 


og p= 


the difficulty introduced by incomplete knowledge of the luminosity- 
curve lies exclusively in the first integral, for the second includes 
only stars brighter than M =+o, for which ¢(M) is known. The 
difficulty cannot, therefore, be very serious because the first integral 
involves only stars relatively near the sun. In fact, as will be 
explained below, the values of A(p) within the limits of the first 
integral may be considered as completely known and its divergence 
from (8), though considerable, still unimportant for the present 
question, by a method wholly independent of (7) for all apparent 
magnitudes to m= 23 (corresponding to the limit 2.8 for log p). 
It will undoubtedly be a long time before we can expect useful data 
for V,, beyond m= 23. 

This being so, it is easy to show that there is an overwhelming 
probability that the first integral in (22) is entirely negligible, at 
least to the limit m=23. We first determine the value of this 
integral, denoted by V,,, on the supposition that for M>+o, 
@(M) is still represented by (1). It is evident that V,, represents 
the number of stars of apparent magnitude m, absolutely fainter 
than M=-+o. : . 








to 


vi 
tN 
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We easily find 


J - a I “P _* 
” ig e~*dz (23) 
‘Vm Vr : 


pao —3(b—g) —cm | (24) 

The constants /, g, r are given by (4); a, 6, ¢ in Contribution 
No. 188, p. 13. 

Table V, computed by the aid of (23), shows that, if the 
luminosity-curve (1) were correct over its whole length, the stars 
fainter than absolute magnitude +9 would not contribute to the 
total number NV,, in the Milky Way more than 1 star in 16,000, 
or to those in galactic latitudes 40° to go° more than 1 star in 300. 
For N.,, N.;, etc., the contribution is much smaller still. 


TABLE V 


VALUES OF V,,/N,, 








m Gal. Lat. o° Gal. Lat. 40°-90 
OY oust suis wales ©.00000050 ©.QOOO0I 13 
| ey te er eee .COOOOT 43 . 0000414 
MR iuie os cele ate eis ; . 0000039 2 .0001 39 
asad ewsateG Metres . OOOOTO4 00043 
Oe ee eee .0000262 .OO121 
DS als Pld mcualerihien .0000634 0032 
ERR Oe as 000148 0076 
Ree ee a 00143 .0635 
eR attez tata Wiaturii Soacaniner coe a 0.0095 . 267 





Therefore, even if beyond M=-+o9, the error-curve (1) should 
begin to fail very seriously to represent the true frequencies of the 
absolute magnitudes, it will still lead to values of V,, which will 
be practically correct to magnitude 23 and probably to even much 
fainter stars. 

Suppose for instance that the frequencies, which up to M@=+9 
are represented with such astonishing closeness by the error-curve 
(1), should, beyond this magnitude, suddenly become tenfold the 
frequencies furnished by (1); even then the representation would 
be quite satisfactory for all values of V,, to m=23. Of course the 
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overwhelming probability is that matters are really much more 
favorable. For we must surely expect that even if the true curve 
begins to diverge from (1), it will not suddenly show tenfold 
values. We should rather expect a deviation which for M = +10 is 
still small and only gradually increases, perhaps to very consider- 
able amounts, as we approach fainter and fainter stars. If so, the 
divergence in the representation of N,, will undoubtedly be still 
much smaller." 

The conclusion from what precedes can only be: as long as obser- 
vations do not extend beyond m=23, formula (7), in which 
¢(M) is given by (1), must be considered as valid. Should they 
ever reach stars fainter than m=23, a direct extension of the 
luminosity-curve beyond 1J=-+9 may seem desirable. It is easy 
to see how this could be accomplished. We should require, how- 
ever, more extensive data on the very faint stars having very 
considerable proper motions. 


POINT III 


It now remains to be seen what result for A(p) can be deduced 
from formula (7) with the aid of our incomplete data for the values of 
N,,. First consider the case: 

a) We assume the form (2) to hold for all magnitudes from the 
very brightest down to m=, referring for less hypothetical con- 
siderations to III and IIIc. The problem evidently is: given, 
o(M) and N,, by the formulae (1) and (2) for all values of M and m; 
required, A (p) for all values of p. 

The solution has been given above in equations (8) and (g). 
Although theoretically this solution leaves nothing to be desired, 
we find in most cases that, for moderately large values of p, it 
becomes practically illusory. The reason for this of course is 
that existing data do not furnish the values of the constants h, k, 1, 
with sufficient precision. 

t Perhaps this will not be granted, for of course it is not impossible, mathematically, 
that after 14=-+0, the frequencies instead of continuing to decrease, may begin to 
increase at such a tremendous rate that the preceding conclusion would no longer hold. 
Even this possibility must probably be denied. If after M@=-+9 the frequencies 
should increase at such a rate, the number of very faintly luminous stars near the 


sun, that is, apparently faint stars with very large proper motion would, I think, 
exceed what is reconcilable with the observations. 
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An example, which has a considerable interest of its own, will 
help to show the difficulty in its true light. As such, we will try 
to derive the arrangement in space of the stars in the star-cloud 
between 6 and y Cygni. We choose this example because we 
happen to possess data for part of this region down to photo- 
graphic magnitude 18.5 (=visual magnitude 17.33). The part in 
question is that covered by Selected Area (syst. plan) No. 64. 
The low limit of magnitude is reached by using unpublished results 
of the Mount Wilson photometric survey together with those in 
the Harvard Durchmusterung.’' Simple countings gave the num- 
bers in Tables VI and VII. 


TABLE VI 


HARVARD DURCHMUSTERUNG 









































Harv. PHoToGRAPHIC Mac. N’ For 

— ——— ‘ m |} COMP. BY Oo-—( 

: 40’ X40’ | (25) | 

Limits Mean 
8.00 to BE  i8. 5 8e wed 8.5 I 2 — | 
9.00 to Saree gious 9.5 7 7 re} 
16.0000. 10:60.....+.<.. 10.5 21 19 | + 2 
11.00 to Rp eraree .3 42 | 49 — 7 
12.00 to 5 Ra eetog? 12.5 143 11g | T2 
ce) ies See 14.5 227 | 269 —42 
14.00 00 4-14.00......... 14.5 | 591 570 +21 
siesta ee ae Se a — 
| 
WE ret cae is Fa csletes aha | 1032 Jesse cece es 
TABLE VII 
Mount WILSON PHOTOMETRIC SURVEY 
M.W. PHotocRrApHic Mac. N’ vor N’ Repucep To 
m m 
Fal | | 15’X15’ 40’ X40 
Limits Mean 
| 
— a : oasieniLigieinadaini etna 
es eB er eee ea 13.0 | 203 146 
oe eee 14.0 433 309 
eS rr 15.0 | 973 693 
15.5 to 16.5 16.0 2003 1426 
ee en 17.0 3633 2584 
et 8 ere 18.0 6773 4817 
BR ee eee Fra oe Gtaeaantacaserae ore oe Eee 





' Harvard Annals tot. 














ARRANGEMENT OF STARS IN SPACE 255 


These tables do not contain even approximately all the stars 
in the two catalogues. The faintest stars measured on the Harvard 
plates are 15.75 (phot.), those on the Mount Wilson plates 18.86 
(phot.). Stars fainter than 14.99 and 18.50, respectively, were 
excluded because we cannot be sure that the catalogues are abso- 
lutely complete below these limits. 

We have first of all to reduce the magnitudes to the same scale 
and to the same area. The latter reduction has already been made 
in the last column of Table VII. To find the scale difference the 
observed numbers Nj, in Table VI were represented by an equation 
of the form of (2). The result may be written 


log Nyy = 2.074+0.370(m—12.5)—0.0145(m—12.5)?. (25) 


The values computed by means of this formula are in the fourth 
column of Table VI; the fifth shows the residuals O—C, whose 
changes of sign are all that can be desired. With the aid of (25) 
the Harvard magnitudes were calculated for which the numbers 
Nj, become equal to the first three values of Vj, in Table VII, viz., 
146, 309, 693. The results are 12.75, 13.68, and 14.78, whence we 
conclude: 


Mt. W.—Harv. 
Harv. mag. 12.75= Mt. W. mag. 13.0 +0725 
Harv. mag. 13.68= Mt. W. mag. 14.0 +0. 32 
Harv. mag. 14.78= Mt. W. mag. 15.0 +0. 22 


There is no indication in the differences Mt.W.—Harv. of change 
with the magnitude, and we adopt 


mag. Mt. W.—mag. Harv. =Const.=+0™26. (26) 


By means of (26) the values in Table VII can now be reduced 
to the scale of Table VI. , 

A further reduction from the photographic to the visual scale 
is necessary because the luminosity-curve defined by (1) and (4) 
gives the frequencies for the visual magnitudes. For this we 
have used the corrections given in Groningen Publications, No. 27, 
p. 42. The first five columns of Table VIII show the complete 
reduction. ; 
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» VIII 


TABLE 


FOR 40’X 40’ 


Ni, 


COMPUTED 


ALI 


Ns 


PHOTOGRAPHI 


AmPI 


HARV 


SCALE 
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The fifth and sixth columns give the observational data fully 
reduced for the present discussion. The numbers in the cloud 
were first compared with the corresponding numbers for the average 
Milky Way. It turned out at once that these numbers are all 
but perfectly proportional, the cloud having about 2.39 times the 
number for the average Milky Way. This is shown by the seventh 
column of the table which gives the carefully interpolated values 
for the Milky Way according to Groningen Publication, No. 27, 
reduced to an area of 40’X4o’ and multiplied by 2.39. The 
agreement with the observed cloud numbers in the sixth column is 
as good as can be expected. The numbers in the Milky Way in 
Groningen Publication, No. 27, do not extend beyond magnitude 16. 
For these very faint magnitudes the table gives in parentheses the 
values furnished by 


log Nin = —4.516-+0.7242m—0.0141m? (27) 


which best fits the numbers of Groningen Publication, No. 27, 
reduced to 40’ X40’ and multiplied by 2.39. 

On the assumption made in this paragraph, that the formula 
holds for all magnitudes, we conclude that at all distances the 
star-densities in the cloud are 2.39 times the corresponding densities 
in the average Milky Way. 

Meanwhile we can represent the observed numbers quite as 
well by a formula somewhat different from (27); for instance, by 


log Ni,= —3.951+0.6334m—0.0106m?, (28) 


The representation by this formula is also shown in Table VIII. 

By means of the constants in (27) and (28) and those of the 
luminosity-curve (5), we can,at once pass from the distribution of 
the values N,, to that of the densities by means of formulae (8) 
and (g). We thus obtain 


Solution (27 Solution (28) 
ee Pee — 8.060 — 3.488 
ee +2.678 (29) 
Bilis asomaaieies —1.375 —o. 882 


The two solutions represent rather widely different arrange- 
ments in space. The best way of showing the difference in 
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compact form will perhaps be to compute for both, by formulae (16) 
and (17), 
N =total number of stars for a field of 40’ X40’, and 
M = median distance below which lie just half of all the stars. 
Introducing the values (29), we find, 


Solution (27) Solution (28 
Be ptalala ik hicsaeiin 39,000 parsecs 272,000 parsecs 
ah (30) 
chica esis 579,000 stars 3,700,000 stars 


Identical data thus lead to two solutions for the arrangement in 
space of the stars in the cloud. In the second, judging by the 
median values, the star-distances are 7 times those of the first 
solution. The total numbers’ are as 1 to 6.5. 

Evidently, therefore, the data at present available do not lead 
to a reliable determination of the real structure of the cloud. The 
reason clearly lies in the fact that the data of Table VIII are insuff- 
cient for a good determination of the constants for V,, in equation 
(2), and it is readily understood that it is mainly the uncertainty 
in the position of the maximum to which the failure of the solution 


is due. 
In the notation of (2), this maximum lies at the magnitude 
b , 
mM. = — (31) 
2¢ 
whence for 
Solution (27) M,= 25.68 
—P : (32 
Solution (28) My = 29.88 


a difference of over 4 magnitudes! 

We are thus led to consider the question as to how far the 
observed data must be extended in order to obtain a determination 
of A(p) which will satisfy moderate demands at least. 

To answer this it will be necessary to determine: (a) how the 
precision of m, (prob. err. p»,) depends on the precision of the 


1 Since the total number of stars resulting from the two solutions differs enor- 
mously, it might seem worth trying to decide between them by measuring the total 
light of the cloud. To test the feasibility of such a plan I compared the values of the 
total light and found that the amounts involved in (27) and (28) are as 1.000 to 1.025. 
Practically, therefore, the idea is worthless, at least with the photometric means at 
present available. 
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counts V,,; (6) how the precision of the median M depends on the 
precision of mp. 

As for (a): Consider first the somewhat easier question: What 
will be the accuracy of m, obtained from observed values of Vp, 
for three, favorably chosen, equidistant magnitudes 


m,—o; My; m,+a? (33) 


For brevity call this the simplified problem. Since a rough esti- 
mate is all that is wanted, assume that the probable error of log Vn 
for these three magnitudes is constant and equal tor. When this 
simplified problem has been solved it will be easy to estimate 
roughly the accuracy obtainable by the use of all values of V,,. 

In the Appendix is derived the formula which shows the relation 
between the probable errors p,,, and r, namely, 


Pm = = a V 12(m,—m,)?+07* (34) 
40°C 

The value of r is readily determined by means of the deviations of 
the observed log N,, from a close-fitting curve. In the case of the 
star-cloud the number of observed stars brighter than 12 is rather 
small. This being the case, consider only the Mount Wilson 
observations in Table VIII, which show a better agreement inter se 
than the Harvard values. For the best-fitting curve I used a 
third solution, intermediate between (27) and (28), for which the 
constants are: 


a=—I0.019 h=—6.038 

b=+1.570 k=+4. 289 (area = 40’ X 40’) (35) 
c= —0.0204 l= —1.165 
My = 26.73 Median M=63,800 parsecs 


The value found for the probable error is r= +o0.008. Experi- 
ence gained in other cases indicates that this value is accidentally 
low and that a better one would be 

r=+0.012. (36) 

For the fundamental magnitudes we choose the values (see 


Table VIII) 
12.18, 14.52, 16.86, (37) 
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then for the computation of p», by (34 


m,=14.52 c¢= , 
I 4-5 34 (28) 
M, = 26.73 C= —0.0204 
which lead to 
Pm, = = 2.0 mag. (29) 


Starting from solution (27) instead of (35), we find nearly the 
same result. 

As for (6): We learn how the accuracy of the median M depends 
on that of m, by comparing the values of m, and M in the three 
solutions (27), (28), (35). Adding a fourth solution, we have 


altogether 
Mo log M 
25.68 4.592 
I.05 0.213 
26.73 4.505 (40) 
I.05 0.203 
27.78 5.008 
2.10 0.420 
29.88 5-434 


Evidently log M changes very nearly proportionally to m,. In 
fact 

log M=—0.544+0.2 mo. (41) 
Therefore 

P log M=O.-2 Pm; (42) 


consequently, adopting the values (39), 


P log M= 0.520 (simplified problem). (43) 


It is now easy to see what can be gained by future extension of the 
counts to fainter magnitudes. 

Table IX shows how the probable errors of both m, and log M 
change with increase in the limiting magnitude. In accordance 
with what precedes, we adopted + 2.6 as the probable error of m, in 
the case of the limiting magnitude 17.5, i.e., for limiting V,,=17.0. 
The remaining probable errors were then computed by formulae 
(34) and (42). The results are summarized in Table IX. 
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Hence, if we were able to photograph the stars to magnitude 
20.5 (visual scale), the probable error of log M in the simplified 
problem would be diminished to 0.176, which, roughly speaking, 
is equivalent to a probable error of about 42 per cent in the numeri- 
cal value of M. I believe the error would be further lowered to 
less than 30 per cent, if the calculation of m, were based on all the 
available data for ,, instead of merely on those for the three 
fundamental magnitudes used for the simplified problem. 

From an absolute point of view this precision may not seem to 
be very satisfactory. Still if it can be reached—and at present 
this is hardly subject to doubt—an enormous step in advance 
will have been made, for it means that the median parallax will 
have been determined with a probable error of about +0”oooor. 


TABLE IX 


PROBABLE ERROR OF mt, AND LOG M (SIMPLIFIED PROBLEM) 





onto Fund. Mags. We - Sol. (3 5) Mo Plog M 
iS Sere 2.0 %4.§ t7.0 14.5 ee 2.00 0.520 
eee ee 12.0 15.0 18.0 15.0 3.0 1.73 346 
19.5. r2.0 15.5 19.0 IS.5 2.5 aad jt.ae 7 242 
a ae I2.0 16.0 20.0 16.0 4.0 ia 0.88 176 
ee, ee 12.0 16.5 21.0 16.5 4.3 0.66 132 
ee 12.0 17.0 22.0 17.0 5.0} 0.51 | 0.102 





Moreover, the accuracy can be greatly increased by taking several 
fields instead of but one. 

The preceding considerations on the structure of one small 
region may of course be applied without change to the investigation 
of larger parts of the sky, for instance, to the Selected Areas for 
different zones of galactic latitude. The only difference is that 
we will obtain a much greater accuracy. 

From the data of Groningen Publication, No. 27, Table V, we 
have found the values of m, for different galactic latitudes (8) 
entered in the second column of Table X. The other quantities 
in the table were obtained in the manner already explained. 

In this table the probable errors are given, in the simplified 
problem, as they would be found for a single region containing the 
same number of stars as that in the ease treated above. In con- 
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formity with (39), p», for 8=o and lim. mag. 17.5 was put at + 2.6. 
For the rest the computation was made with the aid of (34) and (42). 

As was to be expected, the values of m, run somewhat irregularly. 
We therefore adopted the smoothed values of the third column, 
which, with the exception of that for 8 = 20°, represent the observa- 
tions practically as well as those of the second column. 

By using for each region all the material instead of the three 
values of N,, of the simplified problem, the probable error will 
be diminished. If further for each galactic latitude we take the 
greater part or all of the Selected Areas, the probable error will 
again be much reduced. Finally, a further important reduction 





TABLE X 
PROBABLE ERROR IN SIMPLIFIED PROBLEM FOR A REGION 15/X15' 
| | ; 
| | Lim. Mag. 17.5 Lim. Mag. 20.5 
Adopted | ——_—_— — | ———— — 
8 | Mo mMo | . | | 

| | Pme Plog M | Pome Plog M 
eee | 26.7 26.7. |—0.0294 | +2.6 +0.52 +o.88 | +0.176 
20. 23.28 25.6 |— .0297 2.4 .48 | 79 1 158 
40. | 23-1 22.9 [- .O333 1.6 | 32 | 51 .102 
60 19.1 19.9 |— .0407 0.85 | 7 | 24 048 
80 | 1.2 | 17.9 |—0.0486 | #0.45 +0.09 +0O.11 +0.022 








can be made by taking into account the very extensive data for the 
brighter stars. In order to see more clearly what can be gained, 
some estimate must be made as to the effect of all these reducing 
causes. The matter certainly deserves careful study. We have not 
yet made such a study, but estimate provisionally that it should 
be easy to reduce the probable errors of Table X to one-third.’ 
For clearness we may express roughly the accuracy correspond- 
ing to this assumption in percentages of the whole. The results 
for the median M are thus found to be as given in Table XI. 
Granting the suppositions made in IITa, and that we are in 
possession of data for V,, to magnitude 20.5, we conclude that 
we can obtain what may be considered a good insight into the 
t The comparison of the values m, in the second column of Table X with the 


smoothed values in the third leads at once to a precision of about this amount, not- 
withstanding the fact that the data of Groningen Publication, No. 27, extend but little 


beyond magnitude 14.5 visual. 











ARRANGEMENT OF STARS IN SPACE 263 


distribution of the stars in all galactic latitudes. With data to 
17.5 alone (Mount Wilson Survey) this can be said only for the 
regions of high galactic latitude. Below latitude 40° the solution 
becomes rapidly less reliable. ; 
Remark.—Meanwhile we need an extension of star-counts for 
regions near the galactic pole quite as much and probably more 
than for the lower latitudes, for just these observations of W,, 
in the higher latitudes will lead to a judgment on the general 
validity or non-validity of the analytical form (2). This follows 
because of the lower value of m, (Table X). Suppose for instance 
TABLE XI 


PROBABLE ERROR OF MEDIAN M IN PERCENTAGE 
oF Its AMOUNT 














! 
Limiting Magnitude 
| . 





Gal. Lat —- 
17.5 20.5 
Bes odie Srsiaituiw boone 41 
20 | 38 2 
ee ee ee ee 25 8 
bi othe de we aeaane | 13 4 
ir aik nad Oona ek 7 ‘5 





that we obtain counts to apparent magnitude 20.5. In high lati- 
tudes the curve for N,, will then become known considerably 
beyond the maximum, that is, for a very considerable fraction of the 
whole, and a fraction very much greater than that covered in the 
lower latitudes. In case formula (2) is confirmed for this large 
part of the curve in the region of the galactic pole, we will use 
the same analytical form with much greater confidence for the 
lower latitudes also. 

IIIb. We know by observation that as far as magnitude 14 the 
values of NV, are well represented by (2); for the stars fainter than 
14 we know nothing of the sort... What results will it be possible 
to obtain from this knowledge ? 

Schwarzschild’s solution, which is based on the supposition that 
all values of N,,, fall on the curve (2), led to formula (8) which, but 
for the observation errors, would yield a knowledge of A(p) for all 


‘We take the limit m=14 because it represents what at present is really well 
known. The reasoning, of course, holds for any other limit. 
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distances. This solution is not now generally applicable and we 
must look for another. Such a solution, based on quite different 
principles, which I will call solution A, has already been given in 
Groningen Publication, No. 11, and has also been applied in Contribu- 
tion No. 188. It has there been carried through for galactic 
latitude 40° to go, and shows the very gratifying fact, which 
we had no right to expect beforehand, that, but for small distances 
(where Schwarzschild’s solution is confessedly wrong) it agrees 
wonderfully well with the Schwarzschild solution. With the excep- 
tion referred to, therefore, solution (7) represents all the data, those 
of Contribution No. 88 as well as the values of V,, in a way that 
leaves hardly anything to be desired. 

Therefore, notwithstanding what has just been said, the 
Schwarzschild solution (8) must still be considered the most 
plausible to be obtained with the aid of the data at present avail- 
able. This must be owing to the fact that further values of V,,, 
particularly those for the most influential magnitudes 15, 16 

... , do not deviate very much from the curve (2). Still 
this does not prove that formula (8) gives even an approximate 
solution for A(p) at very great distances. We see, for instance, from 
Table II, that according to (8) there are, in high galactic latitudes, 
hardly any stars beyond log p=4.2. As a consequence thereof, 
suppose that between log p=4.0 and log p=4.2, the density is 
a hundred fold the value given by (8); the total numbers of 
stars N,, will be changed by about 1 per cent only. For the 
brighter magnitudes the change would be still smaller. Evi- 
dently, therefore, nothing can be concluded from our data about 
the densities at such great distances. They may or, may not 
agree with formula (8); we have no means of judging. The 
matter would be different if we had good counts (N,,) for very 
much fainter stars. The question thus presents itself: To what 
distance can the densities computed by (8) be accepted with some 
confidence? Or to be more precise: To what distance can we trust 
the values of A(p), found by (8), to be correct within 25 and 4o 
per cent, respectively, of their values ? 

In considering this question we will assume that what was found 
in Contribution No. 188 for galactic latitudes 40° to go° will also be 











ARRANGEMENT OF STARS IN SPACE 265 


found for other galactic latitudes, viz., complete agreement of 
solution K with solution (8). It seems to me very probable that 
such will turn out to be the case. As solution K is independent of 
the values of N,, for m>14 we may thus accept formula (8) for 
log pS2.8' (excepting of course the small distances). 

The question under consideration comes to this: How much 
can we change the densities (8) for log p>2.8, without spoiling 
the representation of the values of V,, to m=14? Suppose that, 
for galactic latitudes 40° to 90°, we multiply the densities for 


log p=2.8 to 3.0 by F,, 


log p=3.0 to 3.2 by F.,, etc. 


Then the numbers in Table II (second part) for log p=2.8 to 3.0 
will all be multiplied by F,; those of the following line by F.,, etc. 
The total number Vs (see Table IT) will thus be changed to 


N3(0.9190+0.0569F',+0.0188F,+0.0044F,+0.0008F,+0.0001F) (44) 


and, in order that the Vs may still be represented, the sum must still 
equal Ns; or better, the expression in parentheses must equal 1.00. 
We also find similar expressions for magnitudes 9, 10, . . . . 14. 

The question thus becomes: What is the smallest distance for 
which we can take for any one of the factors F, a value as high as 
1.25 or as low as 0.75 (or, respectively, 1.40 and 0.60) without mak- 
ing it impossible to determine the remaining factors in such a way 
that the seven equations of condition are still bearably well satisfied. 
For “‘bearably well satisfied’’ we have used the condition that the 
range in the residuals should not exceed to per cent. 

Suppose, therefore, having assumed the constant F for distance 
R to be 1.25, that we solve our equations of condition and obtain 
‘bearable”’ set of residuals. Will this prove that for dis- 
tance R a density deviating 25 per cent from (8) is admissible? By 


‘ 


indeed a 
no means. For as a rule the values found for the other constants 
will be impossible or inadmissible because they lead to negative 


values of A(p) or to values running too wildly for the regularly 
increasing distances. 


In Contribution No. 188 this is the limit to which solution K extends. 
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In order to obtain an acceptable solution we must evidently 
introduce the condition that the values of F shall run smoothly for 
regularly increasing distances. This condition is pretty vague 
and introduces a certain degree of arbitrariness, but I think that 
there is really no serious objection to it. Moderate differences in 
the form of the curve will probably change the results little as long 
as we leave two or three constants indeterminate, to be found in a 
way best satisfying the conditions of the problem. 

Before stating the form adopted for the factors F, consider the 
general course they must follow. From the value 1.00 for log p = 2.8 
they must either increase or decrease. Consider only the first 
case, because by satisfying this we more than satisfy the other. 
The initial increase cannot indefinitely continue because if all the 
factors are greater than 1.00, the total numbers JV,, will all be 
increased. The values of N,, cannot well be represented, therefore, 
unless values of # >1.00 are compensated by later factors< 1.00. 
The factors must therefore pass through a maximum, after which 
decrease sets in and they soon become 1.00. We assume that 
there is but one maximum. Consider now the various curves rep- 
resenting various hypotheses for the factors F. Take more curves 
for which the value of the maximum is respectively Fyg:=1.25, 
Fmaz= 1.30, Fmoz=1.35, .... and let R, R,, RK, .... be the 
values of the distance for which F=1.25. 

Then we think that generally 


4 1 oe 


That is, the factor F=1.25 will be reached sooner by that curve for 
which F=1.25 is a maximum.’ 

Since we wish to find the smallest distance at which F can become 
equal to 1.25, we consider only the latter case, that is, find the 
minimum distance, R, at which F can reach 1.25 as a maximal value, 
without being led to an “unbearable”’ representation of V,,. 

This being premised, we now give the mathematical form 


adopted for F: 


* We have not tried to find the demonstration. We have been content to verify 
its truth further on by means of the data in Table X. 
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A. For the case, maximum factor=1. 25 


Before the maximum, F =0.96+0. 29 e~*”* } (45) 


) 
| 
After the maximum, F=1.25 e~?*” | 


\ 


B. For the case, maximum factor =1.40 
Before the maximum, F =0.96+0.44 e~”*" ; (46) 


> ” 22 
After the maximum, Ff =1.40 e~?* } 


in both of which 


«=log p—log R. (47) 


There remain, therefore, the three constants R, h, p, to be 
determined in such a way that R corresponds to the smallest distance 
for which the values of N,, can be represented with residuals not 
exceeding 10 per cent. This form for the factors F gives perfect 
smoothness throughout. It is that which a priori seems most 
plausible. Beyond the maximum, A(p) is evidently again repre- 
sented by a form like (8), which represented the densities for the 
smaller distances with such perfection. 

It does not, however, start at log p=2.8 with the initial value 
1.00, as premised. The reason is that the small discontinuity thus 
introduced at log p= 2.8 is of no real importance, while the repre- 
sentation of N,, is appreciably better. 

Substituting the expressions for F into (44) and equating the 
result to the successive values of V,, we obtained a number of equa- 
tions of condition. 

In solving them we generally began by assuming two values for 
R, for each of which these equations were solved for # and /p in 
such a way that practically the best possible representation of 
N,,(to m=14) was obtained. The range of the residuals was then 
found for both cases. The value of R for which the range is 10 
per cent could then usually be found by interpolation. If not, the 
computation was repeated for a third value of R. In this way 
were found all the values entered in Table XII. 

In all cases, with the exception of the third for m=14, the 
range of the residuals in N,, is nearly 10 per cent. In the excep- 
tional case it is only 6.5 per cent. Infact, in this case, the knowl- 
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edge of V,, to m=14 takes us little beyond what was reached 
by solution K. 

From Table XII we may be pretty sure, now that we know 
NV, with some precision to m=14, that in the Milky Way there will 
be in the densities as computed by (8), no error of 25 per cent 
up to log p=3.1 and no error of 40 per cent up to log p=3.2. As 
soon as we have reliable knowledge of N,, to m=20 these limits 
will rise to log p=3.8 and 4.0, respectively. We further learn 
from the table that by extending our knowledge of N,, from m=14 


rABLE XII 


Limits OF RELIABILITY OF A(p) BY (8 











Di Oe N ” V ag K N m = 20 
Gat. Lat os GAL. Lat —— ——- 
Fmax r R } p I max rR h ? 
‘e) 1.25 3.2 7.0 1.00 O 1.25 3.0 4.0 1.00 
° are I.40 a 0.0 70 ° I.40 4.0 BaF 3.0 
40 to go...}] 1.25 s.9 2.20 40 tog [25 3.4 4.0 2.05 
40 to go 1.40 2.9 3.0 40 to go I.40 ; 4.0 ee ke 


oon 


to m= 20 we will extend the limits to which we can reliably determine 
the structure of the system: 


JX 


In the Milky Way 5.6 fold (log=o.7 


, (48) 
In Gal. Lat. 40° to go° 2.5 fold (log=o. 5s + 


In the case of the Mount Wilson Survey, which will yield V,, 
to m=17, the gain will be about half this amount. 

IIIc. If finally, having extended the counts NV,, to some higher 
magnitude, we find that the values of N,, no longer fit an error- 
curve like (2), then evidently formula (8) will no longer hold for 
the densities. In this case we must fall back on formula (7). 
We will thus find different values for the densities A(p). We may 
certainly expect, however, that, as long as these values do not 
deviate very widely indeed from the densities furnished by such 
a formula as (8), our considerations on the limits of reliability will 
still practically hold. For the present, therefore, we may safely 
rest content with the results just obtained. 
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POINT IV. EXTENSION OF THE LIMITS OF TABLE XII 


There is, we think, reason to assume that the limits shown in 
Table XII give rather too unfavorable an impression of the real 
reliability. The reason is that the deviation, having once reached 
25 (or 40) per cent, will grow no further, or at least not much further 
before it begins to diminish. This it will continue to do until it 
becomes zero and begins to increase again and, possibly, go on 
increasing indefinitely (on the other side). That is to say, beyond 
the limits assigned in Table XII the deviation from (8) for a certain 
time will still be mostly very moderate or small. Still much more 
important than this is a fact which will appear in a further publica- 
tion, by one of us, in which a formula is derived theoretically which 
in the last analysis is an expression of the differential quotient 
dA(p)/dp as a function of p and A(p). The integration of this dif- 
ferential equation, if the equation could be accepted as definitely 
proved, would give A(p)as a function of p for all distances. As it is, 
the equation is verified by observation as far as observation can 
fairly be relied on. 

We thus have not only reliable values of A(p) up to the limits of 
Table XII, but also to these same limits, reliable values of dA(p) /dp. 
This proves that we have a right to extend the use of the formula 
which yielded A(p) considerably beyond the limits of Table XII. 

There is a third consideration which will enable us to extend 
the limits of Table XII with some probability of reliable results. 
In the publication just alluded to it will be shown that the equi- 
density surfaces are at least approximately similar ellipsoids. If 
we suppose that this form holds somewhat beyond the limits for 
which we have fair certainty, then the limits of Table XII will 
again be extended, at least for the lower galactic latitudes. 

The three considerations (of which the second is the most 
\important) together make us believe, that as soon as we have good 
data for V,, down to magnitude 17, which will be very shortly, we 
shall be able to find fairly reliable values of the A(p) for the entire 
region near the center for which A(p) exceeds one hundredth of 
density near sun, that is, for the whole of the space represented in 
Fig. 2, Contribution No. 188. 
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CONCLUDING REMARK 
The conclusions arrived at in this paper may perhaps best be 
seen in their true bearing by comparing the limits within which we 
can find fairly reliable results for the structure of the stellar system 
by using different sets of data. Of course Table XIII claims no other 
merit than that of furnishing the means of a rough insight into 
what can be obtained at present or in the very near future. 


TABLE XIII 


LIMITS IN PARSECS WITHIN WHICH THE STRUCTURE OF THE STELLAR 
SysTEM CAN BE Founp 


Method $$ —$—$——$_$__—_—_—__—_ 
o'-20 40° -QO Qo" 

Ce ee a - = = - —_ ee - _ — > 
1. Direct parallax determination... . oF = ae ere per 50 
2. Parallactic motion, now well known, to m=10 320 240 300 
ee - (830) | 610) | (720) 
4. Parallactic motion, stars to m=10 and up=0%0!1 400 | 320 | 360 

5. The same together with N,, (the latter to 

M=14)...... ey 1600 | 800 | 1200 
6. The same, V,, to m=17 Se re (4000) | (1600) (2800) 
7. Extension according to Point IV,V,, to m=14 3000—osiY| 1000 2000 

8. The same, N,, to m=17....... pawunnales (8000) 1700 A=A,/ 100 





NOTES TO TABLE XIII 

The limits in parentheses have not yet been reached but may be attained in two 
or fhree years. 

1. The magnificent results lately obtained in direct parallax determinations by 
means of long-focus instruments show probable errors of the order of o’c1. Even if 
there were demonstrably no systematic errors left, a probable error of o’01 would 
imply that we cannot derive from them the structure of the system beyond the limit 
of the table. Since the absence of systematic error cannot be claimed, the limit of 
the table is probably rather high. 

2, 4. The table gives simply the mean distances corresponding to the mean paral- 
laxes in Groningen Publication, No. 29, and in a paper not yet published. 

3. Simple repetition of the parallax plates obtained at the Mount Wilson Observa- 
tory some five or more years after their first exposure, would give excellent values for 
the parallactic motion of the stars of magnitude 13. Similar repetition of the parallax 
work done in America on Boss stars by other observatories with long-focus instruments 
would also contribute to an important degree. 

5,6. From Table XII of the present paper. 

7, 8. According to Point IV of present paper. In the last column of No. 8, 
A=A,/100 means that for all latitudes the densities are fairly reliable for the whole 
of the space where the density is greater than 0.01 of that near the sun (compare 
Fig. 2, Contribution No. 188). 
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All the results will be greatly strengthened as soon as we know 
N,, to m=20. At the same time this will enable us to make good 
estimates for still greater distances. 


APPENDIX 
POINT IIId. ON THE SIMPLIFIED PROBLEM 
Given, r= probable error of each of the quantities log N,, log N,, 
log N;, required, p»,=probable error of the constant m,. 
By formulae (31) and (2) of the text 
log Nin =a — 2umcm+ pcm? (u=mod.). (a) 
Hence 
log Ni =ya—2um,c(m,—o)+yc(m,—a)? 
log N,=ya—2yum,.cm,+ pom? (d) 
log VN;=ya—2umc(m.t+o)+yc(m,+o)? 


from which for m, 











My—m a» tt (c) 
° "T o(H—1) 
in which 
_ log N,—log N, 
a= log V,—log N,° (4) 
Consequently 
o 
Pm" (H—1p "8 (e) 
in which the probable error py must be found from (d). Since we 
assume 
Probable errors of V,, N,, and N, each=r, (f) 


we obtain as in the text 


rv 2 
Pm. = a V 12(m,—m,)?+o? . (g) 


LEIDEN, HOLLAND 
August 1921 














COMPARISON STARS FOR NOVA PERSEI, No. 2 
A CORRECTION 


Father Hagen has called my attention to an error in the identi- 
fication of one of the comparison stars for Nova Persei, No. 2, 
whose magnitudes are given in Mt. Wilson Contr., No. 192; Astro- 
physical Journal, 52, 183, 1920. ‘The star in question is No. 57 of 
his sequence. An uncertainty in the position and brightness was 
noted during the measurement of the photographs, but was over- 
looked when the MS was prepared several months later. 

The confusion arose from the fact that the declination co- 
ordinates of stars 57 and 62 (+12'6 and +11'5, respectively) are 
interchanged in the Chart and Catalogue of Father Hagen. The 
star actually observed here is close to the catalogue position of 
No. 57, but is much fainter than that object. Neither No. 57 nor 
62 has been measured on the Mount Wilson photographs; No. 57, 
in fact, is just outside the limit of distance from the center of the 
plate usually adopted. The .visual magnitudes on the Mount 
Wilson scale are readily found, however, by reducing Father 
Hagen’s values to that system. ‘The results are 


No. 57 a7: No. 62 128 


/ / 
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